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The RICIS Concept 


The University of Houston-Clear Lake established the Research Institute for 
Computing and Information Systems (RICIS) in 1986 to encourage the NASA 
Johnson Space Center (JSC) and local industry to actively support research 
in the computing and information sciences. As part of this endeavor, UHCL 
proposed a partnership with JSC to jointly define and manage an integrated 
program of research in advanced data processing technology needed for JSC’s 
main missions, including administrative, engineering and science responsi- 
bilities. JSC agreed and entered into a continuing cooperative agreement 
with Ui ICL beginning in May 1986, to jointly plan and execute such research 
through RICIS. Additionally, under Cooperative Agreement NCC 9-16, 
computing and educational facilities are shared by the two institutions to 
conduct the research. 

The UHCL/RICTS mission is to conduct, coordinate, and disseminate research 
and professional level education in computing and information systems to 
serve the needs of the government, industry, community and academia. 
RICIS combines resources of UHCL and its gateway affiliates to research and 
develop materials, prototypes and publications on topics of mutual interest 
to its sponsors and researchers. Within UHCL, the mission is being 
implemented through interdisciplinary involvement of faculty and students 
from each of the four schools: Business and Public Administration, Educa- 
tion, Human Sciences and Humanities, and Natural and Applied Sciences. 
RICIS also collaborates with industry in a companion program. This program 
is focused on serving the research and advanced development needs of 
industry. 

Moreover, UHCL established relationships with other universities and re- 
search organizations, having common research interests, to provide addi- 
tional sources of expertise to conduct needed research. For example, UHCL 
has entered into a special partnership with Texas A&M University to help 
oversee RICIS research and education programs, while other research 
organizations are involved via the “gateway” concept. 

A major role of RICIS then is to find the best match of sponsors, researchers 
and research objectives to advance knowledge in the computing and informa- 
tion sciences. RICIS, working jointly with its sponsors, advises on research 
needs, recommends principals for conducting the research, provides tech- 
nical and administrative support to coordinate the research and integrates 
technical results Into the goads of UHCL, NASA/JSC and industiy. 
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Executive Sunmary 

Obiectives: The objectives of this project are to develop a sy: stem for 

displaying computer graphics images of space objects and to demonstrate 
the use of this system as a testbed for evaluating vision systems for 
space applications. 

Overview: In order to evaluate vision systems, it is desirable to be able 
to con trol all factors Involved in creating the images used for processing 
by the vision system. Considerable time and expense is Involved in 
building accurate physical models of space objects. Also, precise 
location of the model relative to the viewer and accurate location of the 
liqht source require additional effort. Although the motion of a small 
satellite model can be controlled by a robotic manipulator, such as a 
PUMA this task requires additional equipment and time. On the otner 
hand’, the position and motion of a computer graphics generated model can 
be controlled accurately with precise light location. 

As part of this project, graphics models of space objects such as the 
Solarmax satellite are created such that the user can control the light 
direction and the relative position of the object and the viewer. The 
work is also aimed at providing control of hue, shading, noise and shadows 
for use in demonstrating and testing imaging processing techniques. The 
simulated camera data can provide XYZ coordinates, pitch, yaw and roll for 
the models. A physical model is also being used to provide comparison of 
camera images with the graphics Images. 
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INTRODUCTION 


A robotics software simulation testbed is being developed by NASA 
Johnson Space Center to identify the enabling automation and robotics 
technologies for space operations into the next century. The project 
described in this report is aimed at evaluating vision systems for space 
applications and providing assistance in the overall effort to develop the 
robotics simulation software testbed. The three complimentary efforts 
involved in this project are graphics modeling, physical modeling, and 
vision system evaluation. Research results for each of these efforts is 
described below. 

Initial work on the graphics modeling was done using APL GRAPHPAK 
software. This permitted rapid conversion of the dataset for a graphics 
model of the Solarmax satellite to be drawn using a dataset obtained from 
Johnson Space Center. A FORTRAN program was written to permit conversion 
of the dataset for use with a Sun worksta t ion working with a network file 
server with the mechanical engineering department Celerity C1200 
computer. Once data had been translated over in a form suitable for use 
with the Sun III workstation, the testbed was developed in a form 
containing menus that permit the user interactive control of the graphics 
modeling. Results using a simple physical model are compared with those 
obtained with the graphics model and both of these models are utilized in 
the evaluation of a computer vision system. The 3M VDL vision system is 
described and examples of its application in image processing are given. 


GRAPHICS MODELING USING APL GRAPHPAK 


One of the objectives of this research Ts to create a computer 
graphics testbed for simulating images of a slowly spinning satellite as 
viewed by an autonomous robot vehicle vision system. In this work it is 
necessary to control the relative location and attitude of a satellite 
with respect to the observer. 

In the preliminary work on this project APL was used to generate 
images of the Solarmax satellite. A brief discussion of image generation 
using APL is presented and results are given for several example positions 
of the satellite relative to the viewer. Computer programs, data and 
additional output images are collected in the Appendices. 


Image Generation Using APL 

Images were produced for different orientations of a satellite with 
hidden line removal. Rotation of the reference frame involves matrix 
manipulation that APL can handle concisely and powerfully. 

The two features of APL that are extremely useful are: 

1) the use of arrays to store all kinds of data 

2) the use of a large number of mathematical functions (each coded 
as a single symbol) to manipulate arrays. 


One main advantage of this approach is the great number of operations 
that can be performed with very short programs. A second advantage is the 
fact that most of the operations are very similar to those seen in 
mathematics. Unfortunately, to be able to use the language, the user must 
learn: 


1) the meaning of quite a few special symbols; 

2) ways to organize data into array format. 

GRAPHPAK is an APL workspace containing predefined functions that 
produce drawings on graphics devices. GRAPHPAK provides us with the 
ability to draw and transform projections of 3-dimensional objects. 

A data set was prepared for the Solarmax satellite and programs or 
functions were written for hidden-line removal to handle data transmission 
between TSO and APL environments. These functions and data sets are given 
in the Appendices. 

Three-dimensional objects are represented by a 4 column array, the 
last three are the x, y and z coordinates of a point and the first column 
contains either a zero or one. APL uses the line method to draw, so a 0 
means "move", and a 1 means "draw". The coordinates used are screen 
coordinates, with positive-x to the right, positive-y up, and positive-z 
out of the screen. 

The original data set "solarmax.vec" listed in the Appendix was 
modified to permit hidden-line removal. For convenience, we separated the 
data set into seven groups, which are ANT, BASE, BODY, HOLE, JET, SPAN and 
TORSO. Also the vertices for each surface must appear in a counter 
clockwise order looking toward the surface. Function "pickup" can create 
the required data from raw data by specifying the surface number only (see 
Appendix). 

Finally, we combined all the subpart data, already processed by the 
function "pickup", into one big data set " sol art". The data set rowse 
is created by the starting and ending vertex number of every surface. 
Both data sets are fed into function "hidden". 


APL does not have an edit function. When editing the data *e 

can transmit the APL data out to TSO mode, edit those data by the TSO 
editor and then transmit back to APL mode. The functions IN and OUT , 
listed in Appendix C, can help us to transmit data back and forth. 


It Is not easy to identify the 3-0 vertices of every surface on a 2 : D 
screen. The function "analyze" draws the picture step by step to permit 
one to determine the coordinates of any vertex in the picture. 


Hidden-Line Removal 

A surface is defined by at least three vertices. The identification 
of vertex 1 is arbitrary, but it is important that the numbering of the 
remaining vertices of the surface continue in a counterclockwise 
direction, as viewed from outside the object, facing the surface. We then 
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identify vector u, directed from vertex number 1 to vertex number 2, and 
vector v, directed from vertex number 1 to vertex number 3. The cross 
product, n = u x v, will be normal to the face of the object and will be 
directed outward. 

With each surface of an object we also associate a second vector w, 
a line-of-sight vector. This Is the vector directed from vertex number l 
of the surface to the viewpoint jthe location of the eye of the viewer). 
The dot prjduc£ of yectjrs w and n has the property 
w • n = Iwl | n j cose + + 

where e is the angle between w and n. For a + visible surface, the angle 

between w and n is between 0° and J0°, an<£ w • n is positive. For a 

hidden surface, the angle between w and n is between 90° and 180°, 

and w • n is negative. 

Function "visible" listed in the Appendix is the visibility filter. 
Those surfaces which face the viewer will pass through the filter, 

eventually to be plotted. Those surfaces which do not face the viewer 
will be rejected. 

The method described above permits drawing of the visible portions of 
a single convex object. If we try to draw the visible portions of several 
objects, it will be necessary to identify those portions which are hidden 
by a closer object. We could consider that the hidden portions of an 
object have been filled, but were covered by the filling of a closer 
object. 

The major part of blackout is to determine the drawing sequence of 
every visible surface of the whole satellite consisting of seven 

subparts. The function "hidden" listed in the Appendix plots every 
visible surface in the order of the magnitude of distance between the 
surface and the viewpoint. 


Graphics Models of Solarmax 

The seven component parts for the Solarmax satellite have been drawn 
with top, front and side views using the APL threeviews function as shown 
in Appendix A. These were plotted on a monochrome Tektronic graphics 
terminal (TEK4010). A plot of the entire Solarmax satellite using the 
threeview function on the data set SOLART is shown in Fig. 1. The x, y, z 
axes for the satellite are drawn on this figure as well as the rotations 
required to show the top and front views. These rotation values are 
angles in degrees for rotations about the x, y and z axes. 

The graphical displays for the satellite with hidden-line removal 
have been drawn using a TEK4105 color graphics terminal. An example of 
the results using this procedure are shown in Fig. 2. The allocated data 
sets are "sol art" and "rowse". The execution procedure is 
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Fig. 1: Threeviews of Satellite 







e e 0, hidden solart 
x y z 

where 0 , 0 and 0 are the rotational angles in degrees about the x, y 
and z a&s.^ It cln be seen from these results that it is possible to 
control the relative attitude of the satellite with respect to the viewer. 


Simulation Model 

We now advance to simulate a slowly spinning satellite as viewed by 
an autonomous robot vehicle vision system. Because a perspective 

projection is needed for this simulation, the function visible is 
changed and listed in the Appendix. 

Two different types of simulated vision systems are studied. One of 
them has the camera view fixed as it moves toward the object along a 
certain desired path. We call this type the "fixed vision system" which 
is easily modeled by experiment. Another one, called the "tracking vision 
system", has the camera view adjusted automatically to focus on the object 
no matter what the path is. We study these two different vision systems 
in the following. 


Fixed Vision System 

In this model, the robot approaches the satellite and follows the 
path which is drawn as threeviews. This figure, containing an extra 
oblique projection, is generated by the new function "threeviews" which is 
modified from the original one and listed in the Appendix. Firstly, we 
create the data set TRACE, listed in the Appendix, whicn contains the 
attitude data from the columns 1 to 3 and 3-0 path data from columns 4 to 
6. A function MOTION, listed in the Appendix, reads the data from TRACE 
and generates the images one by one. The results are plotted on one 
picture to show the changes of attitude and size of satellite when the 
camera approaches the spinning satellite (see Fig. 3). 


Tracking Vision System 

In the tracking vision model, the robot camera always focuses on the 
satellite and is adjusted whenever the path is changed. The path data set 
HELIX, listed in the Appendix, contains the spin angle 0 in degrees in 
column 1, and Y and Z positions in columns 2 and 3. J The helix path 
projection on the X-Z plane is governed by the following formula in polar 
coordinates 

Y = ’ 4 0 y + 60 

where 9 is the spin angle (degree) about the Y axis and y is the distance 
between y the camera and the satellite. 

A new function ROBOT, listed in the Appendix, reads the data from 
HELIX and generates the Images sequentially (see Figures 4 and 5). 
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Results of Fixed Vision System 



Fig. 4 



Satellite at 


-180 

-300 

270 
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Preliminary computer graphics images of the Solarmax satellite as 
described above were generated using APL on the Rice University AS9000 
mainframe computer and TEK4010 and TEK4105 graphics terminals. The next 
step of the process of developing the computer graphics modeling involved 
transferring these data to the SUN- III workstation. 


GRAPHICS MODELING ON THE SUN WORKSTATION 

A complete graphics model of the Solarmax satellite has been 
generated on the SUN- 1 1 1 workstation. The model can be manipulated to 
show arbitrary views and different shading effects. The SUN- III 
workstation supports a SunCore package which implements the ACM SIGGRAPH 
core system that conforms to level 3C (dynamic output with 30 scaling, 
rotation and translation) of the core specification for output primitives, 
and to level 2 (complete input) for input primitives. 


Data Generation 


The satellite data were transferred from the IBM 370 to the 
Mechanical Engineering Celerity C1200 computer. These data originally 
were bulky and complicated, since some vertex coordinates which were 
repeated several times were easily Implemented in the APL graphics 
system. To reduce the size of data, a FORTRAN code (Appendix), "trans.f", 
was implemented to accomplish this job. There are two input files needed 
for this code. One is a vertex coordinates data file, "s.s". Another is 
a file "s.n" containing the index of each vertex for each plane. After 
running the program, the output file solar. s is generated in a concise 
form which covers every vertex only once. Based on the data, a program 
"shad" has again been implemented successfully for displaying graphic 
models with choice of different shading. The data of the solar plate is 
entered once as a no-thickness plane, not a solid object. But in an APL 
system, we enter the data twice to simulate two opposite faces. However, 
the "shad" program does work for a no-thickness panel very well. 


Description of Program 

The C program on the SUN workstation Is described in this section. 
The input data sets are allocated as follows: 

no. of vertices no. of surfaces 

W_ x W x W_ y W y w_ z W z * window coordinates 

x~ y z: coordinates of each vertex 

no. of vertices in each surface every vertex index on each surface 

The C program "shad.c" originally written for demonstration purposes 
by Sun Microsystems, has been modified to provide the following menu 
driven screens: 
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Screen 1: Enter the desired shading style 

1) Wireframe display 

2) Gray shading 

3) Gourand 

4) Phong diffuse [default] 

5) Phong specular 
Enter your choice (1-5) 

Screen 2: 

1) Gray 

2) Red 

3) Green [default] 

4) Blue 

5) Yellow 

Enter your choice (1-5) 

Screen 3: Enter the desired display option 

1) Still frame 

2) Rotate the viewer [default] 

3) Rotate the object 

4) Rotate the light source 

5) Quit 

Enter your choice (1-5) 

Note: "rotate" means rotation about vertical y axis. 

Screen 4: Enter the light source position [default: 0.0, 0.0, -1.0] 

x = 

y 3 

z * 

Enter the viewer position [default: 4000, 8000, 6000] 

x 3 

y 3 

z = 

Screen 5: Enter noise level 

1) White noise 

2) Coherent noise 

3) No noise [default] 

Enter your choice (1-3) 

Note* "The light source position" specifies the direction of the light 
source from the object. The direction is expressed in NDC 
(normalized device coordinates). 

This system permits the viewer to control the shading, color, motion, 
light source position and viewer position. 
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Artificial noise can also be introduced. Gaussian white noise or 
coherent noise can be generated and superimposed on the graphics 
picture. This operation permits the simulation of a picture transmitted 
via a noisy channel. 

Wireframe displays of the Solarmax satellite model and a hexagonal or 
hex model are shown in Fig. 6 and the corresponding gray shading models 
are shown in Fig. 7. The hex model is used to illustrate application of 
the graphics model to evaluation of a commercial vision system. 


PHYSICAL MODELING 

A 6" x 12" (base x height) hexagonal box was assembled to simulate 
the body of the Solarmax satellite. This model is supported by a thin rod 
that is attached to the PUMA 560 Robot. The model can be moved by the 
robot arm to simulate the navigation of a satellite in space. 

A communication channel between the PUMA 560 controller and the IBM- 
PC XT was established and a BASIC program that transfers 6 real numbers 
(dx dy, dz, pitch, yaw, and roll) to the PUMA from the PC was implemented 
and* tested successfully. A VAL II program was written which receives the 
six numbers transferred from the PC and moves the PUMA manipulator arm to 
the corresponding position. This system permits control of position and 
orientation of the model. Considerable time and expense is involved in 

building accurate physical models of space objects and precise location of 
the camera and light source require additional effort. On the other hand, 
the position and motion of the computer grahics model can be controlled 
accurately with precise light location. 

Pictures of the hex model in a room with a single bright light source 
are similar to those produced by the graphics system. Results of the 
image processing procedures are comparable for either the physical or 
graphics models. 


COMPUTER VISION SYSTEM 


512 

for 


A 3M VDL (vision development language) vision system is used In the 
work described In this report. This system Is a high-speed vision 
development workstation with a digitizer that provides a resolution of 
x 512 x 8. This system provides a wide variety of vision algorithms 
picture acquisition including filtering, feature enhancement, image 
segmentation and feature extraction, statistical data manipulation, and 
object data base matching. These algorithms are encoded Into macro forms 
(also called primitives) that can be used as building blocks in a high- 
level language program designed for a specific vision application, These 
macros can be used interactively. In a command interpreter called VDL- 
BASIC or by a C compiler using VDL software in a C command subroutine 

library. 
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Figure 6: 
(a) Solarmax 


Graphics Wireframe Models 
Satellite, (b) Hex Cylinder 



a b 

Figure 7: Computer Graphics Gray Shading Models 

(a) Solarmax Satellite, (b) Hex Cylinder 
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vision Primitives 


The upgraded VDL software provides 151 image processing primitives in 
the areas of control, image I/O, transformation, convolution, graphics, 
color, dilation/erosion, segmentation, statistics, image arithmetic, 
threshold, edge detection, transition detection, and neighborhood 
transformations. All these primitives are implemented such that the 
images are processed in the spatial domain. These spatial domain 
implementations have the advantage of fast processing time. Indeed, most 
of the primitives take only a few seconds of processing time. 

Those primitives applicable to the space applications, such as 
location determinations and object recognition, are especially of 
interest. They were studied extensively and are discussed below. 


Edge Detection 


Ten VDL primitives are available 
gray image. Nine primitives, (BIN(), 
VGR(), VDX(), and VDY()), employ 
gradient. The existence of an edge 
given threshold value, 
dimensional and four are 


for finding edges in a binary or 
KIR(), R0B(), S0B( ) , LGR(), HGR( ) , 
the information of the first order 
is asserted if the gradient exceeds a 
Five of the first order edge detectors are two 
one dimensional. One 2-D second order gradient 


edge detector, GCR(), finds edges by means of zero-crossing, 
processing time is about 5 seconds for an image of size 256 x 256. 
be reduced proportionally by reducing the work window size. 


It 


The 

can 


Dilation/Erosion 

Dilation and erosion commands (DIL(), ER0(), etc.) are available for 
both binary and gray images. The maximum number of pixel neighbors that 
are used in the dilation/erosion procedure is limited to 8 pixels. 


Noise Removal 

Two commands are Implemented for the purpose of relieving noise by 
averaging the pixel values. The AVN() command averages the pixel values 
In the neighborhood of 8 pixels. The GAU() primitive performs a weighted 
averaging process. 


Statistics and Transition Detection 


Several statistics commands and transition detectors are very useful 
In determining the position information of an object. The CEN() command 
returns the center of gravity of a binary image. The MXX() command finds 
the right most edge of a gray object. The LXX() command detects the right 
edge along a line, etc. By using the combination of these primitives, 
position and orientation of simple geometries can be easily determined. 


14 


Programming Structure 


The VDL software provides a very limited programming structure. To 
extend the flexibility of the use of the system, a C compiler was added to 
the system. The VDL software provides a C command subroutine library. 
With the command library, all 151 vision primitives can be called from a C 
program. Several C programs, (Appendix), were implemented to examine the 
performance of the primitives under the C language. The only difference 
observed is that the processing time is slightly increased. 

One C program ' f t . c ' was implemented to compute the 2-D Fourier 
Transformation of an image of size 64 x 64. The image is read from the 
frame buffer 1 and the results are stored in two data files. The 
execution time for this program takes about 20 minutes. This implies that 
if frequency domain processing is required by using this system, a long 
processing time must be taken into consideration. 


Image Processing 


Application of the vision system to image processing of the hex model 
is described in this section. Terminology used is defined in a glossary 
at the end of the report. A digitized picture after use of the Sobel edge 
detector is shown in Fig. 8. The noise in this picture can be eliminated 
by additional image processing. 

A histogram for the digitized image is plotted in Fig. 9a. This 
chart displays the number of pixels (picture elements) at each gray level 
from -123 for black to 127 for white. There should be five clusters of 
gray levels for the four sides of the hex and the backgroun, however, 

there is considerable overlapping of the brightness intensities of the two 
sides on the white end of the chart. The dark background can be seen as a 
sharp spike at the left side of the chart representing the darkest 

pixels. This histogram was segmented interactively by mapping ranges of 
gray values corresponding to each of the faces and the background into 
single gray level values using a partial threshold command. Results are 
shown by the segmented histogram of Fig. 9b. 

The following procedure can be used to display the top surface as 
white and all of the other surfaces and the background as black. the 
partial threshold command PTH(start-value, end-value, replace-value) maps 
all pixels with gray levels between start and end values into a single 
gray levels having replace-value. Thus PTH(-128, -90, -128) maps all 
values below -90 into black (-128) and PTH(-60, 127, -128) maps the range 

from -60 to 127 to black. Then PTH(-90, -60, 127) maps the top surface 

into white (127). ADD() combines the white top surface with the image in 
the second frame buffer. When the image of Figure 8 Is in the second 
frame buffer one obtains the result shown in Figure 10. Each of the faces 
can be Isolated in this manner. 
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Figure 8: Hex Model - Output of Sobel Edge Detector 




Figure 9: (a) Histogram of Hex Model Image, 

(b) Histogram after Segmentation 



Figure 10 Display of Top Surface as White after 
Using Partial Threshold Command 
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By adding all of the surfaces using the segmented histogram of fig. 
6b one obtains the image shown in Fig. 11a. Here each of the surfaces 
contains only a single gray level. Now using the Sobel edge detector on 
this final image one obtains the wireframe model shown in Fig. lib. The 
final result of the image processing procedure is to produce a clean, 
noise-free model and a wireframe model that is essentially identical to 
the graphics wireframe model used in the construction of the image. Such 
a wireframe model can be used by vision algorithms to identify objects 
from a library of properties of objects stored in computer memory. This 
exercise illustrates applications of computer graphics modeling to an 
image processing procedure that utilized several of the macros provided by 
the vision system. 


SUMMARY AND CONCLUSIONS 

Computer graphics modeling has been performed using both the APL 
GRAPHPAK and a C program written for the SUN workstation. These models 
permit user control of model location and orientation, light location, 
viewer position, surface shading parameters, color and introduction of 

artificial noise onto the image. A physical model was constructed and 

image processing techniques were applied to both the graphics and physical 
models. It has been found that the graphics models compare favorably with 
the physical model and offer the advantages of being less expensive and 

easier to provide control of the important parameters, such as light 
source and viewer positions. 

Preliminary evaluations of the 3M VDL vision development language 
system Indicated that serious image processing could not be accomplished 
well without the addition of a C compiler. This has been done and 
examples of applications of the vision system to both the graphically 

generated models and the physical model have been presented. 


GLOSSARY 

DIGITIZATION: the process by which analog image signals are converted 

into intensity values and are assigned to pixels of the digitized image. 


PIXEL: short for "PI(X)cture ELement." The smallest resolvable area in 

an image; an individual photosite on a solid state camera. 


RESOLUTION: the measure of a system's ability to represent Images; the 

number of bits of accuracy or number of gray levels that can be 
represented in a pixel. Also, the number of pixel columns and rows in an 

image raster. 


GITIZER (FRAME GRABBER): a frame grabber; the device that samples 

lalog video input data, converts them to digital values and stores them 

i computer memory. 


FRAME BUFFER: computer memory specially designed to hold image data and 

allow for simultaneous video display. 
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Figure 11: (a) Hex Model Image after Segmentation 

(b) Output of Sobel Edge Detector 
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INTENSITY: the degree of brightness of a pixel, determined by the value 

stored at the pixel's address in memory. 

GRAY LEVEL (GRAY SCALE): the intensity of brightness, stored as a 

discrete value in each pixel. 

HISTOGRAM: a chart showing the frequency of occurence of gray levels in 

an image. 

THRESHOLD: a point at which a pass-fail decision is made. More 

specifically, the intensity value below which a pixel is forced to black, 
and equal to or above which a pixel is forced to black, and equal to or 
above which a pixel is forced to white. 

BINARY IMAGE: an image represented in black and white intensities only. 

GRADIENT: the rate of change in pixel intensity linearly over the raster. 

EDGE: regions of an image where intenstiy levels change rapidly, 

representing a border between distinct areas. 

EDGE ENHANCEMENT: process by which the image intensity changes across an 

edge are heightened. 

EDGE DETECTION: process of locating an edge by inspecting a pixel and its 

neighbors and finding rapid changes in intensities. 

SEGMENTATION: the process of separating objects of interest from their 

image context or background. 

NOISE: meaningless data in an image not related to the image and produced 

by a variety of sources, such as the camera. 


SUMMARY OF 3M VDL MACROS 

BIN() = Binary edge detector 

KIR() * Kirsch edge detector 

R0B() ■ Roberts edge detector 

S0B() 3 Sobel edge detector 

LGR() = Largest gradient edge detector 

HGR() * Horizontal gadient edge detector 

VGR() 3 Vertical gradient edge detector 

VDY() 3 Hough vertical gradient 

DIL() = Dilate gray image 

ER0() a Erode gray image 

AVN() = Average neighbor 

GAU() 3 Gaussian (low pass) filter 

CEN() 3 Center of gravity 

MXX() = Rightmost edge of an object 

LXX( ) 3 Rightmost edge pixel along a line 

PTH() 3 Partial thresholding operation 

THR() 3 Thresholding operation 
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APPENDIX A: APL GRAPHPAK MODELING 

Programs, Data and Examples of Output 


Program Description 


HIDDEN 

[1] problem space window y|_|_ X UR ^UR 

lower left upper right 
corner corner 

[2] number of planes 

[3] Index origin = 1 

[4] map problem space into the scaling viewport 

[5] scales data to fit between -1 and 1 

[6] rotate the object through angles e , e and e The angle is 

x y z 

measured in degrees, and the rotation Is counterclockwise. 

[7] create a perspective projection of an object onto the X,Y plane, 
from a vantage point on the Z-axIs, perpendicular to the plane of 
the paper. 

[8] perform a transformation from the scaling viewport to a window in 
problem space. 

[9] comment 

[11] I - 1 

[12] Initialize matrix "judge" 

[13] create the sequence of Integer numbers from rowse [1,1] to rowse 

[ 1 , 2 ]. 

[14] Is this plane visible? 

[15] increment I 

[16] If I < planeno THEN GO TO [13] 


ELSE CONTINUE 


comment 


[17] 

[18] I = 1 

[19] label 

[20] skip the Invisible plane 

[21] plane number 

[22] vertices numbers of the plane 

[23] fill with a color 

[24] draw the picture edges 

[25] comment 

[26] increment index I 

[27] If I < planeno THEN GO TO [19] 

ELSE EXIT 

[28] comment 
VISIBLE 


11] 

> 

u 


[2] 

V 


[3] 

(u X v) 1? > 0? visible 9 test 


[4] 

comment 


[5] 

Z coordinate of the last vertex of 

the plane 

[6] 

export the visible information and 

maximum Z coordinate of the plane 



Threeviews of Satellite 




Threeviews of Module TORSO- 




Computer Programs and Data 


Hidden Line Removal Function 


v h t ddem rrm v 

v e o r hippeh obj p r*s p or’JDOwii p n p i p rowr p rt p r p r ; ,j ; judge 

cn su* “i “i i J. 

[21 p * - « N E M O f R: o w 5 r |; p ;|. :i 

C31 QiOf-l 

[43 BS«-W INTO 5VP 

[53 OBJDOWMf ”1 1 SCALE OBJ 

[ 6 J RfROT ROTATE OBJC'QWK 

[73 RT«-c[2] perspective r 

[83 Tf.CS KFM RT[P 1 2 3J 

C91 n 

Cl 0.1 \ 900 

Clll i«-l 

C 12 1 JUI'GEf ( PLAHEMO y 3 ) f 0 

[ 133 EVALJROWRf ( ( P:CJWSE[ I ? j. 3.0' ) A (RC!W5E[ I p 21 1 ,! ) ) / It 
[143 JUBGE[ I p 3«-( VISIBLE R[ROWRp ] ) , I 

C 1 51 XRI+1 

C 161 -»< x <PLAHEMO ) /EVAL 

[173 n JUDGEf JUDGE[ ( i JUI'GEf j 21 ) f 1 

C 181 XRl 

C 191 plots; 

C201 ^(0=JUDGE[I p 13) /NEXT 
C211 Jf jupge C x f 3 1 

[221 ROWRf( (ROWSE[ Jp 13 <.»•!) a(ROWSE[ Jp21i { 0 >/N 
C231 FILL T[ROWRP3 

[241 CRAW T[ROWRp3 USING COLOR C[13 

C251 fl <□ 

[261 NEXTJIfl+l 

C271 + ( I <PLAUEHO ) /PLOTS 

C231 n 

v 


V VISIBLE [[]3 <7 

V T Mf VISIBLE P p CR1 p CR2 p VIS pSZ 

[13 CRif P[2? 2 31 -PCI 5 2 31 

[21 CR2«-PC3 f 2 31-PCI? 2 31 

[33 visf( (CRici3xCR2C21)-CF:l[21xCR2r:il ) >0 

C41 n 

[53 2SRPC(fPC?ll) f41 
IT 61 '• VHfVI5,ZZ 


Data Transmit Function 


v in rm y 

V flPI-OfTSOF IN COL; AF-L.S; flFLAJfiOU 

Cl 3 AFLSf GCTFILE T1QF 

[23 AF'LAfREPHS AFL5 

[33 F:OW<-fAFLA 

[43 AF LV<- ( F:OW[ 1 3 , COL ) r t , A FLA 

V 


y GETFILE [03 y 
v AfGETFiLE p;c;d;e;cc 
[13 c«-'c 

[23 TSO 'FREE I'D(C) 1 
[33 *C' ALLOCATE E< 

[43 c<-c, • (192 R ' 

[53 CCfC r ' CTL 1 

[63 -»<2 A .= A <-111 Osvo r.<- 2 2 f'c cc')/n 

[73 'SHARE FAILURE' 

C 83 -»0 

C 93 l 1J-»(0v./ a fC,cc)/l4 

C 103 A <- 0 0 f ' ' 

C113 »-2:-K0=f e<-c)/l 4 
C 123 ■» < 0^ c £ ) / L 3 

[133 Affl vc AT E 

C 143 -»«-2 

C153 L35 'READ ERROR, RETURN CODE} ' , T CC 

C 163 L45DSVR d 

V 


v REF-MS [133 V 
y M«-REFM5 MOfXSjXjJ 

[13 MFMO 
C23 iSffM 
C33 IF1 
C43 lpi:uf1 

C53 LF2:-»5Ki X xM[i;u3 ? i'-' 

[63 M[i;J3«-'-' 

[73 sku+skxxmex JJ3* ■ + ■ 
[93 M[i } J 3 f • 0 ' 

[93 sk { 4 LR 2 X I IS C232 J <- J +1 
[103 -H-Plxxisci3>iFi + l 
- y 



V OUT [Q3 V 

v AF'LV OUT TSOF }DIM(BUK fftDY 
[1] flPLVOftflFLV 
[23 DIMffAFLVA 

[33 FLKf.(r'iM[i] r go~r<iM[23)f ' ' 

[43 RDYf.(DiM[i3,eo)f aflva,[23 »lk 

C53 F ri Y PUTFII.F TRHF 


v FUTFILE [[]] V 

? A«-B F UTFILE CpFJ FC} I pMpF: p SpOIO 

C 1 □ n*°<-l 

[23 »«-<"2* 1 1 

C33 F:«-TSO 1 FREE BP ( CARDS424 > ' 

[43 F: <- T 5 O 1 FREE DB<F)' 

[53 Rf-TSO ' ATTRTEi CARDS424 F:ECFM(F B) LRECL(0O) BL K 5 1 ZE < 4240 ) ' 

[63 £«-' M EW 1 

[73 Kf 1 CATALOG 1 

[S3 c ■ b c r> ■ 

[93 §<- ■ space ( 1 0 » 5 0 ) BtKSiZE ( 4240) ' 

[103 ««-' w 1 

C113 Uf- ' (CARD5424) ' 

C123 E*lll 

C133 T '<- ' < > ' PREPARSE C 

[143 -*(4> lff B )/LA 

[153 ■+('() 1 v,*E.[2 4 Y 13 )/I-A 

[163 Lit- < , »[ I 4 P 3 

[173 B «- 5 0 

[103 ->LB 

[193 lajnk^' ' )/L!<" E '[l i 3 
[203 ®«- 1 0 4>' 

[213 LOJ-KO-ltf 1 0 4® ) / L 4 

[223 LB : -» ( 0« f T«-OE<L. AMK S ( OFISALCD MEMBER D ) /OFT5ALCD ) /!_ 1 0 

[233 ®F' ' f T,' ' 

[243 L 10 t -» < Os f TfOBLAHKS ( OPT5ALCK MEMBER D ) /OFT5ALCK ) /I_l 
[253 £«-' ' » T,' ' 

[263 LJ 5 ■+ ( 0( f TfOBLAHKS ( OPT5COHV MEMBER B ) /OPTSCOMV ) /L2 

[273 £<-' 'fTr' ' 

[283 l 2 ; •+ ( 0= ' USING 1 MEMBER r>)/t_3 

[293 U<-<< T ^' ' ) /T«- f D[ ( ( \ 3 > +*DIO ) + B NAMEINBEX • USING 1 P 3 > f ' ' 

[303 L3‘-+<0=' SPACE' MEMBER D)/L4 
[31] T«- 1 2 3 + D MAMEINDEX 'SPACE' 

[323 *♦<'<> ' >/W«-»®[ ( lfT) 3 >/L7 

[333 §<-• ' f ( ' SPACE ' f ( Wpl ' • )/W<-, D[TJ]) f • ' r OBI- AM K S B[l+-lfTP] 

[343 *( ' > < * v.jCCltw) fltw 4 -<w ? f ■ 1 )/WR,B[4+Tp])/U7 

[353 s<-s,' ' r w, • ' 


C36J Mi-»( 0 - ft «- TCO ' allocate DD(F> da < ■ , n , • > ' , d , « , § , • using • ,<_!>/' 

C37D -»( 0 - fi *- TSO 'allocate i:<d(F) na< ' ,>J, ' ) old keef )/L 5 

C3BD 'ALLOCATE FAILURE, C.'ODE=', T A 

C39] -»0 

C40] L5Ji ' r*- ' ' ' r £>9, ' ' ' ' 

[41] ' ' ,rt ' CTL" - 

C42D -»(2 A .=E Dsvo s<- 2 2 f' F F2')/i-6 

C43] 'SHAPE FAILURE, ' 

C44] » 

C45] L 6 ;->( 0 a,=T(.f,fc)/L 7 

C46D 'OPEN FAILURE* ',y t 

C47] » 

C4Q] L7JMMtfB 

C 492— if-D ' 

C50] lqj-kjk i«-h 1 )/L9 
C 51 D f 7 <- C 1 ? .1 
C52] -*(0 = R«-fe>/l8 

C53] 'WRITE ERROR* ' !t P: 

C54] » 

C55] l 9 jr«-dsvr e 


v Tso l'd:i v 

? C< TSO A I B J |”1 1 O ! D 

cm niof-o 

C 2 U -» < 2 ~ c *-arcmd Qevo 'B')/li 

C 3 ] 'TOO SHARE OFFER FAILURE' 

C 43 ->0 

C 5 D Li;r<fn<.fl 

C 6 ] -»( (OtSTire c), a/r ,c> , c<-&)/o 

171 -»(0- c )/0 

C8] ->(99?<c)/L2 

C9D 'TSO RETURN CODE* ' ,( t C), 1 , ' 

C 103 -»L 3 

Cll] L2J ' TSO ABEND code* ' , ' 0 1 23456789«KeDEF>: ' C 3t ( 6f 16 ) T c] 
C12] L3J 'COMMAND} ' f fD 


Pickup Function 


v pickup an v 

<7 MEWfORPER F- 1 C K Ul" OF: <3 } M ? T 

Cl] IK-fORDERCrlD 
[2] MCW<-(M,4)fO 
Z31 *«-l 

ZA1 RPjMnw[;ii]^ORDEF:ci»13»oi !!: GLORE.EF[;i;2D» 2 3 41 

C53 ->(»OI«-I + l)/FF- 

V 


OTORSO 

0 14 

1 17 

1 4 

1 3 

J 14 

0 17 

1 20 

1 5 

1 4 

1 17 
0 20 
1 23 
1 6 
1 5 

1 20 

0 23 

1 26 

1 7 

1 6 
1 23 
0 26 
1 2 ? 

1 8 
1 7 

1 26 

0 29 

1 10 

1 1 
1 8 
1.29 
0 10 
1 11 
1 2 
1 1 
1 10 
0 11 
1 14 
1 3 

1 2 
1 11 




OTon-'F-o r :i c I'ur 

t n r o 


0 

10.575 

” 2. 4449 

2 . 2303 

1 

10.575 

" 2. 4449 

" 2 . 2303 

1 

3.175 

" 2.4449 

"2 . 2333 

1 

3.175 

" 2.4449 

2 . 2333 

1 

10.575 

" 2.4449 

2.2333 

0 

10.575 

" 2.4449 

" 2.2383 

1 

10.575 

"1 . 1266 

" 3.3375 

1 

3.175 

" 1.1266 

" 3,3375 

1 

3.175 

“ 2.4449 

" 2.2383 

1 

10.575 

" 2.4449 

" 2.2383 

0 

10.575 

" 1.1266 

" 3.3375 

1 

10.575 

1 .1266 

" 3.3375 

1 

3.175 

1 .1266 

" 3.3375 

1 

3.175 

" 1 . 1266 

" 3.3375 

1 

10.575 

" 1.1266 

" 3.3375 

0 

10.575 

1 .1266 

" 3.3375 

1 

10.575 

2.4449 

" 2.2383 

1 

3.175 

2.4449 

" 2.2383 

1 

3.175 

1 . 1266 

" 3.3375 

1 

10.575 

1 .1266 

" 3.3375 

0 

10.575 

2.4449 

" 2.2383 

1 

10. 575 

2.4449 

2. 2383 

1 

3.175 

2.4449 

2.2383 

1 

3.175 

2.4449 

” 2.2383 

1 

10.575 

2,4449 

" 2,2383 

0 

10.575 

2,4449 

2.2383 

1 

10.575 

1 . 1266 

3.3375 

1 

3.175 

1 . 1266 

3.3375 

1 

3.175 

2.4449 

2.2383 

1 

10.575 

2.4449 

2.2383 

o 

10.575 

1,1266 

3.3375 

1 

10.575 

"1 , 1266 

3.3375 

1 

3.175 

" 1.1266 

3.3375 

1 

3.175 

1 .1266 

3.3375 

1 

10.575 

1.1266 

3.3375 

0 

10.575 

"1 . 1266 

3.3375 

1 

10.575 

" 2.4449 

2.2383 

1 

3.175 

" 2.4449 

2.2383 

1 

3.175 

" 1.1266 

3.3375 

1 

10,575 

" 1.1266 

3.3375 


Graphics Analysis Function 


^ AM A I TOC LID V 
? miAijjiotj; it f i ; J ; & ; te:mf- 

C 1 D V.+ "1 “1 1 1 

C2D obj<- "0.95 0.95 c A L r ot.J 

C33 RETICLE 
C4] SKETCH OBJ 

C5D *<-□ 

C6D -»<»=•«*■>/<> 

C7] RETICLE 
COD M«-f OBJC ? 1 D 
C9D i«-l 
C10D p:ef:j<-hi 

CUD TEMPf 2 4 fOE'JCifDfOFJC-JJD 
C12D temp ci ? 1D< 0 

C13D SKETCH TEMF 
C14D *<-□ 

C15D -»(B='Q')/o 
C16D I4-J 
C17D -4(I<M)/F:EF 
V 


V THREEVIEW5 £ [] J ‘>7 
v THREEVIEW5 FpW 

E1D n ORTHOGRAPHIC FRO JECT I OHS 

C 23 Wf ”2 ”2 2 2 

C 3 D ”0 ♦ 95 0 ♦ 95 scale: p 

LAI « FRONT VIEW 

E 53 SKETCH ”1 0 translate f 

161 r side view 

SKETCH 1 ”1 o TRANSLATE Q ”90 

C 8 D r top view 

C9J SKETCH "1 I 0 translate 90 0 0 

V 


0 ROTATE P 
ROTATE F‘ 


APL Function "visible" 


V VISIBLE CDD V 

v VtlfV I S X BLE FJ SUM J I f vis; zz 

C13 *<- 2 3 4 1 

£23 SUM«-I CftnPOT P 

C33 x«- 3 4 2 0 

£4 3 SUM«-<I CRXDOT P ) +SUM 

C53 i«- 2 4 3 0 

££3 SUHfSUH-(I CRXDOT P) 

£73 VIS4-SUM ) 0 

C83 2Z<-FC(fFC»lJ) 543 

£93 YNfVZ S f ZZ 

V 


v CRXDOT £03 v 

V SUMf I CRXDOT P}CR1JCR2»VW 

Cl 3 CR14-PC2J (IC13»IC2D)3-^C1» <XC13»IC23)3 

£23 CR2<- p C3J < I CI3r I C23)3- F 'Cl»< I C13» I C23>3 

C33 vw<-o 

C43 -»<0= I C43 )/coh 

£53 vw<-C£23 

C63 cowj 

£73 suMf(( <CRici3xCR2C23)-CRlC23xCR2C13) >x<vw-F£i;i£333) 

v 



The Modified Function "threeviews", including an 
extra oblique projection 


Zll 

Z21 

C3 1 
Z41 
C53 
Z61 
Z71 
C8T 
C93 
C 103 


V THREEVIEWS [0] V 

V THREEVIEWS PfW 

H ORTHOGRAPHIC PROJECTIONS 

H «- ”2 "2 2 2 

R«- “0.93 0.95 scale p 

n FRONT VIEW 

SKETCH - 1 ~1 0 TRANSLATE P 

n SIZ'E VIEW 

SKETCH 1 “1 0 TRANSLATE 0 “90 0 ROTATE P 
H TOP VIEW 

SKETCH “110 TRANSLATE 90 0 0 ROTATE P 
SKETCH 1,3 1,3 0 TRANSLATE 10 10 10 ROTATE P 
RETICLE 

V 


TEST 

0 60 0 0 
10 0 0 
1 0 40 0 

10 0 0 
1 0 0 20 



The Attitude and Path Data 


TRACE 

0 0 0 

0 10 0 

0 20 0 

0 40 0 

0 60 0 

0 75 0 

0 90 0 


130 

"130 

”200 

98 

”98 

”80 

70 

”85 

”20 

40 

"75 

0 

20 

”50 

20 

10 

”25 

30 

0 

0 

40 


TRACEAX 


0 

"130 

130 

200 

1 

”98 

98 

80 

1 

”70 

85 

20 

1 

”40 

75 

0 

1 

“20 

50 

”20 

1 

"10 

25 

"30 

1 

0 

0 

“40 

0 

0 

0 

0 

1 

0 

0 

0 


C 

2 100 


DTOR 

0.017453292 



Function "motion" for Fixed Vision System 


V MOTION[03V 

? MOTION TRACE $ ROW $ X £ ROT f TRAN 

£13 «<- ”30 ”100 100 30 

[23 PLftl-IENOf f ROWSE[ $13 

£33 D IO «-l 

[43 DS*-W INTO SVP 
[53 ROW«-fTRACE[$13 

£63 *«-0 

C73 ref; 

[33 x <• i + 1 

[93 ROTfTRftCE[I $ 1 2 3D 

[103 TRAN*-TRACE[ X $ 4 5 63 
[113 ROT HID TRAM 
£123 4 ( X <ROW)/REP 

V 

V H I !■ [ □ 3 V 

7 ROT HID TRANJOBJDOWMfNf X $ RT $ R $ J $ JUDGE 
C1D A OBJDOWNf ”1 1 SCALE OBJ 

[23 OBJDOWM*-OBJ 

[33 Rf TRAN TRANSLATE ROT ROTATE OBJDOWN 

[43 RT<-C[23 PERSPECTIVE R 

[53 t 4-ds xfm rt[$ 1 2 33 

£63 n TfRT [$ 1 2 33 
£73 ft 

£83 n<-1?00 

£93 x«-i 

[103 JUDGE*- ( PLANEMO r 3 ) f 0 

[113 EVAL‘ROWRf( (ROWSE[I ; 13iN) A(ROW5E[I ; 23i»-*) >/N 
[123 JUDGE[I $ 3*-(VISIBLE R[ROWR$3),I 

£133 If-I + 1 

[143 4 ( I 1PLANENO ) /EVAL 

[153 A "judge*. JUDGE [ (4 JUDGE [ $23 ) $ 3 

£163 X*-1 

[173 flotsj 

£183 4(0><‘J'JDGE[X}13)/NEXT 

C193 J*-JUDGE[IJ33 

[203 ROWR*-< <rowse[J$ 131N)a(R 0WSE[J|23I n > >/ n 
[213 fill t[ROWR|3 

[223 DRAW T[ROWR|3 USING COLOR C[13 

£233 a «-B 

[243 next:i*-x+i 

[253 •* ( I < PL ANENO ) /PLOTS 

£263 a 

V 




r j w r j 


The HELIX Path Data Set 


HELIX 

300 "200 "340 
70 "180 "300 
40 "160 "260 
10 "140 "220 
80 "120 "180 
150 "100 "140 
120 "80 "100 
90 "60 "60 

60 "40 "20 

30 "20 20 
0 0 60 


HELIX AX 

0. OOOOOOOOOEO "2.199999406*2 

1.000000000*0 5 . 615387607 E ”5 

1 . 000000000*0 1 . 800000389*2 

1.000000000*0 2.771281467*2 

1. OOOOOOOOOEO 2.800000000*2 

1 .OOOOOOOOOEO 2.078460875E2 

1 .OOOOOOOOOEO 9.999998919E1 

1 .OOOOOOOOOEO "7 . 4871 834 95*"6 

1 .OOOOOOOOOEO "6 . 000000324 El 

1. OOOOOOOOOEO "6.928203293E1 
1 .OOOOOOOOOEO "4.000000000E1 
0. OOOOOOOOOEO 0. OOOOOOOOOEO 

1 .OOOOOOOOOEO 0. OOOOOOOOOEO 


2.000000000E2 3.810512120E 

1 .800000000E2 4.000000000E 

1.600000000=2 3.117691229E 

1.400000000E2 1.599999697E 

1.200000000E2 "2. 6205 1421 8 E 

1 .OOOOOOOOOE2 "1.200000162E 
8.000000000E1 "1.732050870E 

6 . 000000000* i -i .600000000* 

4.000000000E1 "1 .039230466E 

2.000000000E1 "3.999999892E 

0. OOOOOOOOOEO 0. OOOOOOOOOEO 

0. OOOOOOOOOEO 0. OOOOOOOOOEO 

0.000000000*0 0. OOOOOOOOOEO 


r j m ro ro i r j ro ro ro 


Function ROBOT of Tracking Vision System 


v robot c 03 v 

V ROBOT HELIX *ROW* I fROTf TRAN} R 

[1] «<- -30 "30 30 30 

C23 A boxcls <- 5 3 r 0 50 50 1 "50 50 1 "50 "50 1 50 "50 1 50 50 

£3J FLflHEIIOff ftOWSE[ » 1 D 

C43 Dio«-l 

C53 E'S <- W INTO 5VP 
£63 ROWff HELXX£ } 1 3 

C73 n-o 
C93 REP! 

C93 K-i+l 

C 103 r«-cc23-helix[;i;3D 

Cll 3 ROT«.-57.29578*<"3o(HELXXCiJ23-i-R> > ,helix£X * j 3 , 0 

C 123 traw<-o»0»helixci;33 
C133 pot hibhel tram 
C 143 ‘ • 

C 153 □ 

£163 A FILL(W INTO SVR ) XFM BOXCLS 

C 173 erase 

C183 -»(I < ROW) /REF 

V 

vnidhelcod v 

S7 P;OT HIDHEL TRAN} OS JDOWM | N f Z { RT $ F: f J f JUDGE 
Ell ft OBJDOWNf “1 1 SCALE OBJ 

121 0£< JDOWNf OBJ 

Z21 OBJDOWN*- (Of F:OT£23 f 0) rotate objdown 

C43 OBjr*owNf.(F:OT[;i] f Of 0) rotate objbowm 

C5D ft R*.<TRAN[ J,] , TRAWC23 f 0) TRANSLATE OBJDOWN 

161 R«-TRAN TRANSLATE OBJDOWN 

[7] F:T«-CC23 PERSPECTIVE R 

L8D ft RT«- 0 0 0 translate rt 
E93 t«-ds XFM RTCf l 2 33 
E 103 ft TfRT Cf 1 2 33 
E 1 1 3 ft 

E 123 HM900 

E 133 U1 

£141 JUDGE*- < PL ANENO r 3 ) f 0 

E153 EVALJROWRf ( (ROWSECX J 1 3i« ) A ( POWSE^ I # 23 > )/M 
E163 JUDGED X r ]^(VX5XBLE R[|ROWR $ 3 ) , I 
E l/3 x*-i + i 

E1S3 *4 < X j^PLANENO ) /E VAL 

E193 ft judge«-judgec<*uudgec*23>*3 

E203 

£213 PLOTS? 

C223 -»(0/JU»gecxU3)/next 
C 233 J«-JUDGECiJ33 

£243 ROWR«-( (ROWSE£J} 13iH)A(ROWSEtJ»231 N > >/ N 
C253 PXLL T£ROWRJ] 

C263 r>RAW TCROWR-n USING COLOR C£l] 

C273 a 

£283 NEXT J 1 4-1 + 1 

C293 -»(Xi.RL ANENO) /PLOTS 

C303 a 


Original Data Set Solarmax.vec 
Provided by NASA/ JSC 
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« 

« 

266 

0 3.1750 

1 3.1750 

1 3.1750 

1 3.1750 

1 3.1750 

1 3.1750 

1 3.1750 

1 3.1750 

1 3.1750 

1 10.5750 

1 10.5750 

1 3.1750 

0 10.5750 

1 10.5750 

1 3.1750 

0 10.5750 

1 10.5750 

1 3.1750 

0 10.5750 

1 10.5750 

l 3.1750 

0 10.5750 

1 10.5750 

1 3.1750 

0 10.5750 

1 10.5750 

1 3.1750 

0 10.5750 

1 10.5750 

1 3.1750 

0 10.5750 

1 10.5750 

0 1.9583 

1 1 . 9583 

1 1.9583 

1 1 . 9583 

1 1.9583 

1 1 . 9583 

1 1.9583 

1: 1.9583 

1 1.9583 

1 1 . 9583 

1 1.9583 

1 1.9583 

1 1.9583 

1 -1.9583 

1 -1.9583 

1 1.9583 

0 -1.9583 

1 -1.9583 

1 1.9583 

0 -1.9583 

1 -1.9583 


1 . 1266 

3.3375 

-1 . 1266 

3.3375 

-2.4449 

2.2383 

-2.4449 

-2.2383 

-1.1266 

-3.3375 

1.1266 

-3.3375 ; 

2.4449 

-2.2383 

2.4449 

2.2383 

1.1266 

3.3375 

1.1266 

3.3375 

-1. 1266 

3.3375 

-1.1266 

3.3375 

-1.1266 

3.3375 

-2.4449 

2.2383 

-2,4449 

2.2383 

-2.4449 

2.2383 

-2.4449 

-2.2383 

-2.4449 

-2.2383 

-2.4449 

-2.2383 

-t. 1266 

-3.3375 

-1.1266 

-3.3375 

-1.1266 

-3.3375 

1 . 1266 

-3.3375 

1.1266 

-3.3375 

1.1266 

-3.3375 

2.4449 

-2.2383 

2.4449 

-2.2383 

2.4449 

-2.2383 

2.4449 

2.2383 

2.4449 

2.2383 

2.4449 

2.2383 

1.1266 

3.3375 

-0.2658 

3.4825 

-1.0158 

2.2150 

-1.4025 

1.9925 

-2.8774 

1.9925 

-2.8774 

-1.9925 

-1.4025 

-1.9925 

-1.0158 

-2.2150 

-0.2658 

-3.4825 

3. 1442 

-1.5000 

2.4183 

-0.2417 

2.4183 

0.2417 

3. 1442 

1.5000 

-0.2658 

3.4825 

-0.2658 

3.4825 

-1.0158 

2.2150 

-1.0158 

2.2150 

-1.0158 

2.2150 

-1.4025 

1.9925 

-1.4025 

1 . 9925 

- 1 . 4025 

1.9925 

-2.8774 

1.9925 


19 

07; 34 1986 

» 

so 1 armax. 

vec Page 2 

1 

1.9583 

-2.8774 

1.9925 

0 

-1.9583 

-2.8774 

1.9925 

1 

-1.9583 

-2.8774 

-1.9925 

1 

1.9533 

-2.8774 

-1.9925 

0 

-1.9583 

-2.8774 

-1.9925 

1 

-1.9583 

-1.4025 

-1.9925 

1 

1.9583 

-1.4025 

-1.9925 

0 

-1.9583 

-1,4025 

-1.9925 

1 

-1.9583 

-1.0158 

-2.2150 

1 

1.9583 

-1.0158 

-2.2150 

0 

-1.9583 

-1.0158 

-2.2150 

1 

-1.9583 

-0.2658 

-3.4825 

1 

1.9588 

-0.2658 

-3.4825 

0 

-1.9583 

-0.2658 

-3.4825 

1 

-1.9583 

3. 1442 

-1.5000 

1 

1.9583 

3. 1442 

-1.5000 

0 

-1.9583 

3. 1442 

-1.5000 

1 

-1.9583 

2.4183 

-0.2417 

1 

1.9583 

2.4183 

-0.2417 

0 

-1.9583 

2.4183 

-0.2417 

1 

-1.9583 

2.4183 

0.2417 

1 

1.9583 

2.4183 

0.2417 

0 

-1.9583 

2.4183 

0.2417 

1 

- 1 . 9583 

3. 1442 

1.5000 

1 

1.9583 

3. 1442 

1.5000 

0 

-1.9583 

3. 1442 

1.5000 

1 

-1.9583 

-0.2658 

3.4825 

0 

2.8233 

-6.3191 

4.7183 

1 

2.3267 

-12.8115 

1 . 8275 

1 

2.3267 

-8.1265 

-8.6941 

1 

2.8233 

-1.6342 

-5.8032 

1 

2.8233 

-6.3191 

4,7183 

0 

2.8233 

-4.7574 

1.2108 

1 

2.3267 

-11.2500 

-1.6800 

0 

2.3267 

-9.6882 

-5.1866 

1 

2.8233 

-3. 1958 

-2.2958 

0 

2.8233 

-3.9091 

-0.6950 

1 

2.9166 

-2.7583 

-0. 1667 

1 

2.9166 

-2.7583 

0. 1667 

1 

2.8233 

-4.0441 

-0.3900 

0 

2.6233 

1 . 6342 

5.8032 

1 

2.8233 

6.3191 

-4.7183 

1. 

2.3267 

12.8115 

- 1 . 8275 

1 

2.3267 

8. 1265 

8.6941 

1 

2.8233 

1 . 6342 

5.8032 

0 

2.8233 

3. 1958 

2.2958 

1 

2.3267 

9.6882 

5. 1866 

0 

2.3267 

11.2500 

1.6800 

1 

2.8233 

. 4.7574 

-1.2108 

0 

2.8233 

4.0441 

0.3900 

1 

2.9166 

2.7583 

-0. 1667 

1 

2.9166 

2.7583 

0. 1667 

1 

2.8233 

3.9091 

0.6950 

0 

-3.4583 

0.0000 

2.2083 

1 

-2.4166 

0.0000 

0.0000 

1 

-3.4583 

1.9125 

1.1042 
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1 

- 3.4583 

1.9125 

- 1 . 1042 

1 

- 2.4166 

0.0000 

0.0000 

1 

- 3.4583 

0.0000 

- 2.2083 

1 

- 3.4583 

- 1.9125 

- 1 . 1042 

1 

- 2.4166 

0.0000 

0.0000 

1 

- 3.4583 

- 1.9125 

1. 1042 

1 

- 2.2083 

0.0000 

0.0000 

1 

• 3.4583 

0.0000 

2.2083 

1 

- 3.4583 

1.9125 

1.1042 

1 

- 2.2083 

0.0000 

0.0000 

1 

- 3.4583 

1.9125 

- 1 . 1042 

1 

- 3.4583 

0.0000 

- 2.2083 

1 

- 2.2083 

0.0000 

0.0000 

1 

- 3.4583 

- 1.9125 

- 1. 1042 

1 

- 3.4583 

- 1.9125 

1 . 1042 

1 

- 3.4583 

0.0000 

2.2083 

1 

- 1 . 2333 

0.0000 

0.0000 

1 

- 3.4583 

1.9125 

1.1042 

0 

- 1 . 2333 

0.0000 

0.0000 

t 

- 3.4583 

1.9125 

- 1 . 1042 

0 

- 1 . 2333 

0.0000 

0.0000 

1 

- 3.4583 

0.0000 

- 2.2083 

0 

- 1.2333 

0.0000 

0.0000 

1 

- 3.4583 

- 1.9125 

- 1.1042 

0 

- 1.2333 

0.0000 

0. 0000 

1 

- 3.4583 

- 1.9125 

1. 1042 

0 

- 1 . 2333 

0.0000 

0.0000 

1 

- 2.2083 

0.0000 

0.0000 

0 

- 3.4583 

0.0208 

0.0000 

1 

- 2.3750 

0.0208 

0.0000 

1 

- 2.3750 

0.0000 

0.0208 

1 

- 3.4583 

0.0000 

0.0208 

1 

- 3.4583 

- 0.0203 

0.0000 

1 

- 2.3750 

- 0.0208 

0.0000 

1 

- 2.3750 

0.0000 

- 0.0208 

1 

- 3.4583 

0.0000 

- 0. 0208 

1 

- 3.4583 

0.0208 

0.0000 

i 

- 3.4583 

0.0000 

0.0208 

0 

- 3.4583 

- 0.0208 

0.0000 

1 

- 3.4583 

0.0000 

- 0.0208 

0 

- 2.3750 

0.0000 

- 0.0208 

1 

- 2.3750 

0.0208 

0.0000 

0 

- 2.3750 

0.0000 

0.0208 

r 

- 2.3750 

0.0000 

- 0.0208 

0 

2.8083 

1.0591 

2.5575 

1 

1 . 9583 

1.0591 

2.5575 

1 

1 . 9583 

- 1.0591 

2.5575 

1 

2.8083 

- 1.0591 

2.5575 

1 

2.8083 

1.0591 

2.5575 

1 

3. 1750 

1.1266 

3.3375 

1 

3. 1750 

- 1.1266 

3.3375 

1 

2.8083 

- 1.0591 

2.5575 

0 

1.9583 

- 1.0591 

2.5575 

1 

1.9583 

- 2.5575 

1.0591 

1 

2.8083 

- 2.5575 

1.0591 

1 

2.8083 

- 1.0591 

2.5575 
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0 3.1730 

1 3.1750 

1 2.8083 

0 1.9583 

1 1.9383 

1 2.8083 

1 2.8083 

0 3.1750 

1 3.1750 

1 2.8083 


0 

1 

1 

1 

0 

1 

1 

0 

1 

1 

1 

0 

1 

1 

0 

1 

1 

1 

0 

1 

1 

0 

1 

1 

1 

0 

1 

1 

0 

1 

0 

1 

0 

r 

0 

1 

1 

1 

1 

1 

0 

1 

1 

0 

1 

1 


1.9583 
1.9583 
2.8083 
2.8083 
3. 1750 
3. 1750 
2.8083 
1.9583 
1.9583 
2.8083 
2.8083 
3. 1750 
3. 1750 
2.8083 
1.9583 
1.9583 
2.8083 
2.8083 
3. 1750 
3. 1750 
2.8083 
1.9583 
1.9583 
2.8083 
2.8083 
3. 1750 
3. 1750 
• 2.8083 
1.9583 
1.9583 
2.8083 
2.8083 
3. 1750 
3. 1750 
0. 1400 
0. 1400 
- 0.2792 
0. 1400 
0.0000 
0. 1400 
- 0.2792 
0.0000 
0.0000 
2.3833 
2.5133 
2.5133 


- 1.1266 

- 2.4449 

- 2.5575 

- 2.5575 

- 2.5575 

- 2.5575 

- 2.5575 

- 2.4449 

- 2.4449 

- 2.5575 

- 2.5575 

- 1.0591 

- 1.0591 

- 2.5575 

- 2.4449 

- 1.1266 

- 1.0591 

- 1.0591 

1.0591 

1.0591 
- 1.0591 
- 1.1266 

1.1266 
1.0591 
1.0591 
2.5575 
2.5575 
1.0591 
1.1266 
2.4449 
2.5575 
2.5575 
2.5575 
2.5575 
2.5575 
2.4449 
2.4449 
2.5575 
2.5575 
1.0591 
1.0591 
2.5575 
2.4449 
1.1266 
2.4183 
2.4183 
2.4183 
2.4183 
2.9166 
2.4183 
2 . 41 B 3 
2.9166 
3.7333 
0. 1500 
0.0750 
- 0.0750 


3.3375 

2.2383 

1.0591 

1.0591 

- 1.0591 

- 1.0591 

1.0591 

2.2383 

- 2.2383 

- 1.0591 

- 1.0591 

- 2.5575 

- 2.5575 

- 1.0591 

- 2.2383 

- 3.3375 

- 2.5575 

- 2.5575 

- 2.5575 

- 2.5575 

- 2.5575 

- 3.3375 

- 3.3375 

- 2.5575 

- 2.5575 

- 1.0591 

- 1.0591 

- 2.5575 

- 3.3375 

- 2.2383 

- 1.0591 

- 1.0591 

1.0591 

1.0591 

- 1.0591 

- 2.2383 

2.2383 

1.0591 

1.0591 

2.5575 

2.5575 

1.0591 

2.2383 

3.3375 

0.2417 

- 0.2417 

0.0000 

0.2417 

0.0000 

- 0.2417 

0.0000 

0.0000 

0.0000 

2.7800 

2.7800 

2.7800 
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1 

2.3833 

-0. 1500 

2.7800 

1 

2.2533 

-0.0750 

2.7800 

1 

2.2533 

0.0750 

2.7800 

1 

2.3833 

0. 1500 

2.7800 

1 

2.3833 

0. 1500 

2.5575 

1 

2.5133 

0.0750 

2.5575 

1 

2.5133 

0.0750 

2.7800 

0 

2.5133 

0.0750 

2.5575 

1 

2.5133 

-0.0750 

2.5575 

1 

2.5133 

-0.0750 

2.7800 

0 

2.5133 

-0.0750 

2.5575 

1 

2.3833 

-0. 1500 

2.5575 

1 

2.3833 

-0. 1500 

2.7800 

0 

2.3833 

-0. 1500 

2.5575 

1 

2.2533 

-0.0750 

2.5575 

1 

2.2533 

-0.0750 

2.7800 

0 

2.2533 

-0.0750 

2.5575 

1 

2.2533 

0.0750 

2.5575 

1 

2.2533 

0.0750 

2.7800 

0 

2,2533 

0.0750 

2.5575 

i 

2.3833 

0. 1500 

2.5575 

0 

2.3833 

0. 1500 

-2.7800 

1 

2.5133 

0.0750 

-2.7800 

1 

2.5133 

-0.0750 

-2.7800 

1 

2.3833 

-0. 1500 

-2.7800 

1 

2.2533 

-0.0750 

-2.7800 

1 

2.2533 

0.0750 

•-2.7800 

1 

2.3833 

0. 1500 

-2.7800 

1 

2.3833 

0. 1500 

-2.5575 

l 

2.5133 

0.0750 

-2.5575 

1 

2.5133 

0.0750 

- 2.7800 

0 

2.5133 

0.0750 

-2.5575 

1 

2.5133 

-0.0750 

-2.5575 

1 

2.5133 

-0.0750 

-2.7B00 

0 

2.5133 

-0.0750 

-2.5575 

1 

2.3833 

-0. 1500 

-2.5575 

1 

2.3833 

-0. 1500 

. -2.7800 

0 

2.3833 

-0. 1500 

-2.5575 

1 

2.2533 

-0.0750 

-2.5575 

1 

2.2533 

-0.0750 

-2.7800 

0 

2.2533 

-0.0750 

-2.5575 

1 

2.2533 

0.0750 

-2.5575 

1, 

2.2533 

0.0750 

-2.7800 

0 

2.2533 

0.0750 

-2.5575 

1 

2.3833 

0. 1500 

-2.5575 
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IBM 370 TO CELERITY C1200 



FORTRAN program trans.f 
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program main 

implicit r eal #8 < a-h , o-z ) 

dimension x (500) , y (500) , z (500) , ibf (500) , node (50) 
open (2, f i le 3 ' 5. n J ) 
open (1, filers. s’) 
k = 1 
i = 1 


10 read 

( 1 , t , err=30, end=30) igar, oar. 

x (i ) 

> y <i ) , 

z (i ) 

do 20 

j = 1 , k 




if 

((x(i> .eq. x(j)) .and. (y(i) 

. eq. 

y(j) ) 


. 

. and . < z ( i ) 

. eq. 

Z ( j > ) ) 

then 


if (i .eq. 1) then 
>: 5 = x ( 1 ) 
x 1 = x ( 1 ) 
ys = y ( 1 ) 
yl = y(l) 
z s = z < 1 ) 
zl * z (1) 
ibf (1) = 1 
end i f 

ibf (i ) = j 
go to 25 
endi f 
20 continue 
• k = k + 1 
ibf(i) = k 
x < k ) = x ( i > 
y ( k > = y ( i ) 
z ( k ) = z ( i ) 


if 

<>' (k) 

.It. 

MS) 


= >: ( k ) 

if 

( x ( k ) 

. gt . 

>; 1 ) 

X 1 

= x ( k ) 

if 

( v < k ) 

.It. 

ys) 

ys 

= y ( k ) 

if 

<y <k> 

. gt . 

yl > 

yl 

= y ( k > 

if 

<z (k> 

. It. 

z s) 

z s 

= z v k ) 

if 

(z (k ) 

. gt . 

zl ) 

zl 

= Z ( k ) 

i = i + 1 
go to 10 
write (*, 

50) 

MS, X 1 

j 

ys ? yi? 


50 f or mat (6g 13. 6) 
do 60 j « 1, k 

write <*,50) x(j>, y(j), z(j) 

60 continue 
kv = k 
k - 0 

70 read <2, *,err=99,end=99) n, (node ( j ) , j = l , n) 
k « k + 1 
do 80 j ■ 1 , n 

node(j)=ibf (node ( j ) > 

80 continue 

write (*,100) n, (node ( j ) , j= 1 , n ) 

1 00 format ( 50 i 4 ) 
go to 70 
99 close(l) 
cl ose (2) 

write (#,*) kv, k , ’ <<<<<===== move this line to first row 

end 


cn (.n ui cn tn c.n c .1 -t* *> -t> -t» •*> -t> 4* -c> c-i C J u <4 •-< u <4 u w u to k 1 

cr cn 4» u i j >- o <i co nj O' cn w * j ►* o -o cd si o lh 4 - w ro •-* o 43 co 


Coordinates vertices (before reduction) 


ar 
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1 

0 . OOOOOOOOOeO 


1 . OOOOOOOOOeO 

3 

1 . OOOOOOOOOeO 

4 

1 • OOOOOOOOOeO 

5 

0 . OOOOOOOOOeO 

6 

1 . OOOOOOOOOeO 

7 

1 . OOOOOOOOOeO 

8 

1 . OOOOOOOOOeO 

9 

(j . o 0 0 0 0 0 0 0 0 e 0 

10 

1 . OOOOOOOOOeO 

1 1 

1 . OOOOOOOOOeO 

12 

1 . OOOOOOOOOeO 

13 

0 - OOOOOOOOOeO 

14 

1 • OOOOOOOOOeO 

1 5 

1 - OOOOOOOOOeO 

1 6 

1 . OOOOOOOOOeO 

17 

0 . 0 0 0 0 0 0 0 0 0 e 0 

18 

1 • OOOOOOOOOeO 

19 

1 • OOOOOOOOOeO 

2 0 

1 . OOOOOOOOOeO 

21 

0 . 0 0 0 0 0 0 0 0 0 e 0 

”, , 

1 . OOOOOOOOOeO 

j— i -r 
.j. 

1 . OOOOOOOOOeO 

24 

1 . OOOOOOOOOeO 

^.S J 

0 . OOOOOOOOOeO 

26 

1 - OOOOOOOOOeO 

.•* f 

1 . OOOOOOOOOeO 


1 . O OO 0 OO O O O e 0 
1 . OO OOOOO OOeO 
0 . 000000000 & 0 
1 . OOOOOOOOOeO 
1 . OOOOOOOOOeO 


s . s F’sge 1 


-2. 1 50794 109e- 16 
1 . 912500000e0 
1.91 2500000e0 
-2. 1 50794 109e- 16 
-2. 1 50794 109e- 16 
-6. 03 10478 13e- 16 
1.91 2500000e0 
-2. 1 50794 109e- 16 
-2. 1 50794 109e- 16 
-1.91 2500000e0 
-6. 03 1 04781 3e- 1 6 
-2. 1 50794 109e- 16 
-2. 150794 109e- 16 
-1.91 2500000e0 
-1.91 2500000e0 
-2. 150794109c— 16 
-2. 1 50794 109e- 16 
-6. 03 10478 13e- 16 
—1.91 2500000e0 
-2. 1 50794 109e- 16 
-2. 150794 109e- 16 
1 . 912500000e0 
-6 . 03 1 0478 1 3e- 1 6 
-2. 1 50794 109e- 16 
1.01 5 800000 eO 
1 . 402500000*0 
1 . 402500000e0 
1 . 015800000e0 
1 . 01 5800000e0 
1 . 4025000 OOeO 
1 . OlSSOOOOOeO 
1.01 5800000e0 


1 . OOOOOOOOOeO 
1 . OOOOOOOOOeO 
O . OOOOOOOOOeO 
1 . OOOOOOOOOeO 
1 . OOOOOOOOOeO 
1 . OOOOOOOOOeO 
1 . OOOOOOOOOeO 
0. OOOOOOOOOeO 
1 . OOOOOOOOOeO 
1 . OOOOOOOOOeO 
1 . OOOOOOOOOeO 
1 . OOOOOOOOOeO 
0. OOOOOOOOOeO 
1 . OOOOOOOOOeO 
1 . OOOOOOOOOeO 
1 . OOOOOOOOOeO 
1 . OOOOOOOOOeO 
O. OOOOOOOOOeO 
1 . OOOOOOOOOeO 
1 . OOOOOOOOOeO 
1 . OOOOOOOOOeO 
1 . OOOOOOOOOeO 
0. OOOOOOOOOeO 
1 . OOOOOOOOOeO 


1 . 4 O 2 5 0 0 0 0 0 e 0 
1 . 4 0250 0000 eO 
-2.41 8300000e0 
-2. 41B300000e0 
-2. 418300000e0 
-2 .41 8300000e0 
-2. 418300000e0 
1 . OlSSOOOOOeO 
1.01 SSOOOOOeO 
2 . 658000000e- 1 
2. 658000000*- 1 
1 .015800000e0 
1 . 402500000e0 
2. B77400000e0 
2.877400000e0 
1 . 402500000*0 
1 . 402500000e0 
1 . 402500000e0 
1 . 402500000*0 
2.877400000e0 
2. 877400000e0 
1 . 402500000*0 
1 . 015800000*0 
2. 658000000e- 1 


— 1 . 2 3 3 3 0 0 0 0 0 e 0 
—3. 458300 0 O 0 e 0 
—3. 458300000e0 
- 1 . 233300000e0 
- 1 . 233300000e0 
-3 . 458300000e0 
-3. 458300000e0 
- 1 . 2333 0 0 0 00 e 0 
- 1 . 23330 00 OOeO 
-3 . 458300000e0 
—3. 453300000e0 
-1 . 233300000e0 
- 1 . 233300000e0 
-3.458300000*0 
-3. 458300000e0 
- 1 . 2333 0 0 0 00 e 0 

— 1 . 23330000 0 e 0 
-3 . 458300000e0 
-3. 45S300000e0 
- 1 . 2333 0 0 0 00 e 0 
- 1 . 233300000e0 
-3. 45S300000e0 
-3. 458300000e0 

— 1 . 233300000e0 
1 . 958300000e0 
1 . 958300000e0 

— 1 . 958300000e0 
- 1 . 958300000*0 
1 . 958300000e0 
1 . 958300000e0 
1 . 958300Q00e0 
- 1 . 958300000e0 
— 1 . 958300000*0 
1 . 958300000e0 
1 . 953300000e0 
1 . 958300000e0 
-1 . 958300000e0 
- 1 . 958300000e0 
1 . 958300000e0 
1 . 958300000e0 
-1 . 958300000e0 
- 1 . 958300000*0 
1 . 95B300000e0 
1 . 95B300000eO 
1 . 958300000e0 
1 . 958300000e0 
- 1 . 958300000e0 
- 1 . 958300000*0 
1 . 958300000e0 
1 . 958300000e0 
- 1 . 958300000*0 
- 1 . 958300000e0 
1 . 958300000*0 
1 . 958300000*0 
1 . 95S300000e0 
1 . 958300000e0 


0 

-1 

1 

0 

0 

— 9 


"1 
0 
0 
-1 
_ 'P 
0 

0 

1 

-i 

0 

0 

r> 

1 

0 

0 

1 

0 

0 

n , 

1 
1 


-1 




2 

1 

1 

1 

1 

1 

-1 

-1 

-1 

-1 

-1 


— 


n U OO u o oO (j #=> o 
1 04200000e0 
1 04200000e0 
OOOOOOOOOeO 
OOOOOOOOOeO 
20S300000e0 
1 04200000e0 
OOOOOOOOOeO 
OOOOOOOOOeO 
1 04200000e0 
2 08 300000 eO 
OOOOOOOOOeO 
OOOOOOOOOeO 
1 04 2000 OOeO 
1 04 20000)0 eO 
OOOOOOOOOeO 
OOOOOOOOOeO 
208 30000 OeO 
104 2 000 OOeO 
OOOOOOOOOeO 
OOOOOOOOOeO 
1 04 200 00 OeO 
2083 OOOOO eO 
OOOOOOOOOeO 
2 1 5000 00 OeO 
9 92500 00 OeO 
992500000e0 
2 1 5 00 000 OeO 
21 500000 OeO 
99250 00 00 eO 
2 1 5000000e0 
2 1 SOOOOOOeO 
99250 0 O 0 0 e O 
99 2 5000 OOeO 
4 1 7000000e- 1 
41 7000000®- 1 
4 1 70000 00 e- 1 
41 7000000®- 1 
41 7000000®- 1 
2 1 SOOOOOOeO 
2 1 SOOOOOOeO 
482500000e0 
482500000e0 
21 SOOOOOOeO 
992500000e0 
992500000e0 
992500000e0 
992500000e0 
992500000e0 
992500000e0 
992500000e0 
992500000e0 
992500000e0 
992500000e0 
2 1 SOOOOOOeO 
482500000e0 
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57 1 - 00000000 OeO 

58 1 . OOOOOOOOOeO 

59 1 . OOOOOOOOOeO 

6 0 0 . 0 0 0 0 0 0 0 0 0 e 0 

61 1 . OOOOOOOOOeO 

62 1 - OOOOOOOOOeO 

63 1 . OOOOOOOOOeO 

64 1 . OOOOOOOOOeO 

65 O . OOOOOOOOOeO 

66 1 . OOOOOOOOOeO 

67 1 . OOOOOOOOOeO 
6B 1 . OOOOOOOOOeO 

69 1 . OOOOOOOOOeO 

7 0 0.0 O 0 0 O 0 0 0 0 e 0 

71 1 * OOOOOOOOOeO 

72 1 • OOOOOOOOOeO 

73 1 . OOOOOOOOOeO 

74 1 . OOOOOOOOOeO 

75 1 . OOOOOOOOOeO 

76 1 . O 0 0 0 0 0 0 0 0 e O 

77 1 . OOOOOOOOOeO 

78 1 . OOOOOOOOOeO 

79 1 * OOOOOOOOOeO 

80 1 . OOOOOOOOOeO 

81 1 - OOOOOOOOOeO 

82 1 . OOOOOOOOOeO 

83 1 . OOOOOOOOOeO 

84 1 , OOOOOOOOOeO 

85 1 . OOOOOOOOOeO 

86 1 . OOOOOOOOOeO 

87 1 . OOOOOOOOOeO 

88 1 . OOOOOOOOOeO 
8 9 O - OOOOOOOOOeO 

90 1 . OOOOOOOOOeO 

91 1 . 0 0 O 0 0 0 00 0 e 0 

92 1 . OOOOOOOOOeO 

93 1 . OOOOOOOOOeO 

94 0 . OOOOOOOOOeO 

95 1 . OOOOOOOOOeO 

. 96 1 . OOOOOOOOOeO 

97 1 . OOOOOOOOOeO 

98 1 . OOOOOOOOOeO 

99 0. OOOOOOOOOeO 

1 00 T. OOOOOOOOOeO 

101 1 . OOOOOOOOOeO 

1 02 1 . OOOOOOOOOeO 

1 03 1 . OOOOOOOOOeO 

1 04 0 . OOOOOOOOOeO 

1 05 1 . OOOOOOOOOeO 

1 06 1 . OOOOOOOOOeO 

1 07 1 . OOOOOOOOOeO 

108 1 . OOOOOOOOOeO 

1 09 1 . OOOOOOOOOeO 

1 1 0 1 . OOOOOOOOOeO 

111 1 . OOOOOOOOOeO 

112 1 . OOOOOOOOOeO 


2 . 658000000*- 1 
1 - 01 580 000 OeO 
1 . 01 580 0000 eO 
-2-41 S300000e0 
-2.41 8 30000 OeO 
-3. 1 4 4 200000 eO 
-3 . 1 4 4 200 000 eO 
-2.41 B300000e0 
-2.41 8300000e0 
-3. 144 200 00 OeO 
— 3 . 144200 00 0 e 0 
-2 . 4 1 8300000 eO 
—2.41 8300000e0 
—2.41830 00 0 0 e 0 
-2 . 4 1 8300000e0 
-2 . 4 1 3 30 000 OeO 
-2.41 S300000e0 

-2.91 6600000*0 
— 2 . 4 1 3 3 < 0 0 0 00 e 0 
-2.91 6600000 ©0 
-2.41 3 30000 OeO 
-2.91660 0 0 0 0 © 0 
- 3 . 733 3 0 0 0 00 e 0 
-2.91660 0 0 0 0 e 0 
- 2 . 41830 0 00 0 e 0 
-2.91 6600000e0 
-2. 4 1 83 0000 OeO 
-2.91 6 600000 eO 
-2- 4 183 0000 OeO 
—2.4183 0 0 0 0 0 e 0 
—2 . 4 1 830 0000 eO 
-2.41 8300000e0 
2. 658000000e-l 
2 . 658000000e- 1 
-3. 1 44200000e0 
-3. 144200000e0 
2. 658000000e- 1 
2. 877400000*0 
2. 877400000*0 
"SVB77 4 <5gQ 50eQ 
2. 877 4 OOOOOeO 
2. 877400000*0 
2 . 65S0T5i500de- 1 
-3. r44200000e0 
-3. 144200000e0 
2. 6580 00000e-l 
22658ooddd0e- 1 
2. 658000000*- 1 
-3. 144200000*0 
—2. 418300000*0 
-2.41 8300000*0 
-3. 144200000*0 
2. 658000000*- 1 
1.015800000*0 
1 . 402500000*0 
2. '877400000*0 


-1 . 9 583000 OOeO 
- 1 . 958300000e0 
1 . 958300000e0 
1 . 958300000e0 
- 1 . 958300000*0 
-1 . 9 583 OOOOOeO 
1 . 958300000e0 
1 . 958300000*0 
1 . 958300000*0 
1 . 958300000e0 
-1 . 958 3 OOOOOeO 
-1 . 9583 0 00 O OeO 
1 . 9 58 30 0000 eO 
1 . 400000000e- 1 
1 . 4 0 0 0 0 0 0 0 Oe- 1 
-2. 792000000e-l 
1 . 400000000e-l 
5 . oB6*j>536v 3 1 e 16 
1 . 4 0000 00 00 e-1 
5. 086358631 e-1 6 
-2. 792 0 0 0 0 0 0 e - 1 
5 . 03635863 1 e— 1 6 
6. 51 0629732e— 16 
5. 086358631 e-1 6 
— 2. . V 9 2 0 0 0 0 0 0 © — 1 
5 . 08635863 1 e— 1 6 
1 . 40 0 0 O O 0 0 0 e — 1 
5. 086358631 e-1 6 
1 . 400000000e-l 
-2. 79200 0 0 O 0 e - 1 
1.4 0 0 00 0 0 O 0 e - 1 
1 . 4 00 000000 e- 1 
1 . 958300000e0 
-1 . 958300000*0 
-1 . 958300000e0 
1 . 958300000*0 
1 . 9 5 83 0 0 0 0 0 e O 
1 . 958300000*0 
1 . 958300000*0 
-1 . 958300000*0 
— 1 • 958300000e0 
1 . 95 83000 00*0 
1 . 958300000*0 
1.958300000*0 
-1.958300000*0 
-1, 95 8300000*0 
1 . 958300000*0 
1 . 958300000e0 
1 . 958300000*0 
1.958300660*0 
1 . 958300000*0 
1 .958300000*0 
1 . 958300000*0 
1 ,958300000*0 
1 . 958300000*0 
1 . 958300000*0 


-3. 482500000*0 
— 2 . 21 50 OOOOOeO 
-2 . 2 1 5 000 00 OeO 
-2.41 7000000©- 1 
-2. 417000000©- 1 
-1 . 5 000 OOOOOeO 
-1 . 50000000 OeO 
-2.41 7000000©- 1 
2.41 7000000©- 1 
1 . 500000000*0 
1 . 5 00000 00 OeO 
2.41 7000000©- 1 
2.41 7000000e- 1 
2. 41 7000000*- 1 
-2.41 7000000©- 1 
0 . 0 0 0 0 0 0 0 0 0 e 0 
2.417 OOOOOOe- 1 
0 . 0 0 O 0 O 0 0 0 O e 0 
-2.41 7000000e-i 
0 . OOOOOOOOOeO 
0 . OOOOOOOOOeO 
0 . 0 00 0 0 0 0 0 0 © 0 
0 . 0 0 0 0 0 0 0 0 0 e 0 
0 . 0 (!) <!> 0 (!) 0 00 0 e 0 
(!) . OOOOOOOOOeO 
(!) • 0 0 0 0 00 00 OeO 
-2 . 4 1 7000000©- 1 
0 . OOOOOOOOOeO 
2.41 7000000©— 1 
0 . 0 0 0 0 0 (!) 0 0 0 e 0 
-2.41 7000000©- 1 
2 . 4170 C) (!> (!) 00 © — 1 
3. 4825 OOOOOeO 
3. 482500o00e0 
1 . 500000000*0 
1 . SOOQOOOOOeO 
3. 4 82500000e0 
1 . 992500000*0 
- 1 . 99 25 0 00 0 Oe 0 
-1 . 992500000*0 
1 . 992500000e0 
1 . 9^56660 OeO 
—3. 482B66600e0 
-nB66006006e0 
-1 . 590000000*0 
-3. 482500000e0 
-3. 482560000*0 
3.482500000*0 
1 . 5 00 (!) 0 0 0 0 (!> e (!) 
2.41 7000 00 Oe- 1 
-2 .41 7000000 e- 1 
- 1 . JSOOOOOOOOeO 
-3. 482500000e0 
-2.21 5000000 eO 
-1 . 992500000*0 
-1 . 992560006*0 
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4 

4 

4 

4 

4 

4 

5 
5 
5 
5 
5 
5 
5 
5 
5 
19 


12 3 4 
5 6 7 8 
9 10 11 


w i 

5 

5 

C_’ 

u 

T 

5 

5 

5 

5 

5 

5 

5 

7 


vJ 

5 


13 

17 

21 

25 

30 

35 

40 

45 

50 

55 

60 

65 

70 

89 

94 


14 

18 

22 

26 

31 

36 

41 

46 

51 

56 

61 

66 

71 

90 

95 


15 

19 

23 

27 

32 

37 

42 

47 

52 

57 

62 

67 

72 

91 

96 


12 

16 

20 

24 

28 

33 

38 

43 

48 

53 

58 

63 

68 

73 

92 

97 


99 1 

L 00 

1 0 1 

1 04 

1 05 

106 

117 

113 

119 

122 

123 

124 

127 

128 

129 

132 

133 

134 

137 

138 

139 

142 

143 

144 

147 

148 

149 

152 

153 

154 

159 

160 

161 

164 

165 

166 

169 

170 

171 

174 

175 

176 

179 

180 

1 B 1 

184 

185 

186 

189 

190 

191 

194 

195 

196 

20 1 

202 

203 

206 

207 

208 

211 

212 

213 

216 

217 

218 

22 1 

222 

223 

226 

227 

228 

231 

232 

233 

236 

237 

238 

241 

242 

243 

246 

247 

248 

251 

252 

253 

256 

257 

258 

261 

262 

263 

266 

267 

268 

271 

271 

273 

276 

277 

278 

281 

282 

283 

296 

287 

288 

291 

292 

293 

296 

297 

298 


29 
34 
39 
44 
49 
54 
5? 
64 
69 
74 
93 
98 
102 
1 07 


76 77 


79 SO SI 82 83 84 85 96 87 89 


1 03 
1 08 


109 110 ill 1.12 113 114 115 116 


1 20 

121 

125 

126 

130 

131 

135 

1 36 

140 

I'll 

145 

146 

150 

151 

155 

156 157 

162 

163 

167 

163 

172 

173 

177 

178 

182 

183 

187 

188 

192 

193 

197 

198 199 

204 

205 

209 

210 

214 

215 

219 

220 

224 

225 

229 

230 

234 

235 

239 

240 

244 

245 

249 

250 

254 

255 

259 

260 

264 

265 

269 

270 

274 

275 

279 

280 

284 

285 

289 

290 

294 

295 

299 

300 


158 
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5 

301 

302 

303 

304 305 

5 

306 

307 

o 8 

309 310 

5 

311 

312 

313 

314 315 

5 

316 

317 

318 

319 320 

5 

321 

322 

323 

324 325 

5 

326 

327 

328 

329 330 

5 

331 

332 

333 

334 335 

5 

336 

337 

338 

339 340 

5 

341 

342 

343 

344 345 

5 

346 

347 

348 

349 350 

i? 

35 1 



7C Jl T-C-ET 

sJ 

wJ-L 

O ■-> 

■JJH .JJJ 

5 

j56 

357 

358 

359 360 

5 

361 

362 

363 

364 365 

5 

366 

367 

368 

369 370 

5 

371 

372 

373 

374 375 

5 

376 

377 

378 

379 380 

5 

381 

332 

383 

384 385 

5 

386 

387 

338 

389 390 

9 

end 

391 

392 

393 

394 395 : 


399 


| 


Sattelite data after reduction 
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116 68 


-15.0 15.0 

-15.0 15.0 

-15.0 15 

-0 . 2 1 5079E- 1 5 

-1 . 23330 

0 . 0000001 

1.91250 

-3.45330 

-1 . 10420 

1.91 250 

-3.45830 

1. 10420 

-0. 603105E- 1 5 

-3.45830 

-2. 20830 

-1.91250 

-3. 45830 

-1. 10420 

-1.91250 

-3.45830 

1. 10420 

-0. 603105E— 15 

-3. 45830 

2. 20830 

1 . 01580 

1 . 95830 

2.21500 

1 . 40250 

1 . 95330 

1 . 99250 

1 . 40250 

-1.95830 

1 . 99250 

1.01 580 

-1.95830 

2.21500 

1 . 40250 

1 . 95330 

-1.99250 

1.01 580 

1 . 95830 

-2.21500 

1.01580 

- 1 . 95830 

-2.21 500 


1 . 40250 - 1 . 95830 - 1 - 99250 
-2.41830 1.95830 -0.241700 
-2. 41330 1 . 95330 0. 24 1 700 


-2.41830 

-1.95830 

0.241700 

-2.41 330 

-1 . 95830 

-0. 241700 

O. 2 £>5800 

- 1 . 95830 

3. 48250 

0 . 265800 

1 . 95830 

3.48250 

2. 87740 

1 . 95830 

1.99250 

2.87740 

- 1 . 95830 

1 . 99250 

2. B7740 

-1.95830 

-1.99250 

2.87740 

1 . 95330 

-1 . 99250 

0. 265300 

1 . 95830 

-3. 48250 

0. 265800 

-1.95830 

-3.48250 

-3. 14420 

- 1 . 95830 

— i . 50000 

-3. 14420 

1 . 95830 

- 1 . 50000 

-3. 14420 

1 . 95830 

1 . 50000 

-3. 14420 

- 1 . 95830 

1 . 50000 

-2.41830 

0. 140000 

0. 241700 

-2. 41830 

0. 140000 

-0. 241700 

-2. 41830 

-0. 279200 

0. OOOOOOE+OO 

-2.91660 

0. 508636E-15 

0. OOOOOOE+OO 

—3. 73330 

0. 651063E-15 

0. OOOOOOE+OO 

- 1 . 059 1 0 

2.80830 

2.55750 

1.05910 

2.80830 

2.55750 

1.05910 

1 . 95830 

2.55750 

-1.05910 

1 . 95830 

2 . 55750 

0. 750000E— 01 

2.51330 

2.55750 

0. 750000E-01 

2.51330 

2.55750 

0. 750000E-01 

2.51330 

2. 78000 

0. 750000E-01 

2.51 330 

2.78000 

0. 150000 

2.38330 

2.55750 

0. 150000 

2. 38330 

2.78000 

0.750000E-01 

2.25330 

2 . 55750 

0. 750000E-01 

2.25330 

2.78000 

0. 750000E-01 

2.25330 

2.55750 

■0. 750000E— 01 

2.25330 

2. 78000 

•0 . 1 50000 

2.38330 

2.55750 

0. 150000 

2.38330 

2. 78000 

1.05910 

2.80830 

-2.55750 

-1.05910 

2. 80830 

-2.55750 
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- 1 . 05? 1 0 

1 . 95830 

-2. 55750 

1 . 059 1 0 

1 . 95830 

-2.55750 

— 0. 750000E-01 

2. 51330 

-2.78000 

0, 750000E-Q1 

2.51330 

-2. 78000 

0. 750000E-01 

2.51330 

-2.55750 

-Op 7 5 00 DOE -01 

2.51330 

-2.55750 

0. 150000 

2.38330 

-2. 78000 

0. 150000 

2.38330 

-2.55750 

0. 750000E-0 1 

2.25330 

-2. 78000 

Op 7 50000 E— 01 

2.25330 

—2 . 55750 

“0. 750000E-01 

2 ■ 25330 

-2 . 78000 

—Op 750000E-01 

2.25330 

-2.55750 

— 0. 150000 

2.38330 

-2. 78000 

—Op 150000 

2. 38330 

—2 . 55750 

2. 55750 

2. 80830 

1 . 059 i 0 

2. 55750 

2. 80330 

-1 . 05910 

2.55750 

1 . 95830 

- 1 . 059 1 0 

2.55750 

1 . 95830 

1 . 059 1 0 

-2. 55750 

2 . 80330 

-1 . 05910 

“2.55750 

1 . 95830 

- 1 . 059 1 0 

-2. 55750 

2.80830 

1 . d)59 1 0 

-2. 55750 

1 . 95830 

1.05910 

1 . 63420 

2. 82330 

-5. 80320 

B. 12650 

2. 32670 

-3.69410 

9. 68820 

2. 32670 

—5 . 1 8660 

3. 19580 

2. 82330 

-2.29580 

1 1 . 2500 

2.32670 

- 1 . 68000 

4 . 757 40 

2.82330 

1 . 21080 

12. 31 15 

2. 32670 

1 . 82750 

6.3191 0 

2. 92330 

4.71830 

2.75330 

2.91 660 

-0. 166700 

3. 90910 

2.92330 

-0. 695000 

4. 04410 

2.82330 

—0. 390000 

2.75830 

2.91660 

0. 166700 

-12.8115 

2.32670 

- 1.82750 

-6.31910 

2. 82330 

- 4.71830 

- 4.75740 

2.82330 

-1 . 21080 

- 11.2500 

2.32670 

1 . 68000 

“ 3. 19580 

2.82330 

2.29580 

“ 9.68820 

2.32670 

5 . 18660 

,- 1.63420 

2.82330 

5.80320 

- 8. 12650 

2. 32670 

8.69410 

- 4.04410 

2.82330 

0. 390000 

- 2 . 75 B 30 

2.91660 

- 0. 166700 

- 2.75830 

2.91660 

0. 166700 

- 3.90910 

2.82330 

0.695000 

2.44490 

3. 17500 

2.23830 

2.44490 

3. 17500 
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APPENDIX C: 


C PROGRAM FOR GRAPHICS MODELING 


ON THE SUN WORKSTATION 


shad • c 


/* 


# in elude <us#rcore.h> 
# include < sun/f b i o . h > 
# in elude <math*h> 

# in elude <stdio*h> 

# include "demolib.h" 


*/ 


#def ine MAX VERT 800 

#def ine MAXPOLV 800 

#define MAXPVERT 6400 

static int nvert, npoly? 

static float *wldx, *wldy; 

static float *ndcx, *ndcy; 

static char string [81 3 ! 

static float bbox[33C23? 

static float planeqCMAXPOLV] [43 ; 

static float verti ces [MAXVERT3 [33 ; 

static f 1 oat normal [MAXVERT3 [33 ? 

static int c index [MAXVERT3 ? 

static short normal count [MAXVERT3 ? 

static int npvert [MAXP0LY3 | 

static short epvertptr [MAXP0LV3 ; 

static short pvert [MAXPVERT3 ? 

static float xlist[103, ylist[103, zlist[103? 

static float dxlist[103, dylist[10], dzlist[103; 

static int indx 1 i st [ 103 ; 

static float red [2563 , grn[2563, blu[2563 ? 
static float dred[2563, dgrn[2563, dblu[2563? 
static int rendersty 1 e=l ; 
static int renderhue=0; 

static float xcut [ 2 ]»< 0 . , 1 . > , z eutO [ 2 ] ={0 . , 0- > , zeut [ 23 =< 0 . , 1 . > ? 


ma in < arge , argv) 
int arge; 
char *argv[]? 

■C 

int i, disopt, length *= 0, visibleO; 
static char str[3 = "Enter your choice <l-5> ?“ ? 
float lx,ly,l 2 ,vx,vy,vz,x, y, z, lxr, lyr, lzr? 
float cO, sO, theta, tempi ,temp2? 

if (arge < 2> < pr^ntf < "Usage * shadeobj objfile\n")? exit<l) 
if <getobjdat<argv[13>> exit<2>? 
initvwO ; 

get__view — surf ace <our_surf ace , argv) ; 
start_up_core< ) ; 
setvwpv< ) | 

cycle* 

style_select(>; 
new_f rani* Of 

create_temporary_segment < ) ? 

move_ > abs_2<300« ,700* > ; » 

text”" Enter the desired display option") ; 
mo ve_ab s_2 ( 330 • , 650 • ) { 
textT" 1 > Still frame")? 
move abs_2<. 330. ,62'0. > ? 

t ex t7"2) Rotate the viewer [default!")? 
move_abs„2<’330* > ? 


/* 


text("3) Rotate the object"); 

mo ve_ab s_2< 330 . ,560. > ; 

text("4) Rotate the light source"); 

mo ve_ab s_2(330 • ,530.); 

text<"5) Quit") ; 

mo ve_ abs_2(300 • ,430.); 

text(str); 

inqui re_te>:t_extent_2( str* , wldx , wldy ) ; 
map_world_to_ndc_2(*wldx ,*wldy ,ndcx 9 ndcy) j 
set^echo j'osi ti on (KEYBOARD , 1 , *nd cx , *ndcy ) ; */ 
set_echo_pos it ion (KEYBOARD, 1,0.5,0.47); 
awai t_ keyboard( 1000000000, 1 , string, 41ength) ; 
disopt = at o i ( string) ; 
c 1 ose_t empor ary_segmen t ( ) ; 
if (disopt =** 0) disopt » 2; 
length “ O; 
switch (disopt) ( 
case 1 : 


new_f rame() ; 

ge t xyz < 41 x , 4 1 y , 4 1 2 , &vx , ^vy , 4vz ) ; 
set_ 1 ight_direct ion( lx ,ly , Iz ); 
setvwpo(vx , vy , vz) ; 
new^framtO; 

create_temporary_segment ( ) j 

set^primi t ive^at tributes ( 4 PR I MATTS) ; 
set_po lygon_int er i or_s ty le ( SHADED) ; 
drawob j ( > ; 

close__temporary_5egment( ) ; 
do { 

awai t_keyboard(0 , 1, string, 41ength>; 

> 

while (length » 0); 
b real--: ; 
case 28 

new_frame<); 

getxyz(41x,41y,41z , 4vx , 4vy , 4vz > ; 

set_l ight_direction( lx , ly ,lz ); 
theta 3 0 • ; 
do { 

new_f rame( ) ; 

cO ■ cos(theta); 

sO ■ sin(theta); 

x » vx * cO + vz * sO; 


y ■ vyi 

Z m — sO * vx + vz * cO; 
setvwpo(x,y ,z> ; 
create_temporary_segment() ; 
set_primi tive_ attributes< 4PRIMATTS) ; 
set_po lygon_ inter i or_s ty le( SHADED) ; 
drawob j < ) ; 

c 1 ose_temporary — segment ( ) ; 
theta * theta + .17453 ; 

awai t ^keyboard (3000000 , 1 , string , 41 ength) 

> 

while (length ■* 0); 
break ; 


new^frameO; 

ge t x y z ( 4 1 x , 4 1 y , 4 1 z , 4 v x , 4 vy f 4 v z ) ; 
?.pf 1 i nh t d i recti on < lx , 1 y , 1 7 > ; 


case 35 



set vwp o(vx,vy,v 2 ); 
theta * • 1 7453 ; 
do { 

new_f r ame < ) ; 

c re at extemporary ^segment < ) ; 
set_p r imitive_attributes< A PR I MATTS) ; 
set_po lygon_ i nter i or_s ty le < SHADED) ; 
dr awob j ( ) ; 

close_temporary_segment<) ; 
cO » cos(theta); 
sO * sin(theta); 
x * vx * cO + v: * sO; 
y ■ vy; 

z s — sO * vx + v: * cO ; 
se tvwpo< x , y , z) ; 

1 x r = lx * c 0 + 1 z * s 0 ; 
lyr = ly; 

1 z r * -sO * lx + lz * cO; 

set_l ight_di recti on<lxr , lyr , lzr ); 

for <i « O? i < nvert; i++) { 

tempi * cO * ver t i ces [ i ] CO] + sO * vert i ces C i 3 C2] ; 
temp2 « -sO * ver t i ces C i 3 CO] + cO * vert i ces [ i 3 [2] ; 
vert i ces C i 3 CO] » .tempi; 
vertices [i] [2] ■ temp2; 

> 

awa i t — key bboard (3000000 , 1 , string, Alength) ; 
theta * theta + *17453 ; 

% > 

while < length «* 0); 
b reak ; 

case 4* 

new_f r ame < ) ; 

g e t x y z ( A 1 x , A 1 y , A 1 2 , A v x , 4 vy , A v z > ; 
setvwpo(vx , vy , vz) ; 
theta * 0* ; 
do { 

new_frafni( ) ; 

cO ■ cos(theta); 

sO - sin(theta); 

x - lx * cO + lz * sO; 

y - ly* 

z * -sO *lx+lz*cO| 
set_l ight^direct ion<x ,y ,z >; 
create-.temporary_segment <) ? 
set_primi t ive_at tributes < APR I MATTS) ; 
**t_polygon_interior_s tyle< SHADED) ; 
drawobjO; 

close_temporary_segment<) ; 
theta - theta ♦ .17453 ; 

await_keyboard<3000000, 1, string, Alength); 

> 

while (length ■« 0)? 
break; 

case 5* 

shut_down_core< ) ; 
exit<); 

def au 1 1 * 

goto cycle; 

> 

goto eye 1 e ; 


9 


> 


style_select< ) 

\ 

static char str[] * "Enter your choice (1-5) ?"y 
static char stri[] * "Enter your choice (1-5) 
int done, segnam, pickid, hutnumy 
int hue f length; 

new_frameO; 

i»tvwpv() ; 

cr eat extemporary ^segment < > ; 
move_abs_2<300* v 700* ) ; 

text (" Enter the desired shading style"); 

move_abs_2<330- ,650. > ; 

textC'l) Wireframe display"); 

move_abs_2(330 • ,620. ) ; 

text( H 2> Gray shading"); 

mo ve_ab s_2 ( 330 • , 590 . ) ; 

text < "3) Gouraud" ) ; 

(nove_ab5> < 330 . , 560 • ) ; 

text ("4) Phong diffuse [default]"); 

move -<i abs_2(330. ,530. ) ; 

text < "5) Phong specular"); 

m o v e _ab s _2 ( 300 . ,480. ) ; 

text(str); 

/* » ‘ inquire -it text_ext«nt_2(5tr,wldj( ,wldy> ; 

map.wor 1 d.tojidc^i^wldx , *wl dy , ndcx , nd cy ) ; 
set^echo^pos i t ion (KEYBOARD, 1 ,*ndcx ,*ndcy) ;*/ 

, set_echo_po5itiop (KEYBOARD ,1,0.5,0-47); 

awai t^keyboard( 1000000000, 1, string, ^length); 
close — t emporary_5egment < ) ; 
renders tyle ■= atoi(string) -1; 
if (renderstyle ■* -1) render-style = 3; 
new_frame O ; 

Greatest emporaryx^egment ( ) ; 
movex®fc'*.2<300. ,700. ) ; 

textC'Enter the desired shading color"); 

move_ab s_2<330 • ,650* > ; 

text ( ” 1 ) Gray") ; 

movex«b»_2(330. ,620. ) $ 

text< H 2> Red"); 

mov«x»t»^2(330. ,590. > y 

text<”3) Green [def aul t] M ) ; 

movex*b*_2<330* ,560. > y 

text < ”4) Blue") y 

move ab s_2(330« ,530- > ; 

text7 M 3> Yell ow H > y 

movex»b»x2<300« ,430. ) y 

text (strl ) y 

/* inquire_textx«xtentx2<strl ,wldx,wldy) ; 

map_wor 1 d_t o — n 4 c_2( *wl dx , *wl dy ,ndcx,ndcy)y 
setx® cho_poi i t ion (KEYBOARD, l,*ndcx,*ndcy)y*/ 
se tx® chOxPOi 1 1 i on (KEYBOARD , 1 ,0*52 ,0-47) y 
awa7t_keyboar dClOOOOOOOOO, 1, string, &length)y 
renderhue * ato i ( s t r ing) — ly 
if (renderhue « -1) renderhue • 2; 
c 1 os ext emporary_segmen t ( ) j 

if (renderhue *■ O) def ine_co 1 or^indi ces( ourj urface, 0,255, r ed,grn,blu) 
else def ine^color^indi ces< our.surface ,0 , 255 , dr ed , dgr n , db lu> ; 


/* amb i ent , di f fuse , spe cul ar* , f 1 ood , bump , hue , sty 1 e */ 

hue ■ render hue 9 


switch 

(render sty 1 e ) { 







case 

1 s 

se t_s had ing_par ameters ( 

■ 01, 

•96, *0, 

0. , 

7 • , h u e , 

0)1 

break ; 

case 

2 s 

set_shading_ parameters ( 

.01 , 

• 96 , • 0 , 

0., 

7 . , hue , 

i); 

break ; 

case 

3s 

set_shading_par amet ers ( 

.01 , 

• 95 , . 0 , 

o.. 

7 • f hue, 

2); 

break ; 

case 

4: 

se t_shading_par ame ters ( 

.05, 

.50, .40 

, 0. 

, 7 • , hue 

* 2 > 

; break; 


defaults break; 

> 

> 

getxyz < 1 x , ly , 1 z , vx , vy , vs > 
float *lx,*ly,*lz,*vx, *vy , *vz ; 
< 


s tat i c 

char- 

strl U 

* 

"x» 

" ? 

static 

char 

str2[] 


"y 3 

1 

s tat i c 

char 

str3[] 

* 

M 2=3 


static 

char 

st*-4[] 

m 

H y^ 


s tat i c 

char 

strSC] 

* 

" y » 

" ; 

static 

char 

strt.[] 

s 

n 2 ** 

11 * 


int length; 

setvwpvO; 
set.jfcex t_index ( 1 ) ; 
new.frameO ; 

creat e_t emporary^segmen t ( ) ; 

mov$_abs_2(300* ,700* ) ; text < "enter the light source position [default * 0*0, 0-0, -1 -0] " 
/♦ set_echo_posit i on (KEYBOARD f 1 v -3 , • 4) ;*/ 
move _abs_2( 300 • ,630. >; 
tex t(5trl); 

/* inqui r e_text_ex tent_2< st r 1 , wl dx , wl dy > ; 

map_wor 1 d_ t o_ndc_2( *wl dx , *w 1 dy , ndc >: , nd cy> ; 
s et_e cho_jpos i t i on (KEYBOARD , 1 , *r«dcx , end cy ) ; */ 
s e t_echo^pos i t i on (KEYBOARD ,1 ,0.31 , 0 • 665 > ; 
await_keyboar d( 1000000000, 1 , string, ^length) ; 

if (length » 0) 

* 1 x * 0*0; 
else 

*lx ■ atof (string; ; 
move_ abs_2<300. ,660- > ; 
text<str2>; 

/* b inquire^ text_extent_ 2( str2, wl dx , wl dy ) ; 

* map^wor 1 d^to^jidc^ 2<*wldx , *wl dy,ndcx,ndcy); 
se t_e cho_posi t i on (KEYBOARD, i,*ndcx,*ndcy);*/ 
set_ ec ho_pot it ion (KEYBOARD, 1 ,0*31 ,0*645) ; 

awai t_keyboard( 1000000000 , 1, string, ^length); 

if (length 0> 

*ly ■ 0-0; 
else 

*ly - atof (string) ; 
move_abs_2(300- ,640. ) ; 
text(str3); 

/* inquire_text_extent_2(str3,wldx,wldy); 

map_world_to_ndc_2(*wldx ,*wldy ,ndcx ,ndcy) ; 

*et_ echo_positi on (KEYBOARD , 1 ,*ndcx , *ndcy) ;*/ 
set_echo w pos i t i on (KEYBOARD ,1,0-31 ,0-625) ; 
ema i t ey boar d ( 1000000000 , 1, string, ^length); 

if (length == 0) 

* 1 z - 0-0; 
else 


* 1 z * atof(string); 

if ((*lx « 0.0) ** (*ly « O.O) && (*lz -= 0.0>) { 

*lx * 0.0; 

* ly - 0.0; 

* 1 z - -1.0; 

> 

move_abs_2(300. ,600- ) ; textC'enter the vitwer position [default * 4000, 5000, 
move_abs_ 2 <300. , 580 . ) ; 
text <str4> ; 

/* inqui r e_tex t_ ex t en t_2< s t r4 , wl dx , wl dy > ; 

map^worl d_to — ndc_2(*wl dx , *wl dy , ndcx ,ndcy) ; 
set_echo_posit ion (KEYBOARD, 1 ,*ndcx ,*ndcy ) ; */ 
set_echo_posit ion (KEYBOARD, 1,0*31 ,0*565) $ 
await_keyboar d< 1000000000, 1 , string, ^length) ; 
if (length «« 0> 

*vx «= 0; 
else 

*vx = atof ( string) ; 
move_abs_ 2(300 . ,560*) ; 
t e x t < s t r 5 ) ; 

/* inqu i r e_t ex t_ex ten t_2( s t r 5 , wl dx , wl dy ) ; 

map — wor 1 d_t o_ndc_2( *wl dx , *wldy , ndcx , ndcy) ; 
set_echo_posi t ion(KEYBOARD, 1 , *ndcx , *ndcy ) ;*/ 
set_echo_pos i t i on (KEYBOARD , 1,0*31 ,0*543) ; 
awai t_keyboar d( 1000000000, 1 , str ing, length) ; 

if (length « 0) 

^vy = 0*0; 
else 

*vy = atof (string); 
move_abs_2(300 • ,540* ) ; 
text ( str 6) ; 

/* inqui re^t ex t_ex ten t_2< str6 , wl dx , wl dy ) ; 

map_wor 1 d_t o_ndc_2( *wl dx ,*wl dy, ndcx , ndcy) ; 
set_e cho jpo» i t ion (KEYBOARD, 1 , *ndcx , *ndcy ) ;*/ 
set_ echojpos it ion( KEYBOARD, 1 ,0*31 ,0*527) ; 
awai t_k ey board ( 1000000000, 1, string, ^length); 

if (length » 0) 

*vz - 0.0; 
else 

*vz ■ atof (string) ; 

- if <<*vx " 0.0) 46 <*vy ** 0*0) 44 (*vz mm 0.0)) < 

.♦vx ■ 4000*0; 

avy ■ 3000*0; 

#vz ■ 6000-0; ~ 

7 

close_ tempo rary_ segment ( > ; 

7 

start up^cortO 
t 

int i ; 

float x , y , z ; 

initial iz encore (BASIC, SYNCHRONOUS, THREED); 
our_iur f ace— >cmaps i ze ■ 256; 
our_surf aee->cmapname [03 * 1 \0 f ; 

* if (initial ize_v iew^surf ace <our_ surface, TRUE)) exi t( 1) ; 

select_view_sur face (our _5ur face) ; 

inqui re_co lor _indi ces ( our^surf ace ,0,255 , dr ed , dgrn , dbl u) ; 

for r i — 1 ; i < *255 ; i + + ) { /* load color LUT */ 


6000 ] " > ; 


red[i] « (float)i * 0.003921563; 
grr,[i] = <float>i * 0-003921568; 
b 1 u [ i ] » (float)i * 0.003921563;; 

> 

redCO] * 0.; grn[0] ■ .7; bluCO] * O. ; 

def ine_co 1 or_indi ces < our_surface , 0 f 255 , red , grn f b lu> ; 

/* ambient , diffuse, specular, flood, bump, hue, style */ 
se t_shading_parameter s < .01, .96, .0, 0., 7., 0,0); 

initial ize_devi ce< KEYBOARD, 1 ) ; 

set_echo< KEYBOARD ,1,1); 

se t_e cho_surf ace < KEYBOARD, 1, our^surface); 

setJceyboardO , 30, 1); 

setvwpvO ; 

se t_f ont < 1 ) ; 

> 

shut_down_coreO 

f 

terminate_devi ce (KEYBOARD, 1 > ; 
dese 1 e c t_vi ew^surf ace < our^iurf ac e ) ; 
terminate^vi ew_surf ace <our_iur face) ; 
t erminat encore < ) ; 

> 

int* get ob j dat ( f i 1 ename > 
char ^filename; 

< 

int i, j, k; 

short vtmp, vl, v2 , v3 ; 

float ftmp, maxd, scale, offset[3]; 

float x , y , z , x0,y0 , zO, length ; 

FILE *fptr; 

if <<fptr « f open(f ilename, M r H >) ■* NULL) { 

pr in tf < "Can 1 1 open filet U\n", filename); 
retur n( 1 ) ; 

> 

fscanf<fptr, "fcdfcd" , Anvert, Anpoly); 

if < (nvt rt > MAXVERT) »! (npoly > HAXPOLY) ) < 

printf< M Too many object vertices or polygons\n M > ; 
return<2); 

> * 

fscanf(fptr, M fcf Xt %t%t%t %t M , AbboxC0]C0], AbboxCO]Ci], Abbox[l][0], 
AhboxCl] C 1 3 t Abbox [2] [0] f Abbox C2] C 1 ]> 1 
maxd » 0-0; 

for <i ■ 0; i < 3; i++> { 

of f set Ci] "■ (bboxCi] CO] + bboxCi]d3> / 2.0; 
bboxCi]C0] -■ offsetCi]; 
bboxCiDCl] — of f set [i] « 
if (bboxCUCO] > bboxCi] Cl] > { 
ftmp ■ bboxCi] [6] ; 
bbox^i] CO] » bboxCi] Cl]; 
bboxCi] El] * ftmp; 

> 

if (maxd < bboxCi] [1]) 

maxd = bbox[i] [1]; 

> . ... ' ' .. - ■ 

scale * 1000.0 / maxd; 


for (i 


■ Oj i < 3} i++> { 
bboxCi]CO] *■ scale; 
bboxCi3Cl] ** scale; 

> 

for <i ■ 0; i < nvert; i++) { 

f scant (fptr, "XiXiXt" % Avert i ces [i 3 [0] , Avert i ces [ i ][ 1 ] , 

Avert i ces [i] C2] > ; 

verti ces[i] CO] = (verti ces[i] CO] - offsetCO]) * scale; 

verticesCi] Cl] - (verticesCi] Cl] - offsetCl]) * scale; 

verticesCi] C2] - <vert i cts C i ] C2] - offsetC2]> * scale; 

normal C i] CO] ■ 0*0; normal Ci3 Cl] ■ 0*0; normalCi]C2] ■ 0*0; 

normal count Ci] “ 0; 

> 

k = 0; 

for <i ■ 0; i < npoly; i++> -C 

f scant <fptr, "Xd " , AnpvertCi])! 
if <<k + npvertCi]) > MAXPV/ERT) { 

printfC’Too many polygon verti ces\n") ; 
return (3> ; 

> 

pvertptrCi] * ApvertCk]; t 

for (j * O; j < npvertCi]; j++> < 

f scant <fptr , "%h.d", Avtmp); 

pvertCk++] ■ vtmp - 1; 

> 

planeqCi] CO] - planeqCUCl] * planeqCi]C2] - planeqCi]C3] => 0.0; 
vl = pvertCk _ 13; v2 = pvertCk “ 2]; v3 * pvertCk — 3]; 

* for <j “0; j < 3; j ++ > ( 

planeqCi]CO] +« verti cesCvl] C13 * 

(verti ces Cv2] C2] - verti ces Cv3] C2] ) ; 
planeqCi]Cl] += verti ces Cvl] CO] * 

(vert i ces Cv3] C2] — verti ces Cv2] C2] > ; 
planeqCi] [2] +■ vert i ces Cvl ] CO] * 

<verticesCv2] C13 - ver t i ces Cv3] C 1] > 5 
planeqCi] C3] +■ vert i ces Cvl] CO] * 

((verti cesCv3] Cl] * verti ces Cv2] C2] > - 
(verti cesCv2] Cl] * verti ces Cv3] C2] )> ; 
vtmp ■ vl; vl ■ v2; v2 * v3; v3 ■ vtmp; 

> 

/* if (planeqCi] C3] > 0.0) 

for (j » 0; j <■ 3; j++> planeqCi]CJ3 ■ -planeqC i 3 C i 3 ; 

*/ 

for (j « 1; j <- npvertCi] » J+-0 < /* accum normls */ 

v^mp ■ pvertCk-j]; 

x ■ planeqCi] CO] ; y - planeqCi] Cl] 1 z - planeqCi] C2] ; 

length ■ sqrt( x*x + y*y + z*z); 

normal [vtmp] CO] +• x/length; 

normal Cvtmp] Cl] +■ y/ length; 

normal Cvtmp] C2] z/ length; 

normal count Cvtmp 3 +♦; 

> 

> 

for (i «0; i < nvert; i-*--*-) < 

normal Ci] CO] /■ normal count Ci] ; 
normal Ci] Cl] /■ normal count C i 3 ; 
normal Ci]C2] /■ normal count C i] » 

> 

f closeC fptr) ; 

return (O) ; 


setvwpo <vx, vy, vz ) 
float v::, vy , vz; 

{ 

int i ; 

float diag, del, objdist, near; 

se t_vi ew.jref er ence_po irit ( vx , vy f vz); 
set_vi ew_p lane .normal (-vx , ~vy , — vz); 

5 et_pro j ect i on ( PERSPECT I VE , 0., 0- , OO; 
set^vi ew_.plane._di s tan ce (256 • 0) ; 
if <(vx 0*0) && (vz 0-0)) 

se t^v i ew_up_3(0 • 0 , 0*0, vy); 

else 

5 e t _v i ew_up_3 ( 0 * 0 , 1«0, 0*0); 

s e t_window(-SO • 0 , 80*0, -80*0, 80*0); 
diag * 0*0; 

for < i « 0| i < 3; i++> < 

del « bboxtncn “ bbo«[i3C03; 
diag += del * del; 

> 

diag * sqrt(diag) / 2*0; 

objdist * s q r t < vx*vx + vy*vy + vz*vz); 

near » (diag >* objdist) ? objdist/2-0 s oh j di s t-diag ; 
set - view - depth( near, objdist + diag); 
se t_ window. c 1 ipp in g< TRUE) ; 

* set^f r or.t_plane_cl ipping(TRUE) ; 

set ba ck_p lant.c lipping (TRUE ) ; 

set~vi ewport_3< • 125, .©74, 0- , -749, 0.0, 1-0); 

' > 

static float irvxt orm[4] C4] ? 


drawob j < ) 

i 

int i ; 

float x,y,z,x0,y0,z0, length? 
char ch? 


if (renderetyl* 44 render»tyle<3> 

map ndc to_world_3< -348*, 34y., -870., 4x,4y,4z)l 
m »p_ndc_to_world_3< 0-0, 0.0, 0.0, AxO,4yO,4zO> 1 

x — xo7 y7— yO| z — z0| • 

length « sqrt< x*x + y*y * z*z)» 
if (length !« 0.0) { 

x /- length ? y /“ length? z /- length? 

> 

for <i ■ 0? i < nvert? i++> < 

.nor-lCimJ.z) • tSi-f 

if ( c index [i] < 4> cindexCi] - 4? 
if < c index [i] > 248> cindexCU * 248? 

/* inquire inver«e_compo»ite_matrix(invxform>? */ 

it JrJnO.rltyl. - «>“• 2) > 

else set 2 buffer_cut( our .surface? xcut, zeuto, <->} 
for < i = 0 ; i < npoly; i++> <■ 

/# if <vi sible(plar.eq[x J > > *' 
drawface(i); 


f 


> 


> 


i n t visible (’ p 1 n > 
float p 1 n [ 3 ; 

■C 

i n t i ; 
f 1 oat r ; 

c = 0.0; 

for <i ■ 0; i < 4; i++> 

c +* invxf orm[2] [ i 3 * pln[i]; 
return <c < 0*0); 

> 


drawface(p) int p; 

•C 

int i, j, k; 
short *ptr; 


ptr = pvertptrtpDs 

for C i = 0; i < npvert [p3 ; i++> t 

j ■ *ptr++; 

1 i s t [ i 3 = vertices [j 3 [03 ; ylist[i] “ vert i ces [ j 3 [ 1 3 ; 
zlist[i3 = verti ces [ j 3 C23 ; 
if < render-tty 1 e < 3> { 

i n d x 1 i s t [ i 3 = Cindex[j3j 

>el se < 

d 1 i s t [ i 3 - nor ma 1 [ j 3 [03 } dylist[i3 “ normal [ j 3 C 1 3 5 
dzlist[i3 =* normal [j 3 [23 ; 

> 

> 

if (render-style < 2) { 

if ( r e n d e r h u e > < 

j = ind:-: 1 i st [03 >>'2. ; if <! j > 62.) j — 62} 
set_f i 1 l_inde;:( j + renderhue*64 -63); 

> 

else set_f i 1 1 _index < indxl ist [03 ; 

else if (render-style == 2) {set_ver tex_indi ces< indxlist, npvert [p3 ); > 
else { s e t_vert ex_norma ls<dx list, dy list, dz list , npvert Cp3 ) ;> 
if ( renderstyle 0) polyline _ab s ( xlist,ylist,zlist, npvert [p3 ) t 

else polygon_abs_3< xlist, ylist, zlist, npvert[p3)} 


static float maxvw, vwpp , maxvdim; 
static float vlx, vby, vfz, vdx , vdy, vdz; 

static float minleft, maxright, minbot, maxtop, minback, maxfront; 


ini tvw( > 

i n t i ; 
float ft.npi 


r 


f tmp * 
f o r < i 


Mil 


::vu 


bho:; [0] Cl]; 

* 1 5 i < 3? * 4 +> 

if * b h a s : C i 3 T 1 1 !'■ ftmp> 

flmo - bbow L i ] L 1 3 
- * ftmp; 

* O f jfir * '•* v w ’ *4 : • * 1 ; 


W|'i 



ma ;:vd i m = ma x vw - c e i 1 ( v wp p > ; 

v 1 x = C b b o >: [03 [ 0 ] + m a x v w > / v w p p ; 

v by = 100-0 + (bbo::[l] [ 0 ] + ma;;vw) / vwpp; 

vfz = 580-0 - <bbox [2] [ 1 ] + maxvw ) / vwpp*; 

vdx = (bbox [0] [ 1 ] - bbox [0] [0] > / vwpp; 

vdy = <bbox [13 [1] - bbox [ 1 3 [0] > / vwpp; 

vds = Cbbox [2] [1] - bbox [2] [0] ) / vwpp; 

mi n left = bbox [03 [0] ~ 5-0; 

max right = b b o x [0] [ 1 ] + 5 • 0 ; 

minbot ■ bbox [13 CO] - 5-0; 

max top = bbox[l][l] + 5-0; 

minback * bbox[2][0] - 5-0; 

maxfront = bbox[23Ci] + 5-0; 

> 

set vwpv < > 

< 

set_vi ew_ref ererice^poirit <0- 0, 0-0, 0-0 > ; 

set _v i ew p< lane ^n o r*m a 1 <. O • 0 , O-O, “ 1 • O ) ; 

s e t v i ewj« 1 an e_di stance <‘ 0 * 0 > ; 

setj-roject ion (PARALLEL, 0-0, 0-0, 1 -0) ; 

set _v i e w u p 3 ( 0 • 0, 1 ■ 0 , O ■ 0 .) ; • 

set _w i n d o w <‘ 0 • 0 , 1 023 • (3 , 0 • 0 , 767 • 0 .) ; 
s e t _ v i ewjdep t h < 0*0, 1-0); 

set i n dow_c 1 i p« p i n g < F A L S E > ; 
set _v lew p o r t _3 ( 0 • 0 , 1 ■ , 0 • 0 , .75, 0 - 0 , 1 • ) 

> 

int insideobj <M f y, z > 
float x, y, z; 

< 

if <<x < minleft) !! <x > maxright.)) 
return ( 0) ; 

if <<y < minbot) !! <y > max top)) 
re turn (0) ; 

if <<s < minback) !! <z > maxfront)) 
ret u r n < 0 > ; 
return < 1 > ; 

> 


Data "torso. s 


4 9 1 0 

-3 3 3.0 11-0 “4 4 
2 . 444900000c O 
2 • 444^00000 eO 
2 . 444900000© O 
2 • 444900000e0 
2 * 444900000 eO 
2 . 444900000 eO 

1 - 1 26600000 eO 
1 . 1 26600000 ©0 
2 • 444900000© 0 
2 • 444900000© 0 
1 • 1 26600000 eO 

-1 • 1 26600000 eO 
-1 • 1 26600000 eO 
1 • 126600000e0 
1 . 1 26600000 eO 
• 1 26600000 eO 
—2 • 444900000 eO 
-2 . 444900000 eO 
-I * i 26600000© O 
-1 - 1 26 600000 eO 
—2 • 44490 0000 eO 
-2 . 444900000 eO 
-2 • 44 4 9 00 000 eO 
—2 • 444900000 ©0 
—2 < 444900000 ©0 
-2 -44490 O 0 0 O e 0 
— 1 - 1 26600000 eO 
-1 • 1 26600000 eO 
-2 • 444900000 ©0 
-2 - 444900000 ©0 
”1 . 1 26600000 oO 
1 . 1 26600000 eO 
1 • 12 6600 000 eO 
• 1 26600000 eO 
-1 • 126600000e0 
1 • 126600000e0 
2 • 444900000 eO 
2 • 444900000© 0 
1 • 1 26600000 eO 
1 . 1 26600000 eO 
2 • 444900000e0 

2- 444900000 eO 
% 1 . 126600000©0 

-1 . 1 26600000 ©0 
-2 - 444900000©0 
-2 • 444900000 eO 
-1 - 1 26600000 eO 
1 - 126600000©0 
2 - 444900000e0 
3 1 2 3 4 3 

5 6 7 9 9 10 

3- 11 12 13 14 

3 16 17 19 19 

3 21 22 23 24 

3 96 27 29 29 

3 ' 1 3;* 7T M 

3 36 37 38 39 

9 41 4* 43 44 


1 - 057500000© 1 
1 • 057500000 e 1 
3 . 1750 O O O O O e O 
3*1 75000000 eO 
1 • 057500000 e 1 
1 -057500000 el 
1 . 057500000© 1 
3. 175000000©0 
3. 1 75000000 eO 
1 • 057500000© 1 
1 - 057500000e 1 
1 .057500000 el 
3. 1 75000000 eO 
3 . 1 75000000 eO 
1 .057500000© 1 
1 - 057500000© 1 
1 - 057500000 e 1 
3. 1 75000000 eO 
3 • 1 75000000e 0 
1 .057500000 el 
1 - 057500000 e 1 
1 - 057500000 e 1 
3. 1 75000000 eO 
3 « 1 75000000 eO 
1 .057500000 el 
1 • 057300000 e 1 
1 • 057500000 e 1 
3* 175000000©0 
3. 1 75000000 ©0 
1 • 057500000 e 1 
1 • 057500000© 1 
1 .057500000 el 
3 - 175000000e0 
3 • 1 75000000 eO 
1 .057500000©! 

1 . 05750i>000© 1 
1 • 057 500000 e 1 
3* 175000000e0 
3 . 175000000e0 
1 . 057500000e 1 
1 .057500000 el 
1 . 057500000© 1 
1 .057500000 el 
1 .057500000© 1 
1 • 057500000 e 1 
1 • 057500000e 1 
1 • 057 500000 e 1 
1 • 057500000e 1 
1 • 057500000 e 1 

15 

20 


40 

45 46 47 48 49 


2 - 23 S 3 0 0 0 O O e 0 
—2 • 2393000 00 ©O 
-2 • 238300000 eO 
2 • 238300000 eO 
2 • 238 3 0 0 00 0 e 0 
—2 • 238300000 eO 
-3 . 337 500000 eO 
-3 • 3 3 7500000 eO 
—2 • 239300000 eO 
—2 • 238300000 eO 
-3 • 337500000e0 
—3 . 337500000e0 
-3 • 337500000 eO 
-3 . 337500000e0 
-3 • 337500000© 0 
-3 . 337500000 ©O 
-2 - 238300000 e 0 
-2 * 238300000 ©O 
- 3 . 337500000 eO 
-3 . 337500000 eO 
- 1 . 238300000 eO 
2 * 238300000 ©0 
2 . 238300000e0 
-2 • 238300000 ©0 
—2 . 23 8300000 eO 
2 .238300000 eO 
3 -337500000e0 
3 . 33 7500000 ©O 
2 - 238300000 eO 

2 - 2 3 8 3 O O 0 0 0 © O 
3-3375 0 0 O O O e 0 
3 • 337500000 ©O 
3 • 337500000e0 
3 • 337500 OOOeO 

3 - 337500000 eO 
3 .337500000©© 
2 • 238300000e0 
2 . 23830 0000 e O 
3 . 337500000e0 
3.33750000090 

-2 . 23S300000e0 
2 . 239300000 ©O 
3 • 337500000e0 
3. 337500000 eO 
2. 238300000 eO 
-2 • 238300000e0 
-3.337500000e0 
-3 - 337500000 ©O 
-2 • 238300000 eO 


Noise Program: temp.c 




************************** 


* Thu Aug 13 12:08:59 1987 

* / julius/ j 2 / .yhl/nasa/temp.c 

************************************************************** 
/ m, i, t ii t nil 1 1i ************************************************** * 

/* temp.c 
/* alias cscr tem 

/* f 77 -£68881 -o !* !*.c -lfortd -lcore -lsunwindow -lpixrect 
/* width - 462(7 sin cycles), height - 325 
/* Output is to 19-inches bit-mapped color display having 
/* 1152 by 900 pixels 

/************************************************************* 


******* 

******/ 

*/ 

*/ 

-lm */ 

*1 

*/ 

*/ 

******/ 


♦include <usercore.h> 

# include "demolib.h" 

♦include <sun/fbio.h> 

♦include <math.h> 

♦include <stdio.h> 

♦include <sys/file.h> 

♦include <pixrect/pixrect_hs . h> 

♦define COLOR_VWSURF (ddname) { 0, ddname, 0, 256, 0, 0 I n- 

struct pixrect *screen; 


♦define MAXVERT 500 
♦define MAXPOLY 500 
♦define MAXPVERT 3000 


static int 
static float 
static float 
static char 
static float 
static float 
static float 
static float 
static float 
static int 
static short 
static int 
static short 
static short 
static float 
static float 
static int 
static float 
static float 
static float 


nvert, npoly; 

*wldx, *wldy; 

*ndcx, *ndcy; 
string[81] ; 

lx, nix, ly, nly, lz, nlz; 

bbox[3] [2] ; 

planeq [MAXPOLY] [4] ; 

vertices [MAXVERT] [3] ; 

normal [MAXVERT] [3]; 

cindex [MAXVERT] , cpindex [MAXPOLY] ; 

normalcount [MAXVERT] ; 

npvert [MAXPOLY] ; 

*pvertptr [MAXPOLY] ; 
pvert [MAXPVERT] ; 

xlist [MAXVERT] , ylist [MAXVERT] , zlist [MAXVERT] ; 

dxlist [MAXVERT] , dylist [MAXVERT] , dzlist [MAXVERT] ; 

indxlist [MAXVERT] ; 

bred [256] , bgrn[256], bblu[256]; 

cred[256], cgrn[256], cblu[256]; 

ambient, diffuse, specular, flood, bump/ 


static int 
static int 
static float 

), 


renderstyle - 1 ; 
renderhue » 0 ; 
xcut[ 2 ] » [ 

0., le 

zcut 0 [ 2 ] ■ { 

0 ., 1 . 



zcut[2] » { 

0 ., 1 - 

} ; 

int bw2dd(); 

^ nt pixel [560] [440] , width, height; 

struct vwsurf vwsurf ■ 

COLOR_VWSURF (bw2dd) ; 

struct suncore_raster raster; 

autodraw, debug, noise; 

void normligO; 


MAIN_() 

{ 

void 

double 

double 

float 

float 

int 

int 

int 

short 

int 

float 

static char 
float 
float 
char 


blkscr ( ) ; 
dseed, dvalue; 
unit, amp - 14 . ; 
ggnqf__(), fvalue; 
factorl, factor- 
cycles - 70; 

1# j; 

axO, ayO, ax, ay; 

temp, vs in, vdsin, vnormal; 

disopt, length * 0, visibleO; 

angle - 0.785398163; /* 45 degree ^*/ 

st r [] ” "Enter your choice (1-5) ?"; 

vx, vy, vz, x, y, z, lxr, nlxO, lyr, nlyO, lzr, 

C 0, ccO, sO, csO, theta, ctheta, tempi, temp2; 

argv(], yn; 


printf ("objfile : ") ; 
scanf ("%s", argv) ? 
printf ("noise ? (1/0); ") ; 
scanf <"%d", snoise) ; 

printf ("autodraw ? (1/0): ") ; 

scanf ("%d", iautodraw) ; 
printf ("debug ? (1/0): "); 
scanf ("%d" , fidebug) ; 


getob jdat (argv) ; 
start_up_core () ; 


inquire_color_ 

indices 

for (i ■ 1; 

i 

'<- 255; 

bred[i] 

- 

(float) 

bgrn[i] 

■c 

(float) 

bblu[i] 

- 

(float) 

bred[i] 

- 

(float) 

bgrnti] 

- 

(float) 

bblu [i] 

m 

(float) 


fi vwsurf , 0 , 255, cred, cgrn, cblu) , 
i++) ( /* load color LUT */ 

i *0 . 00392156B; 
i *0.003921568; 

*0.003921568; */ 


*0.0033064207 + 0.15686272; 
*0.0033064207 + 0.15686272; 
*0.0033064207 + 0.15686272; 


cycle : 

cred[0] ■ 0.6; 


nlzO 


cgrn(0] - 0.0; 
cblu[0] - 0.6; 

defi.ne_color_indi.ces (fivwsurf , 0, 255, cred, cgrn, cblu) ; 
aet_text_index (255) ; 

/* 

* bred[0] - 60; bgrn[0] - 60; bblu[0] - 60; bred(255] - 

* bgrn [255] - 40; bblu[255J - 40; define_color_indices 

* 255, bred, bgrn, bblu) ; 

* ambient, diffuse, specular, flood, bump, hue, style 

* aet_shading_pa rameters (.01,. 96, .0, 0., 7., 0, 0) ; 

*/ 

style_select () ; 
new_frame () ; 

create_temporary_segment () ; 
move_abs_2 (300 . , 700.); 

text ("Enter the desired display option") ; 

move_abs_2 (330., 650.); 

text("l) Still frame"); 

move_abs_2 (330 . , 620 . ) ; 

text ("2) Rotate the viewer"); 

mo ve_abs_2 (330., 590.); 

text ("3) Rotate the object [default]"); 

move_abs_2 (330 . , 560.); 

text ("4) Rotate the light source"); 

move_abs_2 (330 . , 530.); 

text ("5) Quit"); 

move_abs_2 ( 3 0 0 . , 480.); 

text (str) ; 

set_echo_pos it ion (KEYBOARD, 1, 0.5, 0.47); 
await_keyboard(1000000000, 1, string, & length) ; 
disopt - atoi (string) ; 
close_temporary_segment () ; 
if (disopt -- 0) 

disopt “ 3; /* default 00000 */ 

length “0; 
if (debug) ( 

p r ^ nt £ ("\n\n* ******** *********** ************", 

«**************************************\n 

) 

switch (disopt) { 
case 1: /* Still Frame */ 

new_f rame () ; 

getxyz(tlx, tly, *1*, Svx, fivy, 6vz) ; 
set_light_direction(nlx, nly, -nlz) ; 
aetvwpo(vx, vy, vz) ; 
if (debug) ( 

printf ("Still frame\n") ; 

printf ("view ; %f,%f,%f\n", vx, vy, vz); 
printf ("light : %f,%f,%f\n", nix, nly, nlz) ; 

) 

create_temporary_segment 0 ; 


40; 

(fivwsurf, 0, 


/* sets all the primitive attributes to their default values */ 
aet_polygon_interior_style (SHADED) ; 
if (Inoise) { 

new_f rame () ; 
drawob j ( > ; 

) else ( 

screen ■ pr_open("/dev/bwtwol") ; 

blkscr () ; 

drawob j 0 ; 

width ■ 560; 

height - 440; 

axO - (1152 - width) / 2; 

ayO - (900 - height) / 2; 

for (j - 0; j < height; j++) { 

for (i “ 0; i < width; i ■ i++) { 

ax - axO + i; 
ay - ayO + j; 

pixel(i][jj ■ pr_get (screen, ax, ay); 
pr_put (screen, ax, ay, 0); 

) 

) 

blkscr () ; 

unit - 6.283185307 / width * cycles; 
dseed ■ 123457; 

axO - 8; 
ayO - 5; 

for (j - 0; j < height; j++) { 

for (i ■ 0; i < width; i * i++) { 

temp - pixel [i] [j]; 
vsin - temp + (short) (amp * 

sin (unit * (double) (i) ) ) ; 
vdsin ■ temp + (short) (amp * 

sin ((double) i * unit) * 
sin ((double) j * unit)); 
dvalue “ ggnqf_(&dseed) ; 
fvalue ■ *( (float *) fidvalue) ; 
vnormal ” temp + (short) (amp * fvalue) ; 
ax - axO + i; 
ay - ayO + j; 

pr_put (screen, ax + 576, ay, temp); 
pr_put (screen, ax, ay, vsin) ; 
pr_put (screen, ax, ay + 450, vdsin); 
pr_put (screen, ax + 576, ay + 450, vnormal); 

} 

> 

scanf("%d", si); 
pr_close (screen) ; 
free raster (Sraster) ; 


) 


close_temporary_3egment () ; 
do ( 

await_keyboard (0, 1, string, 4 length ) ; 

} 

while (length -- 0); 
break; 

case 2: /* rotate the viewer */ 

new_f rame 0 ; 

getxyz(&lxr, 4lyr, &lzr, Svx, ivy, Svz) ; 
aet_light_direction (nix, nly, -nlz) ; 
if (debug) { 

printf ("light : %f , %f , %f \n", nix, nly, nlz); 

) 

setvwpo(vx, vy, vz) ; 

theta - 0 . ; 

lx “ nix; 

ly - nly; 

lz - nlz; 

x ■ vx; 

y - vy; 

z ■ vz; 

if (debug) { 

printf ("rotate the viewer\n") ; 

) 

do { 

new_f rame ( ) ; 
if (debug) ( 

printf ("view : %f , %f , %f \n", x, y, z) ; 
printf ("light : %f,%f,%f\n", lx, ly, lz) ; 

) 

create_temporary_3egment () ; 
set_polygon_interior_style (SHADED) ; 
drawob j () ; 

close_temporary_segment () ; 
theta * theta + angle; 
cO - cos (theta); 
sO - sin (theta); 
x ■ vx * cO + vz * 30; 
y - vy; 

z - -sO * vx + vz * cO; 
setvwpo(x, y, z) ; 
lx - nix * cO - nlz * sO; 
ly - nly; 

lz - -sO * nix - nlz * cO; 
set_light_direction (lx, ly, lz); 

/* 

* for (i - 0; i < nvert; i++) ( tempi - cO * 

* vertices [i] 10] + aO * vertices (i] [2] ; temp2 - - 

* vertices [ i] (0] + cO * vertices [i] [2] ; 

* vertices I i] (0] - tempi; vertices I i] (2] - temp2; 


30 * 
) 



*/ 

await_keyboard ( 3000000, 1, string, ilength) ; 

} 

while (length -- 0); 
break; 

case 3: /* rotate the object */ 

getxyz (6lxr, Slyr, ilzr, 6vx, ivy, ivz) ; 
set_light_direction (nix, nly, -nlz); 
if (debug) { 

print f (" light : %f,%f,%f\n", nix, nly, nlz); 

) 

theta - 0 . ; 
ctheta - 0 . ; 

nlxO “ nix; 
nlyO - nly; 
nlzO - nlz; 
if (debug) { 

printf ("rotate the object\n") ; 

) 

do { 

new_f rame () ; 
cO - cos (theta); 
sO - sin (theta); 
ccO - cos (ctheta); 
csO “ sin (ctheta); 

x “ vx * cO + vz * sO ; 
y - vy; 

z " -sO * vx + vz * cO; 
setvwpo(x, y, z) ; 
nix - nlxO * ccO + nlzO * csO; 
nly - nlyO; 

nlz - -csO * nlxO + nlzO * ccO; 

/* set_light_direct ion (nix, nly, -nlz); */ 
if (debug) ( 

printf ("light : %f,%f,%f\n", nix, nly, nlz) 
printf ("view : %f,%f,%f\n", x, y, z) ; 

1 

create_temporary_segment () ; 
set_polygon_interior_style (SHADED) ; 
drawobj () ; 

close_texnporary_segment 0 ; 
theta « theta + angle; 
ctheta ■ ctheta + angle; 

await keyboard ( 30 0 00 00 , 1, string, (length) ; 

) 

while (length — 0); 
break; 

case 4: /* rotate the light source */ 

new frame (); 



getxyz (i lxr, ilyr, ilzr, ivx, ivy, ivz); 
setvwpo(vx, vy, vz); 
if (debug) { 

printf ("rotate the light source\n") ; 
printf ("view : %f,%f,%f\n", vx, vy, vz); 

) 

theta - 0 . ; 
nlxO “ nix; 
nlyO - nly; 
nlzO - nlz; 

do { 

new_f rame ( ) ; 

cO - cos (theta); 

sO * sin (theta); 

nix “ nlxO * cO + nlzO * sO; 

nly - nlyO ; 

nlz - -sO * nlxO + nlzO * cO; 
if (debug) { 

printf ( "vertices [ %d] : (%f,%f,%f) : 

indxlist - %d\n", i, xlist[i], ylistfi], zlist(i], 
indxlist [i] ) ; 

) 

create_temporary_segroent () ; 
set_light_direction (nix, nly, -nlz) ; 
set_polygon_interior_style (SHADED) ; 
drawobj () ; 

close_temporary_segment () ; 
theta - theta + angle; 

a wait_keyboard (3000000, 1, string, slength); 

) 

while (length — 0) ; 
break; 
case 5 : 

shut_down_core ( ) ; 
exit () ; 
default : 

goto cycle; 

) 

goto cycle; 

) 


3 tr{] • "Enter your choice (1-6) ?"; 
strl[] ■ "Enter your choice (1-5) ?"; 
done, segnam, pickid, butnum; 
hue, length; 

new_frame 0 ; 
setvwpv O ; 

create_temporary_segment () ; 

set text index(25i5); 


style_select () 

{ 

static char 
static char 
int 
int 


move _ a k s _2 ( 30 0 . , 7 00.); 

text ("Enter the desired shading style"); 

move_abs_2 (330., 650,); 

text("l) Wireframe display"); 

move_abs — 2 (330 . , 620.); 

text ("2) Gray shading [default]"); 

move_abs_2 (330 . , 5 90 . ) ; 

text ("3) Gouraud") ; 

move — abs_2 ( 330 . , 560.); 

text ("4) Phong diffuse"); 

jnove_abs_2 (330., 530.); 

text ("5) Phong specular"); 

move_abs_2 (330 . , 500.); 

text ("6) Quit"); 

move_abs_2 (300 . , 480 . ) ; 

text(str); 

set_echo_posit ion (KEYBOARD, 1, 0.5, 0.47); 
awaitjceyboard(1000000000, 1, string, & length) ; 
close_ temporary_segment {) ; 
renderstyle » atoi (string) - 1; 

/* default */ 
if (renderstyle ” -1) 
renderstyle • 1; 
if (renderstyle >- 5) ( 

shut_down_core <) ; 
exit 0 ; 

) 

new_f rame ( ) ; 

create_temporary_segment () ; 
mo ve _abs_2 (300., 700.); 

text ("Enter the desired shading color"); 

move_abs_2 (330 . , 650 . ) ; 

text("l) Gray [default]"); 

m° v e_abs_2 (330., 620.); 

text ("2) Red"); 

move_abs_2 (330 . , 590 . ) ; 

text ("3) Green"); 

m°ve w abs_2 (330., 560.); 

text ("4) Blue"); 

move_ abs_ 2 (330 . , 530 . ) ; 

text ("5) Yellow") ; 

move_abs_2 (300 . , 480 . ) ; 

text (strl) ; 

set_echo_ position (KEYBOARD, 1, 0.52, 0.47); 
a wa i t_keyboa rd (1000000000, 1, string, ^length); 
move_abs_2 { 30 . , 500 .) ; 

text ("ambient (0.00) , diffuse ( 1 . 00) , specular (0) , flood (0) , bump (7) ?") 
scanf ("%f %f %f %f %f ", & ambient, ^diffuse, aspecular, Gflood, &bump); 
close_temporary_segment () ; 
renderhue ■ atoi (string) - 1; 
if (renderhue — -1) 



renderhue - 0; /* default 00000 */ 

/* ambient , diffuse, specular, flood, bump, hue, style */ 
hue - renderhue; 

if (ambient » 0 LL diffuse -- 0 LL specular — 0 
LL flood -« 0 LL bump *• 0) ( 
ambient - 0.00; 
diffuse - 1.00; 
specular - 0.0; 
flood - 0.0; 
bump - 7.0; 

} 

switch (renderstyle) ( 
case 1: 

set_shading_parameters (ambient, diffuse, specular, flood, 

bump, hue, 0) ; 

break; 
case 2: 

set_shading_parameters (ambient, diffuse, specular, flood, 

bump, hue, 1) ; 

break; 
case 3: 

set^shading^parameters (ambient, diffuse, specular, flood, 

bump, hue, 2) ; 

break; 
case 4: 

set_shading_jparameters (ambient, diffuse, specular, flood, 

bump, hue, 2) ; 

break; 
default : 

break; 

} 

} 

getxyz (lgtx, lgty, lgtz, vx, vy, vz) 

float *lgtx, *lgty, *lgtz, *vx, *vy, *vz; 

{ 

static char strl[] - "x- "; 

static char str2[] • "y- 

static char str3[] - "z- 

static char str4[] - "x« 

static char str5[] - "y m 

static char str6(] - "z- " ; 

int length; 

setvwpv () ; 

new_f rame-0 ; 

create^ temporary_segment 0 ; 
move_abs_2 (300 . , 700 . ) ; 

text ("enter the light source position [default : 0 . 0, 1 . 0, 1 . 0] "} ; 
move -- abs_2 (300 . , 680 . ) ; 
text (strl) ; 

set — echo_ position (KEYBOARD, 1, 0.31, 0.665); 


await_keyboard(1000000000, 1, string, ^length) ; 
if (length — 0) 

*lgtx - 0.0; 

else 

*lgtx » atof (string) ; 
m ° ve _abs_2 (300 . , 660 . ) ; 
text (str2) ; 

set_echo_position (KEYBOARD, 1, 0.31, 0.645); 
await_keyboard (1000000000, 1, string, ^length) ; 

if (length — 0) 

*lgty - 0.0; 

else 

*lgty - atof (string) ; 
move_abs_2 (300 . , 640 . ) ; 
text (str3) ; 

set_ echo^jjosition (KEYBOARD, 1, 0.31, 0.625); 
await_keyboard (1000000000, 1, string, ^length); 
if (length ■■ 0) 

*lgtz - 0.0; 

else 

*lgtz - atof (string) ; 

/* default light direction * / 

if ( (*lgtx — 0.0) && (*lgty — 0.0) (*lgtz — 0.0)) { 

*lgtx - 0.0; 

*lgty * 1.0; 

*lgtz * 1.0; 

} 

normlig (*lgtx, *lgty, *lgtz) ; 
move_ abs__2 (300 . , 600.); 

text ("enter the viewer position [default : -1000, 2200, 5000]") 
move^abs^ (300 . , 580.); 
text (str4) ; 

set_echo_position (KEYBOARD, 1, 0.31, 0.565); 
a wa i t_keyboa rd (1000000000, 1, string, & length ) ; 
if (length — 0) 

*vx -0; 

else 

*vx - atof (string) ; 
move_abs_2 (300 . , 560 . ) ; 
text (str5) ; 

set_echoj>osition (KEYBOARD, 1, 0.31, 0.548); 
await_keyboard (1000000000, 1, string, ^length); 

if (length -- 0) 

*vy - 0.0; 

else 

*vy • atof (st ring) ; 
move_ abs_2 (300 . , 540 . ) ; 
text <str6) ; 

3 et — echo_position (KEYBOARD, 1, 0.31, 0.527); 
await Jceyboard (1000000000, 1, string, ^length); 



if (length -» 0) 

*vz - 0.0; 

else 

*vz - atof (string) ; 

/* default view direction */ 

if ((*vx — 0.0) 4& ( *vy »- 0.0) 44 (*vz -- 0.0)) { 

*vx * -1000.0; 

*vy ■ 2200.0; 

*vz - 5000.0; 

) 

close_temporary_segment () ; 

} 

start_jup_core () 

{ 

int i# 

float X, y r z; 

if ( initial ize_core (DYNAMICC, SYNCHRONOUS, THREED) ) 
exit (1) ; 

if (initialize_view_surface (fivwsurf , TRUE) ) 
exit (2) ; 

if (select_view_surface (4vwsurf ) ) 
exit (3) ; 

initialize_device (KEYBOARD, 1) ; 
set^echo (KEYBOARD, 1, 1) ; 
set_echo_surf ace (KEYBOARD , 1, &vwsurf ) ; 
set_keyboard(l, 80, 1); 

set_font (1) ; 

) 

shut_down_core ( ) 

{ 

bred[0] * 0.8; 
bgrn[0] * 0.8; 
bblu[0] - 0.8; 

def ine_color_indices (fivwsurf , 0, 255, bred, bgrn, bblu) 
terminate_de vice (KEYBOARD, 1); 
deselect_view_surface (fivwsurf ) ; 
terminate_view_surface (fcvwsurf ) ; 
terminate_core () ; 

} 

int 

getob jdat (filename) 


char 

* filename; 

int 

if jf *; 

short 

vtznp, vl, v2, v3; 

float 

ftmp, maxd, scale, offset [3]; 

float 

x, y, z, x0, yO, zO, length; 


FILE 


*fptr ; 


if ( (fptr - f open (filename, "r")) •* NULL) ( 
printf ("Can' t open file: %s\n", filename); 
return (1) ; 

) 

fscanf (fptr, "%d%d", invert, inpoly); 

if ( (nvert > MAXVERT) II (npoly > MAXPOLY) ) { 

printf ("Too many object vertices or polygonsXn" ) ; 
return (2) ; 

fscanf (fptr, "%f %f %f %f %f %f ", ibbox[0][0], tbbox[0][l], tbbox(l](0], 
ibboxfl] [1] , t bbox [2] [0] , ibbox[2][l]); 
maxd « 0.0; 

for (i “ 0; i < 3; i++) ( 

offsetfi] - (bbox[i] [0] +bbox[i][l]) / 2.0; 
bbox[i] [0] — offset!!]; 
bbox [i] [1] — offset l i] ; 
if <bbox(i] [0] > bbox [ i ] [ 1 ] ) { 

ftmp - bbox [i] [0] ; 
bbox[i] [0] - bbox(i]tl]; 
bbox[i]ll] - ftmp; 

> 

if (maxd < bbox[i][l]) 
maxd - bbox[i] [1]; 

} 

scale - 1000.0 / maxd; 
for (i “ 0; i < 3; i++) ( 
bbox[i][0] *■ scale; 
bbox[i](l] *• scale; 

) 

for (i - 0; i < nvert; i++) ( 

fscanf (fptr, "%f%f%f", ivertices [i] [0] , tvertices [i] [1] , 
ivertices [i] [2] ) ; 

vertices (i] [0] “ (vertices ti] [0] - offset[0]) * scale; 

vertices [i] [1] - (vertices [ i] [1] - offset[l]) * scale; 

vertices [i] (2] - (vertices [i] [2] - offset [2]) * scale; 

normal(i](0] ■ 0.0; 

normal[i][l] “ 0.0; 

normal[i][2] “ 0.0; 

normalcount I i] *0; 

) 

k - 0; 

for (i - 0; i < npoly; i++) ( 

fscanf (fptr, "%d", inpvert[i]); 
if <(k + npvert(il) > MAXPVERT)- ( 

printf ("Too many polygon vertices\n") ; 
return (3) ; 

) 

pvertptr [ i] - ipvert(k]; 
for (j - 0; j < npvert (i] ; j++) ( 

fscanf (fptr, "%hd", ivtmp) ; 



pvert [k++] “ vtmp - 1; 

} 

planeq[i] [0] - planeq[i][l] - planeqti] [2] * planeq(i][3] - 0.0 

vl * pvert [k - 1] ; 

v2 ” pvert [k - 2] ; 

v3 ” pvert (k - 3] ; 

for (j ■ 0; j < 3; j++) ( 

planeqti] [0] +” -vertices [vl] [1] * 

(vertices [v2] [2] - vertices [v3] [2] ) ; 
planeq[i] [1] +■ -vertices [vl] [0] * 

(vertices [v3] [2] - vertices [v2] [2] ) ; 
planeq[i] [2] +■ -vertices [vl] [0] * 

(vertices [v2] [1] - vertices [v3] [1] ) ; 
planeq[i] [3] +- -vertices [vl] [0] * 

( (vertices [v3] [1] * vertices [v2] [2] ) - 
(vertices [v2] [1] * vertices [v3] [2] )) ; 
vtmp - vl; 
vl - v2; 
v2 - v3; 
v3 “ vtmp; 

) 

/* The normal of this face points to the center of the object. */ 

/* if (planeqti] [3] > 0.0) 

for (j - 0; j < 3; j++> 
planeqti] [j] m -planeqti] [j] ; */ 

x - planeqti] [0] ; 
y - planeq[i] [1] ; 
z - planeqti] [2] ; 

length - sqrt(x *x + y*y + z*z); 
planeqti] [0] /■ length; 
planeqti] [1] /“ length; 
planeqti] [2] /- length; 

planeqti] [3] /- length; 

if (debug) { 

printf ("planeq[%d] : %f , %f , %f , %f \n", i f planeqti] [0] , 

planeqti] [1], planeqti] 12] , planeqti] [3] ) ; 

for (j - 1; j <- npvert [i J ; j++) t /* accum normls */ 
vtmp ” pvert [k - j]; 
normal [vtmp] [0] +- planeqti] 1 0 ] ; 
normal [vtnp] [1] +- planeqti] [1] ; 
normal [vtnq>] [2] +- planeqti] [2] ; 
normalcount [vtnp] ++; 

) 

) 

for <i - 0; i < nvert; i++) t 

normal [i] [0] /« normalcount [i] ; 
normal [i] [ 1] /- normalcount ti] ; 
normal [i] [2] /- normalcount [i] ; 


if (debug) { 

printf ("normal [%d] : %f , %f , %f \n", i, normal ( i ][ 0 ] , 

normal [i] [ 1] , normal [i)[2] ) ; 

) 

) 

fclose (fptr) ; 
return (0) ; 


setvwpo(vx, vy, vz) 

float vx, vy, vz; 

{ 

int i; 

float diag, del, objdist, near; 


set_view_reference_point (vx, vy, vz); 
set_view_plane_normal (-vx, -vy, -vz); 
set_projection (PERSPECTIVE, 0., 0., 0.); 
set_view_plane_distance (256 . 0) ; 
if ( (vx -- 0.0) && (vz ■■ 0.0)) 
set_view_up_3 (0 . 0, 0.0, vy) ; 

else 

set_view_up_3 (0.0, 1.0, 0.0); 
set_window (-80 . 0, 80.0, -80.0, 80.0); 
diag - 0.0; 

for (i ■ 0; i < 3; i++) { 

del ■ bbox(i] [1] - bbox[i](0]; 
diag +■ del * del; 

) 

diag ■ aqrt (diag) / 2.0; 

objdist - sqrt(vx * vx + vy * vy + vz * vz) ; 

near - (diag >- objdist) ? objdist / 2.0 : objdist - diag; 

set_view_depth (near, objdist + diag) ; 

i* 

* set_window_clipping (FALSE) ; set_f ront_plane_clipping (FALSE) 

* set_backjplane_clipping (FALSE) ; 

*/ 

aet_viewport_3 ( . 125, . 874, 0.,.749, 0.0, 1.0); 

) 

static float invxform(4] [4] ; 


drawobj () 
( 


int 


i; 

float 


x. 

char 


ch 

bred [0] 

- 0.6; 


bgrn [0] 

- 0.0; 


bblu [0] 

- 0.6; 


cred[0] 

- 0.6; 


cgrn[0] 

o 

o 

R 



y, z, xO, yO, zO, length; 



cblu[0] - 0.6; 
if (renderhue ■- 0) ( 

set_text_index ( 1 ) ; 

def ine_color_indices (ivwsurf , 0, 255, bred, bgrn, bblu) ; 

} else { 

set_J:ext_index(l) ; 

def ine_color_indices (&vwsurf , 0, 255, cred, cgrn, cblu) ; 

} 

/* if (renderstyle 6S renderstyle < 3) */ 
for (i - 0; i < npoly; i++) ( 

cpindex[i] ■ 

(ambient + diffuse * (planeqfi] [0] * nix + planeqfi] [1] * 
nly + planeq[i] [2] * nlz) ) * 254.; 

} 

for (i * 0; i < nvert; i++) ( 

cindex(i] ■ 

(normal (i] [0] * nix + normal(i][l] * nly + 

normal[i][2] * nlz) * 254.; 
if (cindex(ij < 4) 
cindex[i] - 4; 
if (cindex[i] > 248) 
cindexti] “ 

248; 

) 

/* inquire_inverse_composite_matrix (invxform) ; */ 

/* if (renderstyle -■ 3) 

set_zbuffer_cut (fivwsurf , xcut, zcut, 2); 

else 

set_zbuf fer_cut (fivwsurf , xcut, zcutO, 2); */ 
for (i - 0; i < npoly; i++> ( 

/* if ( (visible (planeqli] ) ) II (renderstyle -- 0)) */ 
drawface (i) ; 


) 

) 

void 

normlig (lgtx, lgty, lgtz) 

float lgtx, lgty, lgtz; 

( 

float x0, yO, zO, length; 

/* map ndc_to_world_3 (lgtx, lgty, lgtz, £nlx, £nly, &nlz) ; 
map_ndc_to_world_3 (0 . 0, 0.0, 0.0, fixO, &y0, £z0) ; 

printf ( M lgtx - %f : nix - %f\n", lgtx, nix - xO) ; 

printf ("lgty - %f : nly - %f\n", lgty, nly - yO) ; 

printf ("lgtz - %f : nlz - %f\n", lgtz, nlz -z0); */ 

xO - 0.0; 
yO - 0.0; 
zO - 0.0; 


lgt x — xO; 
lgty -• yO; 
lgt 2 — zO ; 

length - sqrt(lgtx * lgtx + lgty * lgty + lgtz * lgtz); 
if (length !« 0.0) { 

nix ” lgtx / length; 
nly - lgty / length; 
nlz “ lgtz / length; 

) 

) 

int 

visible (pin) 


float 

pin [ ] ; 

int 

i; 

float 

c; 

0 

H 

o 

o 

for (i - 0; 

i < 4; i++) 


c +* invxform[2) (i) * pln[i]; 
return (c > 0.0); 

) 

drawface(p) int 
p; 

int i, j, k; 

short *ptr; 

float para ” 215.0 / 255.0; 

if (debug) { 

printf ("face : %d\n", p) ; 

) 

ptr - pvertptr [p] ; 
for (i - 0; i < npvert [p] ; i++) ( 

j - *ptr++; 

xlist t i ] - vertices I j] [0] ; 
ylist [i] - vertices [j] [1] ; 
zlist [i] - vertices [j] [2] ; 
if (renderstyle < 3) { 

indxlist [i] - cindextjl; 
if (debug) { 

printf ("vertices ( %d] : <%f,%f,%f> : indxlist - %d\n", 

i, xlist f i] , ylist ti], zlist [i], indxlist [i]) 

) 

) else ( 

dxlist [i] - normal tj] (0) ; 

dylistti] ■ normal [j] [1] ; /* Phong */ 

dzlist [i] - normal (j) [2] ; 

) 


) 


if (renderstyle < 2) { 

j - (cpindex [p] > 0) ? cpindex(p) : 1; 
j - (j > 255) ? 255 : j; 

/* j - (int) (para * (float) j) + 40; */ 

if (debug) ( 

printf ("indxlist [%d] : %d\n", p, j); 

) 

if (renderhue) { 

j . j » 2; t* if (j >62) j - 62; other colars */ 
if (! autodraw) { 
if (debug) { 

scanf("%d", 4j); 

) else 
do { 

awa it_keyboard ( 0 , 1, string, 4i) ; 

) 

while (i ■* 0); 

> 

set_fill_index( j + renderhue * 64 - 63); 

) else { 

if <! autodraw) { 
if (debug) ( 

scanf ("%d", 4j); 

) else 
do { 

await_keyboard (0, 1, string, 4i) ; 

) 

while (i ■■ 0) ; 

) 

set_fill_index( j) ; /* gray colar */ 

) 

) else if (renderstyle “2) { 

set_vertex_indices (indxlist, npvert [p] ) ; 

) 

/* Gouraud */ 
else { 

aet^ertex^normals (dxlist, dylist, dzlist, npvert [p]); 

} /* Phong */ 

if (renderstyle — 0) 

polyline_abs_3(xlist, ylist, zlist, npvertlpl); /* wireframe / 

else 

polygon abs__3 (xlist, ylist, zlist, npvert [p] ) ; 

) 

static float maxvw, vwpp, maxvdim; 
static float vlx, vby, vfz, vdx, vdy, vdz; 

static float minleft, maxright, minbot, maxtop, minback, maxfront; 
setvwpvO 

set view reference j oint (0.0# 0.0, 0.0); 

set viewjlane ^normal {0.0, 0.0, -1.0); 


set_view_plane_distance (0.0) ; 
set_project ion (PARALLEL# 0.0# 0.0# 1.0); 
set_view_up_3 (0 . 0# 1.0# 0.0); 
set_window (0 . 0# 1023.0# 0.0# 767.0); 
set_ view__depth (0.0# 1.0); 
set_> wind°w_( clipping (FALSE) ; 
set__viewport_3 (0 . 0# 1.# 0.0 # .75# 0.0# 1.) 
set_text_index (255) ; 

) 

void 
blkscr () 

{ 


int 


i» j; 



static 

float 

x [4] - 

i 




O 

CD 

1 

«» 

o 

CD 

o 

CD 

1 

o 

CD 

); 

static 

float 

y [4] - 

{ 




O 

CD 

* 

O 

00 

o 

CD 

O 

CO 

J ; 

static 

float 

z [4] - 

( 




o 

o 

o 

, 0. 


} ; 

static 

int 

n « 4; 




for (j - 0; j < 900; j++) { 

for (i - 0; i < 1152; i - i++) { 
pr_put (screen# i# j, 1); 

) 

} 

set_line_index (1) ; 
set_f ill_index (1) ; 
set_linewidth (0 . ) ; 
polygon_abs_2 (x# y# n) ; 


) 


APPENDIX D: 


C PROGRAMS FOR IMAGE PROCESSING USING 


THE 3M VDL VISION SYSTEM 


/* r ROC R AM FT C - This program calculates the Fourier Transformati 
of the input image. 

Vrittern by Chi rs K Vu, Jun, 1987. 

* / 


^include 
iinclude 
^include 
# i n c I u d e 


<stdioh> 

< ma t h . h > 
(string s.h> 

< c i h > 


^define SIZE 6 4 
^define M 6 


extern vdlinit(),vdlterm(); 
extern long sci < ) ; 

ma i n < a c , a v > 
i n t ic i char * av C 3 ; 

( 

FILE *fopen(),*fpl,*fp2,*fp3,*fp4,*fp5; 

char cmd [ 5 0 3 , fi!ereC503, fiIeimC503; 

int u 1 x 1 , u l y i ; 

int i,j,ti # tj,nt; 

float rebufCSIZEl , i mb ufCSIZEl ; 


i f ( ac 1 * 2 ) ( 

pr int f ( " Usage: %s filename\n",avCO]>; 

printf<"This program calculates the FT of the i npu t \ n" ) ; 
p r i n t f < " i ma g e . The results are stored in the fol lowing\n" ) ; 
pr int f ( " f i 1 es : \ n" ) ; 

p r i n t f ( " r e_f i I ename - the real part of the FT\n">; 
p r i n t f < " i ro_f i 1 ename - the imaginary part of the FT\n M ); 
ex i t ( - l ) ; 

} 

if ( v d 1 i n i t ( > < 0) C 

pr intf ( "umble to initialize VDL\n”>; 
e x i t ( - 1 ) , 

> 

sc i < p i np , (C00RD>1 *0) , 
sci (pcop, < COORD ) 2 , 0 , 1 ) ; 
sci(pbuf; (COORD) 2/ -2, 1) ; 
sci ( pw i n , (COORD) 1 , - 1 ) ; 
s c i ( pwi n, (COORD) 1 , 0) ; 

5 c i (pbuf t <COORD)2,-2,0) ; 

printf ("Due to the stack size limitation, the size of the image to 
pr int f ("processed is restricted to 64x64. \n\n">; 

.printf ("Enter the coordinates of the upper left corner of picture 
scanf < "%d %d" , 6u 1 x i , &u I y i ) , 
f p 1 3 f o pen ( 11 f t temp .re", M w" ) ; 
fp2»fopen( M ft temp . im" , ff w M > ; 
f pS-f open< "wu p i c ” , H w M ) ; 

/•perform FFT to each row*/ 

for (i=ulyl;(ti«i-ulyi><SIZE;i++> ( 

for (j*uIxl,(tj*j-uIxl)<SIZE;j++> { 
nt*ti+t j , 
imbufCtj]*0; 

rebuf C tj 3 « < f loat )pw(-l ,nt) 

*<float)(sci(pget,( COORD )3,(long>J,(Iong)i,i)+128); 


on 


b«\n" 

" ) i 



fwrite<rebuf,4,SIZE,£p5> , 

( f 12 ( rebuff imbuf ( SIZE ,M) ; 
fwri ttUebuf , 4 , SI ZE , fpl> ; 
fwri teC imbuf / 4 , S I Z E , fp2) , 

) 

f c I c 5 e < fpl) , fclosef fp2) ; f c 1 o s e < f p 5 > ; 

sprintf ■ cmd , " /dO /user / f lipreal ftt emp re %d M , SI ZE > ; 
system(cml) ; 

sp r i n ♦. f < cmd , " / d 0 / use r / f 1 i p r • a 1 f tt«mp . im Sd" .SIZE); 
syst em< end ) } 

/* perform F FT to each column */ 

fpl=fopenC H fttemp. re" x " ) ; 
f p 2 = f open ( " f t t emp im" , " r " > ; 

sprintf(filere,"r e_" * * 
strnciK f i lere , avC 1 1 , 10) ; 
f p3 = f open ( fi Iere ( "w") ; 
sprint ft fileim," i m_" > ; 
strncat < fi I e i m, av C 1 ] , 10) ; 
fp4sfopen( f i 1 e im, "w M ) ; 
for f isO, i^SIZE, i++) f 

f read ( rebuf , 4 , 5 1 Z E , fpl) ; 
f r e a d ( i mb u f ; 4 , 5 I Z E , fp2) , 
f f t. 2 < rebuf , i mb u f , S I Z E # M > ; 
fwri t e < rebuf , 4 , S I Z E , fp3) ; 
fwriteti mb uf , 4 , S I ZE , fp4> ; 

) 

f-Io5e(f?l);fclo»e<fp2);fcIos«<fp3>;£close(fp4), 

sprintf(cmd,"/dO/us«r/£Iipr«»l %s %d " , t i I e r • , S I ZE ) ; 

5 y 3 t e m < cmd ) ; 

sprintf f cmd , " / dO/ustr/f liprtal *hs # /*d " , f i 1 e i m , S I ZE > ; 

s y s t em ( cmd ) ; 

uni i nk C 14 f t t emp r e " ) ; 

un 1 i n Jc < M f t t emp . i m" ) ; 

if ( vdl term( M 0 ) C 

printf ("unable tp terminate VDL\n">; 
exit (-1 ) i 
) 

) 

pw ( x , n ) 
i n t x , n , 
i 

int i , p * 1 ; 

..if (nrsO) return(l) ; 
else ( 

for < i »0 ; i < n ; i ++ ) p»p*«; 

return(p) ; 

> 

) 


0-4 



/•SUBROUTINE F FTZ - C o o 1 e y -Tu r k e y radix-2 Ff'i algorithm 

Chris K Vu 
* I 

lincludt Cstdioh) 

•include <raa t h h > 

(f tz<x,y,n,n) 

float * x , * y ; int n,m; 

C 

int i , ie, ia i 3 iki 1 i 1 1 i "1 ,n2 ! 

float c , s , x t , y t , 

double pi , wr C 5 l 2 3 , wi C 5 1 2 3 ; 

pi=6 . 28319 / (double)n; 

for( i«0; i<n, i + +) I 

wrC i ]acos(pi* (double) i > i 
wiCil=sin(pi*<double>i> ; 

> 

n2 =n ; 

for ( i«0 ; i <n; i++) C 

nl =n2 ; 
n 2 an 2 / 2 ; 
i e=*n / n 1 ; 
i a = 0 ; 

f or( j=0 , j <n2 , j++) < 

c»( float )wrC ill ; 
s a ( f 1 oa t ) wi C i a 3 ; 
i a * i a + i e ; 

forCkaj ;k<n;kak+nl> f 

Uk + n2 ; 
xt*xCk3-xtll; 

xCk3axCk3+xC13; 

y t = yCk3-y C 1 3 ; 
y Ck3*yCk3+y[ 1 3 ; 
xC 1 l = c*»t + i*yt ; 
yC 1 ]ac*yt-«*xt ; 

3 

) 

) 

/•digit reverse counter*/ 
j * 1 ; 

f or( i*l ; i <n; I 
if ( j > i > C 

1 a i - 1 ; 

1 1 « j - 1 ; 
xt-xC 113; 
xC 1 1 3axC 1 3 ; 
x C 1 3 * x t ; 
x t ay C 1 1 3 ; 
y C 1 1 3ay t 1 3 ; 
y C 1 3 a x t ; 

1 

kan / 2 ; 

while < i > k ) C 
j a i-k ; 
k-k/2 ; 

1 

j * j «-k ; 

} 

return<0) ; 

1 


/ # 

SUBROUTINE F L I PRE AL -Th i s routine flips the input floating point 

ras ter along the diagonal axis for the 2-D 
FFT algorit hm . 
written by Chris K Vu 

* / 


einilude <stdioh> 

ma i n ( a c , a v ) 
int a c , ch a r * a v t 3 , 

{ 

char t emph C 2 0 3 , cmd t 8 0 3 ; 

: n t size; 
if ( a c ! = 3 ) C 

print! (" Usage: •/•* raster siz«\n" , ivC03) ; 

ex i t ( - i > ; 

} 

size=atoi(iv:21)/2; 
sprintf ( temph/'f rt emp " ) , 
dorotateCavtll , t«mph, site) ; 
uni i nk (av[ i ] ) , 

sprintf (cmd/'renime f r temp •/•* " , a v C 1 3 > > 
ay s tem( end ) , 
e x i t < 0 > , 


t + 

SUBROUTINE DOROTATEO -- Hindi* rotation* of the raster file 
» l 

dorotatef raster , t tmph , size) 
char *raster,*t emp h ; int size; 

( 

int i , j ; 

float buflC323C443 ,buf2C323C323i t emp C 6 4 1 , d t ; 

FILE * f open < > , * f p , * f p 1 ; 
fp=fopen(raster,"r" ) ; 

for <i=0;i<si*e;i+*> C 

fread<buflCi3 , 4, size, fp) i 
fread(buf2C i 3 , 4 , s i X e , fp) i 
) 

f p 1 * f o p en ( " r o t emp " « "w" > ; 
f o r ( i * 0 ; i < * i * • ; i ♦ ♦ > 

dt>bufi nuiii 

buf 1 C i 1 C j 3»buf 1 C 3 3 C i 3 ; 

buf 1C j 3 C l 3«d t ; 

) 

for ( i = 0 , i < size; i ♦ ♦ > C 

fwrite<bufiCi3,4,*ize, fpl) ; 
fwrite(buf2ti3,4,size,fpl>; 

) 

for (i«0,i<size;i++> f 

f r ead( buf 1 C i 3 ,4 , s i z# , f p > ; 
fread(buf2Ci3,4,size,fp>; 

3 

fclose(fp); 
for(i*0,i;siz*« i ♦♦ ) 

fot C 

d t *bu f 2 C i 3 C j 3 ; 


bu £ 2 C i][j3-buf2Lj]Ci3 , 
buf2Cj3Ei3=dt, 

} 

for (i=0;i<si:e;i++) C 

fwr i te^ buf 1 C i M , si :e , f pl^ » 

fwr i t • < buf 2 C i 3 , A , si z« < tpi > * 

) 

f clo»«v f pi ) ; 
fp=fopen<" rot emp" , M r") ; 
for<i=Q,i<size;i++> C 

f reed ( temp , 4 , s i le , f p > ; 
f read< buf 2C i 3 , 4 , s i t e , f p > ; 

} 

f or ( i=0; i < s i ze ; i+ + > 

f o : ( 3 = 0 ; j ( s i ze ; j + + ) C 
dt=buflCi3£j3; 
buflCilC }]=buf2Cj]C i 3 ; 
buf2Cj3Ci3*dt; 

} 

f dost ( f p> ; 

fp=fopen("rctemp" r"> , 
fpl = fopen( temph; "w" ) ; 
for ( i«0, i (si it; f 

fread( temp , 8 , size, f?> ; 
fwrite 1 ' temp, 4, size, fpl) ; 
fwr i te( buf 2 l i ] , 4, size, fpl) ; 

) 

for(i=D,i<size,i++> C 

f r e a d < temp, 8, size, fp) , 

for (j=0,:<size;j++) tempCj3=buflCi3Cj3 
fwri te< temp, 9, size, fpl > ; 

) 

f c I o s e ( f p ) , 
fclose(fpl); 
un 1 i nk ( " r o temp " > ; 
return(O) ; 

) 


/ * PROGRAM SPECTRUM. C - This program tales the outputs from 
the FTC program and c r e a t s a data file containing the 
spectrum, of the image 

Vrittern by Chris K Vu, Jun, 198? 

* / 

^include <stdioh> 

^include (raath.h) 

#include<strings.h> 

^define SIZE 64 

ma in(ic, a v ) 

int ac;char *av[) ; 

{ 

char fiIereC503 # fileimC503; 

FILE f open ( ) , * f p 1 ,*fp2, * f p 3 ; 

int i , j ; 

long iadC5IZE3, 

float scale, £ ma x ; 

float rdCSIZEl , idCSIZE] , fadCSIZE] , 
double dad, 

;f (ac>=2) ( 

p r int f ( “Usage : %s filename \n" ,avC03 ) ; 

p r i n t f ( " f i I e n ame - o u t p u t data filename from the FT program. \n" ) , 
priatf ("The output data is stored in file s p_f i I enam* . \ n" ) ; 
ex i t ( -1 ) ; 

) 

f j* a x = ( - t 2 8 . ) , 
sprintf(filere," r e_" ) ; 
strncatC filere,avC13> ; 
sp r int f ( f i 1 e in, " i m_" ) ; 
strncat(fi 1 e im, a v C 1 ] ) , 
f p 1 = f open( fi iere, “r") ; 
fp2=£open* fi l * im, " r " ) ; 

fp3=fopen( , '5ptemp","w") , 

for(i= 0 ,i<SIZE,i++> C 

f read ( td , 4 , 3IZE , f p 1 ) , 
f r e a d ( i d , 4 , S I ZE , fp 2 ) ; 
f or < 3 =0 , j < SI 2E ; j ♦♦ > t 

dad=(double)(rdCj3*rdCj3+idCj]*idCj3>; 

f a d t j 3 * < float)sqrt(dad) , 
if<fadCj3> f ma x ) fmax* f ad [ 3 3 , 

) 

fwrite( fad, 4 ,SIZE,fp3> ; 

) 

1 dad=(double) fmai+1 ; 

scale*256 /(float)log(dad), 
f dost ( f pi ) ; fclo»e(fp2> , fclose<fp3> , 
sprintf(filere," s p_" > ; 
strncat(fiIere,avC13) j 
f pl*f open( "sptemp" , "r " ) ; 
fp2«fopenf filer#, "w" ) ; 
f or( i*0; i ( S I 2 E ; i ♦ ♦ ) £ 

fread<fad,4,SIZE, fpl> ; 
f or ( 3*0; j < S I Z E ; j ♦ ♦ ) ^ 

dad*(doub!e) fadC j 3 + 1 , 

iadC 3 3*( long) (scaled float ) log(dad) )-12Q, 


'write* i*d,4,SIZ£, fp2) , 
3 


uni 1 nk ( M *p t erap " ) , 
exit(G), 

) 



/* PROGRAM CREATPIC.C - This program plots the spectrum of the 
image on the monitor 
Vr i t tern t y Chris K Vu Jun, i ° 8 7 

* / 

^include ( c i h > 

^include <stdioh> 
define SIZE 6 4 

extern v d 1 i n i t ( ) ; 
extern long sci ( ) ; 
extern vdl tirn( ) ; 

nai n < a r g c , argv, envp) 
int argc; 

char *argv[ ] , *envp[ 3 ; 

r 

\ 

int i,j,ti,tj,xs,ys, 
long buf CSIZE] ; 

FILE *fopen(),*fp, 
if (argc ' = 2 ) ( 

pr int f ( "Usage: °/os s p_£ i 1 e r. a me ■ n " , a r g v Z 0 1 ) , 

ex i t ( - 1 ) ; 

) 

if (vdl;nat<> ( 0) C 

fpr int f ( stderr , " unable to initialize VDLVn" ) ; 
r e t u r n ( - i > , 

} 

x s * 1 2 1 ; 
y s = 1 2 i 

s c i ( pw : n , (COORD) 1, 0 > ; 
sci ( p c I e , < COCRD • 1 ,0) ; 
ici ( pbuf , <CCORD)2,-2,0> ; 
f p=f openc argvC i ] , n r M ) ; 

for( ] sy s ; < t j = < j - y s ) / 2 > < S I Z E , j = j + 2 ) { 

£read(buf,4,5IZZ, fp) , 

for ( isxs , (tis(i-xsW2)(SIZE; i = i + 2 > C 
i f ( bu f C t i 3 > 1 2 7 ) b u f C t i 3 « 1 2 7 ; 
if(bufCtij<-128) bufC ti3»(-128> ; 
sci (pput , < COORD > 4 , ( long) i # ( long) j , b u f C til , 1) ; 
sc i ( pput , ( COORD ) 4 , (icng)i , <long>j+l,bufCti],l) ; 
sci (pput , (COORD) 4 , < long) i + 1 , ( long) j , buf C t i 1 , 1 ) ; 
sci (pput , < COORD >4, ( long ) i + 1 , ( Iong)j+l,bufCti3,l> ; 

} 

) 

f c 1 o s • < fp) ; 

sci(pbu£,<COORD)2,-2,0> f 
sc i < pwi n , ( COORD) 4,100,100,263,263); 
sci (pout, (COO RDM ,0) , 
if (vdlterm( ) < 0) ( 

fprintf ( stderr , "unable to terminate VDL\n">; 
return(-i) , 


) 


/* PROGRAM PATTERN C - This program identify a 2-D geometry from 
a library of geometries. The geometric properties, area, 
perimeter, and maximum radius, of the viewed geometry are 
measured and used to check against the library set 

Writ tern by Chris K. V u , May, 1987. 

* / 

# include < s td i o . h> 

^include <math h> 

^include <ci.h> 

extern vdlinit(),vdlterm(>; 
extern long sci ( ) ; 

i n t w i nd ow[ 4 ] , threshold, count ; 
i n t are a_a rrCIO] , peri^irrC 101 , mr_a rrC 101 ; 
char name.i rrClOlCSO], 

ma i n ( ) 

( 

long v ; 

int area, peri, mr , c ; 
int cx,cy,mx,my,t imp ; 
double distsq; 

if (vdlinitO < 0 > C 

print f ("unablt to initialize VDL* 1 ); 
ex i t ( - 1 ) ; 

} 

printf C'pUct any object inside view area to " ) ; 
printf ("deterisine the view window\n u ) ; 
s c i C p d i g , ( COORD ) l , 1 ) ; 
s c i ( pw in, (COORD) 1,-1); 

printf ("perform learning process <y/n)? " ) ; 

while ( ( c - getcharO) !* 1 y * && c i- ’Y' && c !* ‘n’ && c ■= 'N* 

; 

ge t cha r < ) ; 

if ( c •• 1 y 1 ii c *■ 1 Y 1 ) 

1 ea rn^sub < ) ; 

el se 

g e t d a t a < ) , 

printf ("pUce the unknown geometry inside the view window\n"); 
start : 

5 c i < pd i g , ( COORD >1,1); 

sci (pthr , ( COORD > 1 , ( long) threshold) ; 

sci ( p f r q , ( COORD ) 0 ) ; 

sci (prng , ( COORD ) 2 , - 1 2 8 , 0 ) , 

v » sci (prpc i < COORD ) 0 ) ; 

area* < i nt > v ; 

pr int f ( H the area of the viewed geometry is %d u , eree) ; 

p r i n t f ( M pixels \n\n H > ; 

sci(pbin,<COORD)0); 

sci < p f r q , < COORD > 0 ) ; 

sci (prng, ( COORD ) 2 ,0, 127) ; 

v ■ sci (pcen, (COORD)O) ; 

cy = (int)v / 512, 

cx * (int)v - ( c y * 5 t 2 > ; 

v * s c i ( p r pc , ( COORD ) 0 ) ; 



peri * ( int ) v; 

pr intf (" the perimeter of the viewed geometry is %d" ( per i ) , 

p r i n t f ( " piKtIi\n\n") i 

v a sc i ( pmnx ,< COORD > 0 ) ; 

m-y ■ ( i n t ) v / 5 1 2 ; 

mx a ( int ) v - my*5l2; 

distsq » ( c x -mx ) * ( c x -mx ) ( cy-ray ) * ( cy-my ) ; 

me = <ir*t> sq r t < d i s t sq ) ; 

v = s c i ( pmny , ( COORD > 0 ) , 

my » ( int ) v/512 ; 

mx * ( int ) v - my * 5 1 2 ; 

distsq « ( c x -mx Y* < c x -mx ) + < cy-my ) * (cy-my ) ; 

temp « (int) s q r t ( d i s t s q ) ; 

if (temp > mr ) mr«temp; 

v a s c i ( pmx x ,( COORD ) 0 ) ; 

my a <int)v/512; 

mx a (int)v - my*512; 

distsq * ( c x-mx ) * ( c x -mx ) + ( c y -my ) * ( c y -my ) ; 

temp a (int) sqrt(distsq) , 

if (temp > mr ) mratemp; 

v a sci (pmxy , (COORD) 0) ; 

my a ( int) v/512; 

mx a (int)v - my*512; 

distsq a ( c x -mx ) * < c x -mx ) + (cy-my ) * ( cy-my ) ; 

temp a (int) gqrt(distsq) ; 
if (temp > mr ) mr « t emp ; 

pr int f ( M the maximum radius of the viewed geometry is %d M ,mr); 
printff" pixels\n\n u ); 
recog(area,peri , mr > , 

p r i n t f ( " i d e n t i f y next unknown geometry (y/n>? 11 ) ; 

while ( < c a ge t cha r ( ) ) !« 'y' && c !« * Y 1 && c !» 'n' && c !a 4 N‘ ) 
/ 

get cha r ( > ; 

if ( c aa ' y 1 : i cam 1 Y 1 ) 

goto start; 

e I se 

if ( vdl term( ) < 0) C 

printf ("unable to terminate VDLVn"); 
ex i t ( -1 ) ; 

) 

} 

/* subroutine to learn the geometric properties of 
the target objects */ 

lea rn_sub < ) 

{ 

int i,c,cic,cy,mx,my,t emp ; 
long v ; 

double distsq; 

print! ("uie the numerical key pad to move the " > ; 
printff" cursor to determint the threshold vilu«\n") ; 
p r i n t f ( " press (return) key to •*it\n" ); 
sci (pers, (COORD)O) ; 

print! ( "enter the threshold value "), 
scanf ( "Sd" ^threshold) ; 

printf< "threshold value is set to %d\n" , threshold) ; 
count a i ; 


loop 

printf <"plice the geometry inside the view window and\n M ) , 
pr i nt f ( "input the given name of the viewed g e ome t r y \ n " ) , 
scant ( M %s M ,na me _arrCcountl) , 
sci ( p d i g , < COORD ) 1 , 1 ) ; 

sex (pthr , (COORD) l , ( long) threshold) ; 

sci ( p f r q , (COORD)O) ; 

sc i ( prng , (COORD)2, -123,0) ; 

v = sc i ( p r pc , ( COORD ) 0 ) i 

a r e a_ a r r C c o un t 3 a (int) v; 

printf ( M th* area of the viewed geometry is %d " , a r e a — a r r [ c o un t ] ) ; 

p r i n t f ( " pixels\n\n" >; 

sci ( p b i n , (COORD)O) ; 

sci (pf rq , (COORD)O) ; 

sci (prng, ( COORD ) 2 , 0 , 1 2 7 ) ; 

v * sc i ( pcen ,< COORD ) 0 ) ; 

cy a (int)v / 512; 

cx = (int)v - ( c y * 5 1 2 ) ; 

v * s c i ( p r p c ,( COORD ) 0 > ; 

peri^arrCcount] « (int) v; 

print! ("the perimeter of the viewed geometry is %d " , p e r i _a r r C c o un t 3 ) ; 
pr int f ( " pixels\n\n" ) ; 
v « s c i ( pmnx ,( COORD ) 0 ) ; 
my a (int)v/512, 
mx a (int)v - my * 5 1 2 ; 

distsq a ( c x -mx ) * ( c x -mx ) + ( c y -my ) * ( c y -my ) , 

mr_a r r C count ] a (int) sqrt(distsq); 
v a s c i ( pmny ,( COORD ) 0 ) , 
my a<int)v/512; 
mx = (int)v - my*512; 

distsq a ( c x -mx ) * ( c x -mx ) + ( c y -my ) * ( cy -my > ; 

temp a (int) iqrt(distsq) ; 

if (temp > mr_a r r C count ] ) ror_a r r C c o un t ] a t emp ; 
v a 5 c i ( pmx x , ( COORD ) 0 > ; 
my a ( int ) v / 5 1 2 ; 
mx a (int)v - roy*512; 

distsq = (cH-mx)Mcx-nix) + ( c y -my ) * ( c y -my ) ; 

temp a (int) sqrt (distsq) , 

if (temp > mr^arr (count ] ) mr_a r r C co un t 3 a t emp ; . 
v a sci ( pmxy , ( COORD) 0 > ; 
my a < int ) v / 5 1 2 ; 
mx a (int)v - my*512; 

distsq a ( cx— mx ) * ( cx— mx ) ♦ ( c y -my ) * < c y -my ) ; 

temp a (int) sqrt(diitsq) ; 

if (temp > mr^trrCcount] ) mr_ a r r C count ] a t emp ; 

print! ("the maximum radius of the viewed geometry is " , mr_a r r C c o un t 3 ) ; 

pr int f ( " pixels \n\n" ) ; 

print! ("next geometry (y/n)? ">; 

while < (c « getchxr ( ) ) !* 'y' && c !a 'Y' && c !a 'n' && c !■ ’N' ) 

i 

getchar ( ) ; 

if ( C a » 'y* I! caa * Y ' ) { 

++count ; 
goto loop; 

) 

else 

savedata< ) ; 
return(O) ; 

) 



/ * subroutine to save data to file "pattern d a t " * / 

saveda ta ( ) 

C 

i n t i ; 

FILE * f open ( > , * f p , 

fp = fopenC" /dO/user/pat tern.dat" ,"w" > ; 
fprintflfp/'^dXn", threshold) ; 
fpr int f ( fp, "%d\n" , count ) ; 
for < i s 1 t i < count 4- l t i + + ) ( 

fprintf<fp, 1,| As\n" , narae^a r r C i 3 ) ; 

f p r i n t f < f p , ‘ ,f /od B /*d %d\n u , area^arrC i 1 , per i^ar r C i 3 ,mr_arrC i ] > ; 
) 

f c I o s e ( fp) ; 

} 

/* subroutine to read the geometric properties of the 
learned objects from the data file * / 

g e t d a t a ( > 

{ 

i n t i ; 

FILE *fopen<),*fp, 

fp = f o p e n ( " I d 0 / u 5 e r / p a t t e r n . d a t " , “ r 11 ) , 
f scan f < f p , "•Ad" ^threshold) ; 
f scanf ( f p , '"/•d" ,&count) ; 
for ( i»i ; i <count4-l , ;+♦) { 

f scan f < f p , "Vos" , name_a r r C i 1 ) ; 

f s c a n f < f p , "Sd ®^d %d" , & arc a_a rrC i 3 , &ptr i_a r r C i 3 , Siar^a r r C i 3 ) ; 
} 

fclose(fp) ; 

) 

/* subroutine to identify the viewed geometry against 
the library geometry */ 

recog ( area, peri , mr ) 
int a r ea , pe r i , mr ; 

{ 

int i , c , rec_f 1 ag»Q ; 

char * rec_name l , * r§c_na»t2 ; 

double cdsq , cd , ir r , ri c_e r r 1 * i , r e c_e r r 2 ; 

'for <i»i;i<count+l,i++> C 

cdsq * <double)area/<doubie)a rta_a rrC i 1 ; 
cd « sqrt(cdfq) ; 

err * 05 * fabs<peri-cd*peri_arrti3)/<cd*peri_arr[i3>; 

irr « err ♦ 0.5 * f a b s ( mr -mr_ a r r C i 3 * cd ) / ( mr_a r r t i 3 * cd ) ; 
print f( n match the viewed geometry to %s \ n " # niBt^a r r C i 3 ) i 
print If" %72f percent in difference\n\n",err*100.>; 
i f ( er r < 0 1 ) 

if ( r i c_ flag 1 * 0 ) C 

if (err < rec.ir r 1 ) < 

rec_er r2 * rec_errl ; 
rec.naraeZ * rec — namel; 
rec err 1 ■ err; 



name_a r r t i 1 


r e c_name 1 * 

) 

else C 

rec_er r2 - err; 
rec_name2 a name_arrCi3, 

> 

r ec_£ lag a 2 ; 

} 

else C 

r e c_e r r 1 a err; 
rec^namel * name_irr C i 1 ; 
r ec_f legal; 

) 

3 

if ( re c_£ lag ! = 0 ) f 

pr int f ( " \ n\ ntht geometry is recognized as \n\n *** * * * " * r e c__name l ) , 

printfC" with confidence of %5 . 1 f percent \n\n M ,( 1 . -rec_«r rl )* i 00 ) ; 

i f ( rec_f leg a= 2 ) 

pr int f ( " *** %s *** with confidence of l f percent\n\n” , 

r e c_ name 2 , ( 1 . - r e c_e r r 2 ) * 1 0 0 ) ; 

) 

else { 

printf ("can* t recognize the geometry. \n”) ; 

printf ("want to store the geometry into library (y/n)? " > ; 

while ( (c * getcharO) !a ‘ y' && c !* ' Y ' && c != ’ n' & & c • 1 N * ) 

i 

if < C as ’y' !! C sa ' Y 1 ) C 

count ++; 

printf ("input the given name of the g e ome t r y \ n" ) , 

scanf ("%s" f name_a rr [count 3 ) ; 

aree^arrCcountl = area; 

peri.arrfcount] a peri; 

mr^arrtcountl a mr ; 

savedataC ) ; 

3 

) 

return(O) ; 

3 



Tne macro 'XYANGLE.M' 



wi n( 50 ,70 ,350 ,450) 

divj( ) 

thr (-20) 

bin ( ) 

rng (110,127) 

v=mxx( ) 

ry s v/512 

rx=v-( ry*512) 

v=mxy( ) 

by = v/512 

bx=v-(by*512) 

rng( 110,127) 

v=cen( ) 


/* set up the processing window */ 

/* acqui re image */ 

/* convert into binary image */ 

/* binary edge */ 

/* specify the intensity range */ 

/* the rightmost pixel within the intensity range */ 


/* the bottommost pixel */ 


/* locate t he gravity center */ 


cy=v/512 

cx=v-(cy*512) 

print "the gravity center is located at\n " 

print "x = " 

print cx 

print " \ n" 

print "y =" 

print cy 

print "\n“ 

den=rx-bx /* compute the orientation angle */ 

num=by-ry 

print "tan(theta) = " ( 

print num 

print 

print den 

print “ \ n" 


Digital picture 
output from the Of 1 
system 

(Black hexagox bo> 
with 0 = 60°) 


Binary image 
resulted from the 
thresholding proce 
(tnr()) on picture 

a. 


Wireframe resulted 
from the binary ed 
detection process 
(bin()) on picture 




Figure 


Image space and the principal axes of area moment 
of inertia about the gravity center 
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Executive Summary 


Objectives : 

The objectives of this project have shifted from computer graphics and vision systems 
to the broader scope of applying concepts of artificial intelligence to robotics. 
Specifically, the research is directed toward developing artificial neural networks, 
expert systems and Laser Imaging techniques for autonomous space robots. 

The statement of work is 

1. Develop a Computer Graphics laser range finder simulator 

2. Use laser imaging simulator to study use of artificial neural 
networks for antonomous robotic navigation. 


Overview : 

Primarily as a result of our contacts with Dr. Timothy Cleghom through this RICIS 
project the emphasis of much of our robotics research has changed to research 
directed toward artificial intelligence applications in robotics. We have become 
interested in applications of CLIPS, NETS and Fuzzy Control. 

A laser range finder simulator has been developed and used to study use of artificial 
neural nets for robot navigation. CLIPS and NETS have worked their way into our 
robotics courses and we are currently investigating an application of NETS in aero- 
dynamics. During the past two years there have been 5 MS and 5 PhDs degrees 
awarded to our students in robotics. Three MS and three PhD theses are directly 
related to this project and are summarized in this report. Also a number of technical 
papers and reports that have been written during the duration of this project are listed 
at the end of this report. 
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I. Introduction 


This report describes recent robotic activities in Mechanical Engineering at Rice that 
have been stimulated by RICIS Research Activity AI.2. This work has been directed 
toward developing artificial neural networks, expert systems and Laser Imaging 
techniques for autonomous space robotics applications. A computer graphics laser 
range finder simulator developed by Wu has been used by Weiland and Norwood to 
study use of artificial neural networks for path planning and obstacle avoidance. 

The support provided by NASA/JSC Software Technology Branch through RICIS has 
been leveraged significantly through support of 5 graduate students by the U. S. Army 
(Captains Norwood, Weiland, Schuster, Atkinson and Hanusa). 

Theses research by Wu, Norwood, Weiland and Schuster is summarized below. An 
application of NETS in the advanced robotics lab course is described and the potential 
for NETS in aerodynamics is discussed. Our recent interest in fuzzy logic control for 
redundant manipulators and mobile robot navigation is mentioned. Finally a listing of 
recent technical papers, reports and theses is given. 

II. Laser Imaging System Simulator 

A computer graphics simulation was developed to permit emulation of laser range 
scanners for use in our robotics research. Details of the simulation are contained in 
the PhD thesis by Chris Wu (1990). A brief summary of his thesis is given in this 
section. 

The Laser Imaging Simulation Algorithm (LISA) provides simulated laser images 
using the system parameters entered by the user. All parameters (such as hardware 
specifications, output data size, and viewing parameters) are changeable to permit 
emulation of the ERIM and Odetics laser scanners as well as future scanners. LISA 
computes system characteristics such as signal-noise ratio, ambiguity interval and 
viewing volume and then generates appropriate images. 

Traditional computer graphics techniques that extract - Z-buffer data for distance 
information would give an orthogonal view of the scene whereas a laser scanner 
produces a radial measure. A straight line in the real world will appear curved in the 
spherical world. Thus the simulation program emulates this radial measure by 
emitting a ray corresponding to each rangel position. Each ray’s intersection with the 
image is calculated and returned as a proportion of the ambiguity interval. 

Reflectance values are calculated in the simulation program based on the object’s 
orientation relative to the scanner. LISA generates artificial noise to simulate sources 
of noise in real laser systems such as transmission noise, photon noise generated by 
the photodetector, ambient noise, and noise in subsequent amplifiers. 

After generating the range and reflectance data of a scene, the ambiguity interval 
effect is obtained by dividing each range value by the ambiguity interval and storing 
the remainder in the depth buffer. The depth and reflectance buffers are then 
combined and saved to a file for further processing. The output data are arranged in 
standard 16-bit format where the higher 8- bit stores the range data and the lower 8-bit 
stores its corresponding reflectance value. 

LISA has been used in robotics research by Wu, Norwood and Weiland. Their theses 
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are listed at the end of this report. 

III. Application of Laser Range Finder to Robotic Navigation 

This section summarizes research conducted by Peter Weiland (1989) for his MS 
thesis. The laser ranger finder simulator discussed above was used by Weiland to 
study the application of laser range data to autonomous robotic navigation. 

A scanning laser was chosen for the sensory task because of its ability to readily 
determine the spatial relationship of objects in the scene. It is not hampered by 
illumination or the correspondence problem of vision systems. It readily displays 
distance measurements to each point in the scene. In the navigation system proposed 
by Weiland a preprocessor is used to enhance the laser data. It first screens the laser 
data for missing points caused by specular reflection. It then scans the image for the 
existence of ambiguity intervals and noise. Missing points and noise are corrected by 
using a ’donut filter’. Ambiguity interval problems are corrected by the addition of 
one cycle to the range value. After preprocessing, the data are sent to the processor 
for conversion. 

Actual conversion to the Cartesian representation is done in the processor stage. 
Knowing the geometric configuration of the laser scanner, data points are converted to 
an (x,y) plane. A corresponding elevation, z value, is entered into the (x,y) plane at 
the appropriate position. The converted image requires additional procesing. There 
are numerous missing data points caused by the varying data point density across the 
image. It also has missing point regions caused by object occlusion. These missing 
points must be filled in by the post processor before the image is usable by the 
reasoning system. 

Conversion enhancement is accomplished by the post processor. It fills in missing 
points and frames the image. Shadow points are determined by finding the edges of 
obstacles from the laser image. Edges are determined by range discontinuities. 
Missing points caused by data density problems are filled in using the locus 
algorithm. 

Weiland’s research produced a perception process that effectively represents 
simulated laser data in a form that permits cognition. This process is the conversion 
of laser data into a top down, Cartesian elevation map of the sensed environment. 
System applicability was demonstrated by navigating through the artificial robot 
world using simulated laser data and the connectionist network navigation system 
developed by Norwood and described below. 

IV. Applications of Connectionist Networks 

Some applications of connectionist networks for robotic path planning, autonomous 
robot navigation, redundant manipulator control, autonomous insertion task and in 
aerodynamics are described in this section. 

IV -1. Path Planning and Obstacle Avoidance 

This subsection summarizes the MS Thesis research of Chris Schuster (1990). 
Automated path planning and obstacle avoidance have been the subject of intensive 
research in recent times. Many efforts in the field of semiautonomous mobile-robotic 
navigation involve using Artificial Intelligence search algorithms on a structured 
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environment to achieve either good or optimal paths. Other approaches, such as 
incorporating Artificial Neural Networks, have also been explored. By implementing 
a hybrid system using the parallel processing features of connectionist networks and 
simple localized search techniques, good paths can be generated using only low-level 
environmental sensory data. This system can negotiate structured two- and three- 
dimensional grid environments from a start position to a goal, while avoiding all 
obstacles. Major advantages of this method are that solution paths are good in a 
global sense and path planning can be accomplished in real time if the system is 
implemented in customized parallel-processing hardware. 

The electronic hybrid network system developed by Schuster represents an original 
method for constrained semi-autonomous robotic navigation control. It takes simple 
binary environmental input, along with a start and goal location, and processes the 
data through a connectionist network which provides a nodal ‘voltage potential’ look- 
up table. The voltage potentials are analyzed by a second network which determines 
the move direction based on an examination of the neighborhood around a given 
current node (beginning with the start position). Finally, the system presents a 
solution path based on a set of locally optimal steps. 

This system exhibits some distinct advantages over the traditional approaches. Due to 
the parallel architecture of the connectionist network, the system can be expected to 
be much faster (and possibly more damage resistant) than Artificial Intelligence 
search algorithms. This navigation system also has the flexibility to account for both 
moving obstacles and a moving goal. This is accomplished by simply applying new 
inputs to the connectionist network and reevaluating the path problem. Also, 
assuming the hybrid network system is implemented in hardware, it does not require a 
dedicated CPU/microcomputer to make it work and could conceivably be built right 
into the sensory system and servomotors of a robot. 

An example of a near term use for an expanded ‘Maze Machine’-type navigation 
system is the control of a robotic delivery vehicle in a large factory. The 
environment/floorplan would be fairly stable, however local sensors in the factory 
could easily update the ‘maze’ database on board the robot through radio 
communications. A human or central computer would assign the robot a task, 
probably also by radio link. An example task could be to "take pallet #196 from point 
A to point B". The navigation system on board the vehicle would then plan the path 
from the current location to point A (the pick up point for the pallet) and then plan a 
second path from A to B (the pallet drop off point). Other routines on board would 
handle the pallet upload/download procedure, emergency stop procedures, etc. 

There are disadvantages as well. First, the traditional AI procedure, regardless of 
whether they find an optimal path or just feasible paths, are proven methods that can 
be implemented relatively cheaply on microcomputers. The network system 
presented here can not be implemented on microcomputers while keeping its parallel 
architecture, however the sequential simulation AMAZ3D can be a valuable 
alternative method, especially in feasibility tests to determine if a custom VLSI 
network setup is warranted for a particular application. Since the AI approaches are 
software based, they can be modified much more easily to represent different type/size 
systems and environments. Also, because they are tried and proven procedures, 
software is readily available for their implementation. 

In summary, this electronic hybrid connectionist network system solves path 
planning/obstacle avoidance problems for a grid-structured, two- or three- 
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dimensional environment. This navigation system provides good (often optimal) path 
solutions based on a collection of locally optimal steps/moves found using the output 
of the sensory analysis connectionist network. It operates using only low-level 
(binary) descriptions of the environment which can be provided by a variety of current 
and experimental sensory systems. The navigation machine would best be used in 
combination with a computer and global sensory input systems. Possible applications 
of this hybrid network system span from the home to industry to outer space. 

IV -2. Autonomous Robot Navigation 

This subsection is from the MS thesis of John Norwood (1989) and the Ph.D thesis of 
Peter Weiland (1991). This research relates to that of Weiland that was discussed 
earlier. 

Robotic path planning and obstacle avoidance have been the subjects of intensive 
research in recent years. Most solutions to these problems have been reached through 
the use of traditional Artificial Intelligence search techniques. However, these 
methods have proven inadequate when applied to highly unstructured or unknown 
environments. By using an Artificial Neural Network, one can get near optimal paths 
using only low level information about the scene. In this way, it is possible to 
navigate from a start position to a goal position while avoiding all obstacles. Major 
advantages of the method developed by Norwood are that the solution is very fast and 
does not rely on any a priori knowledge of the robot’s environment. The system was 
shown to be very effective for path generation when used in conjunction with the 
simulated Laser Imaging System. 

Norwood’s connectionist network model consists of a collection of individual nodes 
each having N-l inputs and two outputs. The first of these outputs is a transient 
output and connects only to the inputs of the other nodes in the network. The second 
output is the one of primary interests in the navigation problem. This output is an 
artificial potential that is induced on a given node by both the goal node and the 
obstacle nodes in the network, the former having a high positive potential and the 
latter having a smaller negative potential. Once inputs, representing the robot’s 
environment, are applied to the network, the output potentials emerge as the transient 
outputs are propagated through the network. Navigation is then a matter of local 
processing in the neighborhood of the robot until such time as the current node is 
equivalent to the goal node. 

This navigation scheme has been shown to be well suited to simulated unstructured 
environments where precise information about the terrain and obstacles therein is not 
generally available. This characteristic of the network is developed from the low- 
level representation of the environment on which the network operates. The network 
model can be applied wherever a top-down view can be easily generated from the data 
provided by the Laser Imaging System Algorithm (LISA). 

The navigation network is a general method. It has the flexibility to account for both 
a moving goal and moving obstacles. This is done simply in hardware by applying 
new inputs to the network. In software, one need only set a variable which is checked 
at certain intervals to allow the looping structure to change. In this way, the new 
environment can be evaluated, and the potentials changed to generate an appropriate 
new path. Additionally, the network should generalize rather easily to a higher 
dimension to allow one to generate obstacle free paths in a robot workspace. 
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Robotic navigation has been an area of intense research since the onset of mobile 
robot development. The usefulness of mobile robots ultimately reside in their ability 
to move and interact with the environment. Current approaches to robotic navigation 
are primarily based on simulating intelligent, human-like behavior through the 
intelligent system model processing cycle; sense, perceive, reason, act. Unlike these 
methods, Weiland’s PhD thesis (1991) presents a navigation system based on 
biological and behavioral principles which function in a stimulus-response manner. 
Using connectionist architectures, a relationship between stimulus and response is 
acquired through the learning of conceptual information pertaining to navigation. In 
this research, the mammalian visual system provides a guide for the processing of 
environmental stimulus. Simulated laser range data are processed in retinal patch size 
elements by a cellular neural network. This network is designed to detect obstacle 
existence for each path segment based on an invariant feature of range discontinuity. 
Obstacle information is then mapped in binary format, indicating the traversable state 
of the patch, to the system’s visual cortex. Response to this mapping is derived from 
a hierarchical structure of back error propagation neural networks in which each 
network has learned a particular navigational behavior, obstacle avoidance, wander, 
and goal seeking. Outputs from these networks indicate appropriate motor response 
for the environmental stimulus. 

A simulation was develolped to evaluate the performance of this system by having a 
simulated robot traverse an environment. The connectionist approach was verified 
through system display of human-like navigational behavior for the simulation’s 
environment Advantages of the neural network approach were also demonstrated by 
its processing speed and adaptability. Procedures are discussed for actual system 
implementation in which cycle times of under one second are completely feasible. 
Proposals for unsupervised learning of navigational responses for environmental 
stimulus are also made. Weiland’s research provides a foundation for our continuing 
study of the connectionist approach to the problem of autonomous robot navigation. 

IV -3. Redundant Manipulator Control 

Norwood’s PhD thesis (1990) is concerned with the application of Artificial Neural 
Networks to the problem of controlling a redundant manipulator. 

Redudancy in robots is very much an open research area in the field of robotics. As 
the tasks required of robots become more and more complex, the ability of robots to 
perform satisfactorily in these applications must increase accordingly. Redundant 
manipulators have a greater ability to perform difficult tasks, such as obstacle 
avoidance, than non-redundant ones. In order to make use of this extra ability of 
redundant robots, more effective control schemes must continue to be developed and 
to this end, more and more researchers are looking to expand the body of knowledge 
in this area. Norwood’s thesis addresses the problem of moving a redundant robot 
within a defined workspace in the presence of obstacles. Additionally, criteria are 
developed that may be applied to the robot to constrain the redundant equations. 
Finally, a neural network solution to the redundant inverse kinematic problem is 
presented. It is shown that the inverse kinematics can be developed through a network 
architecture which provides accurate and fast solutions to a problem that is 
computationally and structurally complex. 

IV -4. Insertion Task 

The insertion task, described here, was performed by students in our Advanced 
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Robotics Laboratory course last semester. It illustrates how results of our research 
efforts are being transferred to the class room. Two-dimensional geometrical objects 
consisting of a triangle, a rectangle and a circle are positioned at random locations in 
the workspace of our RTX robot. The objective is to have the robotic system locate 
the object specified by a voice command, retrieve the object and insert it into the 
appropriate receptacle. 

This exercise involves use of the TI Speech System, the 3M VDL Vision System, 
NETS and a JR3 force/torque sensor. The perimeter and area data from the digitized 
images is sent to the artificial neural network that classifies the objects. A small 
circular rod located at the center of the area of the object to be retrieved, is grasped by 
the two parallel jaws of the gripper and force control is used during actual insertion. 

Currently this task is being extended to the use of our three fingered hand for grasping 
the objects by their edges rather than using the rod at the C. G. Then insertion will be 
attempted using finger control instead of arm motion only. This work is part of the 
current thesis research by Paul Hanusa and Joe Grisham. 

IV -5. Application of NETS in Aerodynamics 

Presently wind tunnel experiments and computational fluid dynamics (CFD) 
numerical modeling are the tools used by researchers and designers for understanding 
fluid flow phenomena and exploring various aerodynamic designs. For example, in 
order to predict the pressure distribution across an airfoil, one can use experimental 
data or obtain numerical solutions of the Navier-Stokes equations that depend on 
Mach number, Reynolds number, and angle of attack. 

We believe that Artificial Neural Networks (ANN) have potential applications in 
aerodynamics to complement and extend the traditional tools. As an example of such 
an application, four sets of experimental data for chordwise coefficient of pressure 
(Cp) distribution on the NACA 0012 airfoil were used to train a backpropagation 
network using NETS with Mach number as input and Cp at 38 points on the airfoil as 
output. The network was then tested at 5 other Mach numbers for which data were 
available with excellent agreement between the network output and the experiments. 
This work is being don e in c ollaboration with Dr. Andrew Meade, a member of our 
faculty who works with CFD research at Rice and had a NASA Faculty Fellowship at 
NASA Langley during the summer of 1991. 


V. Preliminary Studies of Fuzzy Logic Control Applications in Robotics 

After discussions with Mr. Bob Lea and Dr. Jani of the Software Technology Branch 
as NASA/JSC we have been convinced that there are opportunities for applying fuzzy 
logic control in robotics. William Atkinson and Kevin Magee are two graduate 
students in our program who are currently investigating ways to apply both artificial 
neural nets and fuzzy control to redundant manipulators. Sarmad Adnan, a PhD 
candidate, with the help of Alex Stewart, an undergraduate, is working toward 
developing a camera tracking system that can have applications for use with fuzzy 
control of a mobile robot. 


6 


Technical Papers, Reports, and Theses 


Cheatham, J. B., and Norwood, J. D., "Robotic Path Planning and Obstacle Avoidance: A Neural 
Network Approach", 5th LASTED Int. Conf. Expert Systems and Neural Networks, Honolulu, 
Aug. 16-18, 1989. 

Cheatham, J. B., Weiland, P. L. and Norwood, J. D., "Framework for a Laser Navigation Sys- 
tem", RICIS ’88 Symposium, Nov. 7-9, 1989. 

Regalbuto, M. A., Cheatham, J. B. and Krouskop, T. A., "A Model-Based Graphics Interface for 
Controlling a Semi- Autonomous Mobile Robot", 11th Annual Conference EEEE-EMBS ’89, 
Seattle, Nov. 9-12, 1989. 

Cheatham, J. B. and Adnan, S., "Kinematic Analysis and Trajectory Control of a Mobile Omni- 
Directional Robot", First National Conference on Applied Mechanisms and Robotics, Cincin- 
nati, Nov. 5-8, 1989. 

Wu, C. K., Weiland, P. L., and Cheatham, J. B., "A Computer Graphics Testbed for Developing 
and Testing Laser Imaging Algorithms," SPIE/SPSE Symposium on Electronic Imaging, Santa 
Clara, CA., Feb. 11-16, 1990. 


Krouskop, T. A., Cheatham, J. B., Kachroo, P., and Barry, P. A., "Non-Invasive Measurement of 
the Stiffness of Soft Tissue Using an Ultrasonic Perturbator," 8th Annual Conference on 
Biomedical Engineering Research in Houston, Feb. 15-16, 1990. 

Norwood, J., Weiland, P. L., and Cheatham, J. B., "Robot Navigation From Local Sensory Data 
Using Neural Networks," 5th Int. Service Robot Congress, Detroit, June 6, 1990. 

Adnan, S., and Cheatham, J. B., "An Omnidirectional Platform to Simulate a Free-Flying 
Robot," ISMCR ’90 - 1st Int. Sym. Measurement and Control in Robotics, NASA-JSC, June 21, 
1990. 

Chen, Y. C., Walker, I.D., and Cheatham, J. B., "Grasping Analysis for a Multifingered Hand," 
Proceedings International Symposium on Measurement and Control in Robotics, NASA/JSC, 
June 1990, Vol 3, pp J2-3- 1-6. 

Kachroo, P., Krouskop, T. A., Kachroo, A., Cheatham, J. B., and Barry, P., "Ultrasonic Tech- 
nique and Artificial Intelligence: Differentiation of Tissue Types," Annual Fall Meeting of the 
Biomedical Engineering Society, Blacksburg, VA, Oct 23, 1990. 

Regalbuto, M. A., Cheatham, J. B., and Krouskop, T. A., "A Framework for a Practical Mobile 
Robotic Aid for the Severely Physically Disabled," RESNA 13th Annual Conference, Washing- 
ton, D.C., 1990. 

Walker, I. D., Cheatham, J. B. and Chen, Y. C., "An Efficient Method for Computing the Force 
Distribution of a Three-Fingered Grasp," Proceedings SPIE Conference on Cooperative Intelli- 
gent Robotics in Space, Boston, MA, November 1990. 

Wu, C. K., Lin, Y. H., and Cheatham, J. B., "Computer Graphics Modelling for Simulating and 
Testing Robot Vision Systems," Journal of Modelling and Simulation, 10 (2), 67-70, 1990. 

Chen, Y. C., Walker, I. D., and Cheatham, J. B., "A New Approach to Force Distribution and 
Planning for Multifingered Grasp of Planar Objects," submitted to Journal of Robotic Systems, 
1990 

Walker, I. D., Cheatham, J. B., and Chen, Y. C., "An Intelligent Grasp Planning Strategy for 
Robotic Hands," Proceedings SPIE International Symposium on Optical Engineering and Pho- 
tonics in Aerospace Engineering, Orlando, FL, April 1991. 


- 1 - 



Chen, Y. C„ Walker, I. D., and Cheatham, J. B., "A New Approach to Force Distribution and 
Planning for Multifingered Grasps of Solid Objects," Proceedings of the IEEE International 
Conference on Robotics and Automation, Sacramento, CA, April 1991. 

Chen, Y. C., Walker, I. D., and Cheatham, J. B., "A New Vectorized Approach to Visualization 
and Full Use of Kinematic Redundancy," submitted to 1991 IEEE International Conference on 
Robotics and Automation. 

Chen, Y. C., Walker, I. D., and Cheatham, J. B., "Grasping Technology for Dextrous Remote 
Manipulation," submitted to 1DEEA One - The First International Design for Extreme Environ- 
ments Assembly, University of Houston, November 12-15, 1991. 

Adnan, S., Cheatham, J., Galicki, P., and Brock, J. "Telepresence for Remote Robotics Opera- 
tions," submitted to IDEEA One - The First International Design for Extreme Environments 
Assembly, University of Houston, November 12-15, 1991. 

Oldham, T„ Stokes, S., Venverlok, D., and Zearfoss, J. "Wheel for Omnidirectional Platform," 
Final Report, Senior Design Project, Mech 408, Rice University, April 26, 1990 

Fessler, J., Jenkins, E„ Miller, R., Peiffer, T., and Ragan, E., "Seven Degree of Freedom Robotic 
Arm," Final Report, Senior Design Project, Mech 408, Rice University, April 20, 1990. 

Ducceschi, L., Matthews, M., O’Connell, S., and Yuan, L. "The Rice Four-Fingered Dextrous 
Robotic Hand," Final Report, Senior Design Project, Mech 408, Rice University, April 19, 1991. 

Bartosh, B., Boncimino, G., Caldinell, T., Kirkpatrick, D., and Shampine, R., "Four Degree of 
Freedom Robotic Wrist," Final Report, Senior Design Project, Mech 408, Rice University, April 
19, 1991. 

Recent Theses Supervised 

Paul B. Fisher. "Development of Sensors and Algorithms for Automating Robotic Grasping," 
M.S. Thesis in Mechanical Engineering, June 1988. 

Sarmad Adnan. "Kinematic Analysis and Trajectory Control of the Rice Omni-directional 
Mobile Robot," M.S. Thesis in Mechanical Engineering, April 1989. 

Peter Weiland. "Use of Laser Scanning Rangefinders for Autonomous Robotic Navigation," 
M.S. Thesis in Mechanical Engineering, May 1989. 

John Norwood. "Robotic Path Planning and Obstacle Avoidance: A Neural Network Approach," 
M.S. Thesis in Mechanical Engineering, May 1989. 

Chris K. Wu, "The Use of Laser Imaging System for Automated Vehicle Guidance and Space 
Servicing Tasks," PhD Thesis in Mechanical Engineering, April 1990. 

Michael A. Regalbuto, "A Semi-Autonomous Mobile Robot/Teleoperator with Applications as 
an Aid for Severely Handicapped People," PhD Thesis in Mechanical Engineering, January 
1990. 

Pushkin Kachroo, "Ultrasonic Technique and Artificial Intelligence: Differentiation of Tissue 
Types," M.S. Thesis in Mechanical Engineering, May 1990. 
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RICIS RESEARCH ACTIVITY AI.2 

Computer Graphics Testbed to Simulate and 
Test Vision Systems for Space Applications 

by 

J. B. Cheatham 

The major objective of this research activity has shifted from 
computer graphics and vision systems to the broader scope of applying 
concepts of artificial intelligence to robotics. Specificially, the 
research is directed toward developing Artificial Neural Networks, 
Expert Systems and Laser Imaging Techniques for Autonomous Space Robots. 

This activity is being conducted by Wu, Chen, Norwood, Weiland, 
Schuster and Atkinson, who are Mechanical Engineering graduate students 
at Rice, and it is directed by Professor John Cheatham. Major 
accomplishments have been reported in technical papers and two M.S. and 
one Ph.D. thesis. A computer graphics simulator for laser imaging 
systems has been developed by Chris Wu as part of his Ph.D. research. 
This simulator has been utilized in research by Peter Weiland to provide 
a top-down map of the environment from laser range data. John Norwood 
has applied potential theory for obstacle avoidance and developed a 
network for path planning using the representation of the environment 

provided by Welland's research. 

Weiland and Norwood are extending their research while working 
toward Ph.D. degrees. This effort Is now being directed toward direct 
interpretation of laser Imaging data using an Artificial Neural Network 
and extension of the robot navigation to three dimensions. Work is also 
aimed at using an Artificial Neural Network for control of a 7dof 

redundant manipulator with three 3dof fingers. It is hoped that this 

research will be of some assistance to the NASA robotic activities 

related to the Space Station and to exploration of the moon and Mars. 
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Chris K. Wu. "The Use of Laser Imaging System for Automated Vehicle 
Guidance and Space Servicing Tasks", Ph.D. Thesis in Mechanical 
Engineering, Rice University, January 1990. 
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International Automation Conference, Detroit, Michigan, June 1990. 

C. K. Wu, J. B. Cheatham, Y. H. Lin and T. F. Cleghorn. "Computer 

Graphics Modelling for Simulating and Testing Robot Vision Systems", 
to be published in Journal of Simulation and Modelling . 
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